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EDITORIAL 

El último número de 2013 aparece en el contexto de un gran debate nacional sobre las 
reformas estructurales promovidas por el gobierno de la República que se han presentado 
como un avance para el país y que son de gran trascendencia. En realidad, las reformas 
planteadas, algunas de ellas todavía en discusión en el Congreso de la Unión, sí tienen el 
potencial de transformar la manera de abordar algunos de los temas más relevantes para 
la nación, como el laboral -reforma elaborada desde principios del año-, el educativo -de 
promulgación muy reciente-, la reforma hacendarla -apenas aprobada- y la energética 
-que todavía se analiza en las cámaras legislativas federales. De su Impacto y eficiente 
aplicación dependerá el avance del modelo de desarrollo adoptados en este país. Se 
podrán llamar exitosas si logran que el estándar de vida de los habitantes se Incremente, 
que la sociedad en su conjunto tenga un desarrollo y crecimiento más Igualitario, que se 
abatan los índices de pobreza y se reduzcan de manera significativa los niveles de violen- 
cia que aquejan a amplias secciones del territorio mexicano. En caso contrario, no pasarán 
de ser meros Intentos de progreso con muy buenas Intenciones pero con escasos resulta- 
dos prácticos. Habrá que seguir su desarrollo con atención. 

En este número se Incluyen seis temas de Investigación de gran Interés e Importancia. El 
primero trata sobre el aprovechamiento del recurso maderable tropical nativo en una 
comunidad ubicada en la Reserva de la Biosfera Selva El Ocote, en Chlapas, México. En 
dicho estudio se realizaron encuestas seml-estructuradas a los ejldatarlos y se colectaron 
ejemplares para la Identificación de las especies. Se determinó un total de 35 especies 
maderables nativas, Incluidas en 31 géneros y 24 familias. Entre los principales usos que 
la comunidad le da a las especies maderables está extracción de leña y postes para cer- 
cado, la elaboración de herramientas de trabajo y la elaboración de muebles, construcción 
y reparación de casas, todas ellas de Importancia en la vida de la comunidad. 

El segundo trabajo describe la respuesta de Cedrela odorata L. a diversos Inoculantes 
mlcorrízlcos procedentes de dos ecosistemas tropicales. Aquí se realizó un bloensayo con 
diferentes fuentes de Inoculo procedentes de dos ecosistemas tropicales donde se distri- 
buye esta especie en el estado de Veracruz. Las variables evaluadas fueron: blomasa 
seca aérea y radicular, altura, diámetro, número de folíolos, número de esporas y porcen- 
taje de colonización. Los resultados mostraron diferencia estadística en las variables de 
crecimiento de los Inoculantes provenientes de la selva mediana subperennlfolla (sms). 
Aquí se demuestra la fuerte relación micorrízica entre los hongos y C. odorata; ante ello, 
se sugiere reproducir y utilizarlos en la producción de plántulas a fin de lograr una mejor 
adaptación al establecerlas en condiciones naturales. 

El tercer artículo presenta la relación de variables físicas en la productividad de Pinus ari- 
zonica y Pinus engelmannii en el sur de Chihuahua, Méx. Se utilizó el índice de sitio (is) 
para Identificar la Influencia de variables de suelo, clima y relieve sobre el crecimiento de 
las dos especies. Con los modelos de is se Identificó la calidad de estación de cada par- 
cela muestreada. Se realizó un análisis de discriminante canónico (adc) para Identificar 
diferencias multlvarladas entre calidades de estación y un análisis de varianza para Iden- 
tificar la Influencia de las variables físicas en el crecimiento de las especies. El método 
empleado predijo con precisión el crecimiento de ambas especies. El adc mostró que 
existen diferencias multlvarladas entre las calidades de estación, con alta Influencia de la 
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pendiente, disección vertical, profundidad de suelo y temperaturas. De acuerdo con el 
análisis de varianza, los valores para las variables de mayor influencia sobre la mejor cali- 
dad de estación fueron profundidad de suelo, pendiente, temperatura mínima y tempera- 
tura máxima promedio. Estos valores pueden ser de gran utilidad para construir modelos 
espaciales de productividad en los estados de Chihuahua y Durango, México. 

El cuarto es un trabajo sobre el establecimiento de Pinus pseudostrobus Lindl. en tres con- 
diciones de propagación y trasplante en cárcavas desprovistas de vegetación en una locali- 
dad del estado de Michoacán, México. La primera condición corresponde a plantas de 
procedencia conocida, germinadas y propagadas en condiciones controladas, con y sin ino- 
culo con el hongo ectomicorrízico Pisolithus tinctorius (Pers.) Coker et Couch. Como una 
segunda condición, ambos lotes fueron trasplantados con adición de sustrato fértil (fibra de 
coco, corteza de árbol y agrolita) en las cepas. La tercera condición fue una reforestación 
convencional con plantas de la misma especie. Los resultados muestran que los pinos que 
fueron plantados bajo las primeras condiciones presentaron mayor crecimiento y supervi- 
vencia que los pinos de la segunda reforestación. La supervivencia de los pinos inoculados 
fue la más alta, seguida por los de la misma plantación sin inoculo y los de la segunda refo- 
restación. La inoculación se plantea como una estrategia útil en proyectos de restauración 
de sitios severamente degradados y erosionados, con buenas expectativas de incrementar 
el éxito en los programas de reforestación. 

El quinto manuscrito se enfoca a tableros de partículas fabricados con residuos industria- 
les de la madera de Pinus patula. Los objetivos son manufacturar tableros aglomerados de 
partículas homogéneos a dos niveles de densidad, determinar sus propiedades físico- 
mecánicas y evaluar los resultados bajo las normas alemana y venezolana. La absorción 
de agua promedio de los tableros a ambas densidades cumple con lo estipulado en las 
normas. La variación del espesor (hinchamiento) promedio por inmersión en agua para los 
tableros con mayor densidad cumple por lo exigido en ambas normas. Respecto de las 
propiedades mecánicas, tanto el módulo de ruptura en flexión estática como la tensión 
perpendicular (adhesión interna) para ambos niveles de densidad superan el valor mínimo 
aceptado en las normas, demostrando que los residuos industriales de la madera de P 
patula presentan buenas características tecnológicas para la fabricación de tableros de 
partículas de densidad media. 

Finalmente, el sexto trabajo plantea las opciones reales para determinare I turno óptimo en 
sistemas silvo-pastoriles a través de un análisis de inversión. El trabajo propone analizar y 
valorar desde la perspectiva económica-financiera el momento óptimo de cosecha en un 
sistema silvo-pastoril utilizando las Opciones Reales (or) como herramienta para la toma 
de decisiones estratégicas. Esta se basa en un proceso estocástico para la proyección de 
precios mientras que el Valor Actual Neto (van) utiliza una tasa fija en todo el proyecto. Los 
resultados mediante or indican cuando es rentable talar para obtener el mayor beneficio 
desde el punto de vista económico. Estos trabajos ilustran de manera muy adecuada 
diversos problemas a los que se enfrentan quienes hacen el manejo y aprovechamiento de 
recursos forestales en Latinoamérica y ofrecen interesantes y útiles soluciones a dichos 
problemas. Estemos seguros que serán de valor para los lectores. 

Raymundo Dávalos Sotelo 

Editor 


EDITORIAL 

The last issue of the year 2013 will come to light in the context of a great natío nal debate 
on structural reforms promoted by the Government of the Republic. These reforms have 
been presented as a breakthrough for the country, and which are of a great importance. In 
reality such reforms, some of them still under discussion in the Congress of the Union, 
have indeed the potential to transform the way of addressing some of the most relevant 
issues for the nation, as is the case with the labor issue, reform elaborated since the begin- 
ning of the year, the educational theme, of very recent enactment, the taxes on revenue, 
just adopted and, the reform about energy, which still is reviewed and under analysis in the 
federal legislative chambers. The progress ofthe development model adopted in this coun- 
try will depend largely on their impact and efficient implementation. We may cali them 
successful ifthey manage to increase the standard ofliving ofthe inhabitants, ifthe society 
as a whole has more egalitarian growth, ifpoverty rafes are diminished and ifthe levels of 
violence that plague broad sections of the Mexican territory are reduced significantly. 
Otherwise, they will not pass from being mere attempts to progress with very good inten- 
tions but with few practical results. It will be necessary to follow its development closely. 

Six themes of research of great interest and importance are inciuded in this number. The 
first one deais with the use ofthe native tropical timber resource in a community located on 
the reserve of the biosphere Selva El Ocote in Chiapas, México. In this study, semi-struc- 
tured surveys were conducted with the ejidatarios and specimens were collected for the 
identification ofthe species. A total of 35 native timber species were identified, including 31 
genera and 24 families. Among the main uses that the community gives to the timber 
resources is extraction offirewood and fence poles, the manufacturing ofworking tools and 
the production of furniture, construction and house-repairs, all of which are of importance 
in the Ufe of the community 

The second work describes the response of Cedrela odorata L. to different mycorrhizal 
inoculators from two tropical ecosystems. Here, a bioassay was conducted with different 
sources of inoculum from two tropical ecosystems where this species is distributed in the 
State of Veracruz. The variables evaluated were aerial and root dry biomass, height, and 
diameter, number of leaflets, spore number and percentage of colonization. The results 
showed statistical difference in the growth variables of inoculants from médium rainforest. 
In this study, the strong relationship between fungí mycorrhizal and C. odorata is shown; in 
view of that, it is suggested to reproduce and use them in the production of seedlings in 
order to achieve a better adaptation to the establishment of trees in natural conditions. 

The third article presents the relation of physical variables in the productivity of Pinus ari- 
zonica and Pinus engelmannii in the south of Chihuahua, Mex. The site Índex (si) i/i/as used 
to identify the influence of variables of soil, climate and surface reliefon the growth oftwo 
species. Is models identified the quality of each sampled plot station. A canonical discrimi- 
nant analysis (cda) u/as performed to identify multivariate differences between qualities of 
station, and an analysis of vari anee to identify the influence of the physical variables in the 
growth ofthe species. The method predicted precisely the growth ofboth species. The cda 
showed multivariate differences between the qualities of station, with high influence of the 
slope, vertical dissection, depth of soil and temperatures. According to the analysis of 
vari anee, the valúes for the variables of greater influence on the best quality of station were 
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depth of soil, slope, minimum temperature and máximum average temperature. These 
valúes can be useful for building spatial models of productivity in the States of Chihuahua 
and Durango, México. 

The fourth manuscript is a work on the establishment of Pinus pseudostrobus Lindl. in 
three conditions of propagation and transplant, on gullies devoid of vegetation in a town in 
the state of Michoacán, México. The first condition corresponds to plants of known prove- 
nance, germinated and propagated in control led conditions, with and without inoculum with 
the ectomycorrhizic fungus Pisolithus tinctorius (Pers.) Coker et Couch. As a second con- 
dition, both batches were transplanted with addition offertile substrate (coconut fibre, bark, 
tree and agrolite) strains. A third condition i/i/as a conventional reforestation with plants of 
the same species. The authors evaluated the height, diameter at the height of the base, 
coverage, numberof branches and percentage ofsurvival. The results showthat pine trees 
that were pianted under the first conditions showed better growth and survival than the 
second reforestation pines. The survival ofthe inoculated pine ivas the highest, followed by 
the same plantation without inoculum and the second reforestation. Inoculation is presen- 
ted as a useful strategy in projects of restoration of sites severely degraded and eroded, 
with good expectations of increasing success in reforestation programs. 

The fifth manuscript f ocuses on particleboards made of industrial i /vaste from wood of Pinus 
patula. This study, conducted in Venezuela and Costa Rica, is a proposal to use industrial 
waste ofsawmilTs that use of Pinus patula, for the production ofchipboard. The objectives 
are manufacturing agglomerate boards with homogeneous partióles at two levels of den- 
sity, to determine their physical-mechanical properties and to evalúate the results under 
Germán and Venezuelan standards. Average water absorption for boards at both densities 
compiies with provisions ofthe rules. The average variations ofthe thickness (swelling) by 
immersion in water for boards with higher density, comply with requirements by both stan- 
dards. Mechanical properties, with regard to both the modulus of rupture in static bending 
and tensión perpendicular to the grain (internal bond) for both levels of density, exceed the 
minimum valué accepted by the standards, demonstrating that industrial waste from wood 
ofP. patula has good technological characteristics for médium density partióle board manu- 
facturing. 

Finally, the sixth paper presents the real options to determine the optimal turn in silvo- 
pastoral systems through an analysis of investment. The work proposes to analyze and 
assess from the perspective of economic-financial the optimum time of harvest in a silvo- 
pastoral system using real options (ro) as a tool for strategic decision-making. It is based 
on a stochastic price projection process while the net present valué (npv) uses a fixed rate 
throughout the project. Results through ro indicate when it is profitable to cut to get the 
most benefit from the economic point of view. 

These works il lústrate very properly various problems facing those who carry out the mana- 
gement and exploitation of forest resources in Latin America and provide interesting and 
useful Solutions to these problems. We are sure that they will be of valué to the readers of 
the Journal. 


Raymundo Davalos-Sotelo 

Editor 


Madera y Bosques 19(1), 2013:7-21 
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ARTÍCULO DE INVESTIGACIÓN 


Aprovechamiento del recurso maderable 
tropical nativo en la comunidad 
de Emilio Rabasa, Reserva de la Biosfera 
Selva El Ocote, Chiapas, México 

Native tropical timber resource use in the community 
of Emilio Rabasa, at the Selva El Ocote Biosphere Reserve 

in Chiapas, México 

Carolina Orantes-García , 1 Miguel Ángel Pérez-Farrera , 1 
Carlos Uriel del Carpio-Penagos 2 y César Tejeda-Cruz 3 

RESUMEN 

El presente estudio tuvo como finalidad conocer el aprovechamiento de los recursos maderables tropi- 
cales nativos por parte de la comunidad Emilio Rabasa, ubicada en la Reserva de la Biosfera Selva El Ocote, 
en el estado de Chiapas, México. Se realizaron encuestas sem ¡estructuradas a 30 ejidatarios y se colectaron 
ejemplares para la identificación de las especies. Se determinó un total de 35 especies maderables nativas, 
incluidas en 31 géneros y 24 familias. Las principales especies maderables multipropósito mencionadas fueron 
el bojón (Cordia alliodora (R. & P.) Oken), el palo amarillo (Terminalia obovata (R. & P.) Steud) la caoba (Swie- 
tenia macrophylla G. King), el jolocín ( Hellocarpus donnell-smithii Rose), el cedro ( Cedrela odorata M. Roem. 
King) y el copalchi ( Croton guatemalensis Lotsy). Entre los principales usos que la comunidad le da a las 
especies maderables está, 38% la extracción de leña y postes para cercado, 17% para la elaboración de 
herramientas de trabajo y 16% para la elaboración de muebles, construcción y reparación de casas. 

PALABRAS CLAVE 

Árboles multipropósito, flora nativa, recurso forestal, selva Zoque, trópico de México. 


ABSTRACT 

This study aimed to identify the use of native tropical timber resources in the Emilio Rabasa commu- 
nity located in the Ocote Forest Biosphere Reserve, in Chiapas State, México. 30 semi-structured surveys 
were conducted with ejidatarios “land-owners”, and several specimens were collected fortheir identification. 
A total of 35 native timber species were determined, included in 31 genera and 24 families.The main multi- 
purpose timber species mentioned were bojón (Cordia alliodora (R. & P.) Oken), palo amarillo ( Terminalia 
obovata (R. & P.) Steud), mahogany ( Swietenia macrophylla G. King), the jolocín ( Heliocarpus donnell- 
smithii Rose), cedar ( Cedrela odorata M. Roem. King) and copalchi ( Croton guatemalensis Lotsy). Among 
the main uses that the community gives these timber species are 38% for extraction of firewood and fence 
posts, 17% for making tools and 16% for making furniture, construction and houses repair. 

KEY WORDS 

Multipurpose trees, native flora, resources, Zoque forest, tropic of México. 


1 Universidad de Ciencias y Artes de Chiapas (UNICACH)-Facultad de Ciencias Biológicas, Chiapas, México. 
C.e.: c_orantes@hotmail.com 

UNICACH-Centro de Estudios Superiores de México y Centroamérica. 

UNICACH-Facultad de Ingenierías. 
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Aprovechamiento del recurso maderable tropical nativo 


INTRODUCCIÓN 

México es uno de los países del mundo 
con mayor diversidad biológica y cultural; 
se ubica entre los cinco primeros países 
llamados “megadiversos”, los cuales 
albergan entre 60% y 70% de la diversi- 
dad biológica conocida del planeta (Mitter- 
meier et al., 2000; Francis y Currie 2003; 
Conabio, 2006). Posee 10% de la flora del 
mundo, con un estimativo de 30 000 espe- 
cies de plantas vasculares, ocupando el 
cuarto lugar a nivel mundial en riqueza de 
plantas endémicas y nativas (Villaseñor y 
Espinosa-García, 2004; Conabio, 2006). 
De acuerdo con Rzedowski (1992), en el 
territorio mexicano se encuentran alrede- 
dor de 2500 especies de árboles nativos, 
que constituyen un reservorio de recursos 
bióticos de valor actual y potencial por los 
múltiples productos y beneficios ambien- 
tales que de ellos se derivan, motivo por el 
cual su presencia y permanencia es fun- 
damental para contrarrestar los efectos 
de la deforestación así como para garanti- 
zar la sustentabilidad de la tierra (Gligo, 
1993; Vázquez et al., 1999; Mardocheo et 
al., 2001). 

Aunque en la actualidad la gestión 
sostenible de los bosques ha sido amplia- 
mente adoptada como una política y un 
objetivo de la administración pública, las 
altas tasas de pérdida y degradación del 
bosque se siguen produciendo en muchas 
áreas. Las actividades humanas como la 
tala, la ganadería, los cultivos, el desarro- 
llo de infraestructura y el uso del fuego, 
están causando la pérdida generalizada 
de biodiversidad, restringiendo el pro- 
greso del manejo sostenible de los bos- 
ques. Sólo en México, la tasa anual de 
deforestación (237 000 ha) es la causa 
principal de la pérdida de vegetación en 
los ecosistemas tropicales (Viana et al., 
1997; Witmore, 1997; Dirzo, 2001). Ante 
tal problemática, existe una urgente nece- 
sidad de herramientas que logren propor- 
cionar una evaluación integrada de los 


impactos humanos sobre la diversidad 
biológica forestal y que puedan apoyar la 
toma de decisiones relacionadas con el 
uso del bosque (Newton et al., 2009). 

La Reserva de la Biosfera Selva El 
Ocote (Rebiso) en el estado de Chiapas 
(estado cuya riqueza arbórea es de 1400 
a 1500 taxa; González-Espinosa et al., 
2004), con una superficie total de 101 288 
- 15 - 12,5 ha (Semarnat, 2000), es consi- 
derada de especial importancia biológica 
y cultural por su cercanía con las selvas 
de los Chimalapas en Oaxaca y Uxpa- 
napa en Veracruz, que en su conjunto for- 
man un corredor biológico que constituye 
la selva zoque (Gómez-Pompa y Dirzo, 
1995). La Rebiso protege uno de los cen- 
tros de diversidad biológica más importan- 
tes de México y del mundo, ya que se 
encuentra ubicada en una zona de transi- 
ción de dos provincias neotropicales, la 
Pacifiqúense y la Tehuantepecana (Vás- 
quez y March, 1996). Además, es un 
espacio donde se satisfacen necesidades 
básicas de comunidades campesinas 
locales, para lo cual hacen uso de los 
recursos forestales nativos, lo que puede 
ocasionar un deterioro de dichos recursos 
(Moreno, 2009). En la Rebiso están pre- 
sentes numerosas especies forestales de 
importancia económica (maderable, medi- 
cinal, comestible, forrajera y ornamental), 
no obstante, hoy en día no se cuenta con 
estudios exhaustivos al respecto. Las con- 
tribuciones al conocimiento de dicha flora 
son todavía incipientes, pocos son los 
autores (Miranda, 1952; Vásquez-Sán- 
chez, 1988 y Ochoa-Gaona, 1996) que 
han llevado a cabo algunos trabajos, prin- 
cipalmente estudios florísticos. El pre- 
sente trabajo es el primer estudio que 
registra el aprovechamiento del recurso 
maderable nativo dentro de la Rebiso, 
como un primer paso para establecer pro- 
gramas y evaluaciones integrales del 
manejo de estas especies por los habitan- 
tes de las comunidades asentadas en el 
área de amortiguamiento de la misma. 
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OBJETIVO 

El presente estudio tuvo como finalidad 
conocer el aprovechamiento de los recur- 
sos maderables tropicales nativos por 
parte de la comunidad de Emilio Rabasa, 
localizada en la Reserva de la Biosfera 
Selva El Ocote, en el estado de Chiapas, 
México. 


MATERIALES Y MÉTODO 
Área de estudio 

La Rebiso se ubica dentro de los munici- 
pios de Cintalapa, Ocozocoautla, Tecpa- 
tán y Jiquipilas, en el estado de Chiapas, 
México, entre los paralelos 16°45'42” y 
17°09'00” de latitud norte y entre los meri- 
dianos 93°21'20” y 93°47'00” de longitud 
oeste (Semarnat, 2000). Dentro de la 
zona de amortiguamiento de la Rebiso se 
encuentra la localidad de estudio, Emilio 
Rabasa, la cual pertenece al municipio de 
Ocozocoautla (Fig. 1). El paisaje predomi- 
nante en la zona corresponde al de las 
montañas del norte, constituidas por sie- 
rras y serranías de altitud variable entre 
los 750 msnm y 1500 msnm. De acuerdo 
con la clasificación de Kóppen, modificada 
por García (1980), el grupo climático pre- 
sente en la zona corresponde al cálido 
subhúmedo con lluvias abundantes en 
verano (Aw^ y con una precipitación total 
anual que varía entre 1500 mm y 2500 
mm, con una temperatura media anual de 
22 °C. 

Según Rzedowski (1994), la vegeta- 
ción presente en la Rebiso corresponde a 
cuatro asociaciones vegetales: 

a) Selva alta perennifolia, predomi- 
nan las especies de chicozapote ( Manil - 
kara zapota (L.) Van Royen), caoba 
( Swietenia macrophylla G. King), cedro 
( Cedrela odorata M. Roem. King) y mojú 
(Brosimum alicastrum Sw.). 


b) Selva alta o mediana subperenni- 
folia, se encuentran especies como el 
molinillo (Quararibea gentlei Lundell), el 
copalillo ( Bursera bipinnata (Sessé & 
Mociño) Engl.), el palo mulato ( Bursera 
simaruba (L.) Sarg.) y el frijolillo ( Pithece - 
llobium arboreum (L.) Urban). 

c) Selva mediana o baja perennifolia, 
las especies representativas son: palo de 
coleto ( Oreopanax peltatus Linden ex 
Regel) y memelita ( Clusia flava Jacq.). 

d) Vegetación secundaria, en donde 
las especies representativas son el corcho 
colorado ( Belotia mexicana (DC.) Schum.), 
la yerba buena de monte ( Pilea pubes- 
cens Liebm.), el guarumbo ( Cecropia pel- 
tata L.) y el macús ( Calathea allouia 
Aubl.). 

De acuerdo con la información pro- 
porcionada por los habitantes de la comu- 
nidad de Emilio Rabasa, los primeros 
pobladores llegaron en 1970 en busca de 
terrenos para la agricultura, provenientes 
de los municipios de Villaflores, Ocozo- 
coautla, Jiquipilas y Tecpatán, Chiapas. En 
1981 se fundó el ejido Emilio Rabasa a tra- 
vés de un decreto federal, con una exten- 
sión de 840 ha. La comunidad está 
integrada por 150 habitantes 82 mujeres y 
68 hombres. Los habitantes son mestizos, 
hablan castellano y 70% saben leer y escri- 
bir. Noventa por ciento de los habitantes de 
la comunidad son ejidatarios que tienen en 
promedio de 20 ha a 30 ha de terreno cada 
uno, mientras que el 10% restante son per- 
sonas que viven en la comunidad que sólo 
poseen una casa y un solar. 

Dentro de las actividades económi- 
cas de la comunidad se encuentra la agri- 
cultura, principalmente la siembra de maíz 
(Zea mays L.) y frijol ( Phaseolus vulgaris 
L.), que es practicada por 37% de los pro- 
ductores; 32% se dedica a la producción 
de café (Coffea arabica L.), 22% a la api- 
cultura y únicamente 9% practica la gana- 


MACROLOCALIZACIÓN 


10 


Aprovechamiento del recurso maderable tropical nativo 



0 

£ * 
<D 00 

> CD 
CO T- 


CT3 


o D 

o E 
.Sí, 5 
2 n> 

O CL D 


0000061. 0000681. 0000881. 0000/81. 000098L 


~ I I I I 

0000061- 0000681 0000881- 0000/81- 0000981 


O 

m 


I ¿-0 <& ^ 


« 


Figura 1. Comunidad de Emilio Rabasa en la zona de amortiguamiento Reserva de la Biosfera Selva el Ocote, Ocozocoautla, Chiapas 

(Semarnat, 2000). 
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dería ovina; tales actividades son para 
autoconsumo. La cantidad de terreno des- 
tinada a la producción por cada productor 
oscila entre las 2 ha y 3 ha con un rendi- 
miento de 800 kg/ha para el maíz y 500 
kg/ha para el frijol. El cultivo del café 
(Coffea arabica) es la segunda actividad 
productiva comercial en importancia den- 
tro de la comunidad por la superficie ocu- 
pada, el rendimiento promedio de este 
cultivo es de 700 kg/ha a 1000 kg/ha de 
café convencional. El rendimiento ha 
aumentado con la utilización de técnicas 
orgánicas, la comercialización la realizan 
con introductores externos. 

La apicultura constituye una alterna- 
tiva productiva complementaria a la agricul- 
tura, al mismo tiempo que es una actividad 
compatible con la conservación de los 
recursos forestales de la comunidad. El 
producto principal es la miel, que es ven- 
dida a productores externos de la cabecera 
municipal de Ocozocoautla, Chiapas. 
Actualmente los productores están en un 
proceso de capacitación y ampliación. 


Método 

Se aplicaron encuestas semiestructuradas 
para la recopilación de la información (Geil- 
fus, 1997). Las encuestas se realizaron a 
los 30 ejidatarios de Emilio Rabasa, con 
edades entre los 20 y 70 años. Esta herra- 
mienta permitió indagar con respecto a las 
formas de uso del recurso forestal madera- 
ble, especies utilizadas, técnicas emplea- 
das y percepción. El trabajo se acompañó 
de la observación directa así como de con- 
versaciones y cuatro recorridos por las 
zonas de aprovechamiento junto con los 
informantes. Los datos recopilados se 
registraron en una hoja de cálculo de 
Microsoft Excel para facilitar su análisis, 
mediante una estadística descriptiva. 

De acuerdo con la clasificación de 
Rzedowski (1994), se determinó el tipo 


de vegetación donde se encontraron las 
especies maderables nativas. Siguiendo 
la metodología de Lot y Chiang (1986), se 
realizaron colectas de ejemplares de las 
especies arbóreas nativas para su poste- 
rior identificación, con ayuda de bibliogra- 
fía especializada (Martínez, 1979; 
Breedlove, 1981; Niembro, 1986; Rze- 
dowski, 1992; Gentry, 1996; Miranda, 
1998; Ricker y Douglas, 1998; Penning- 
ton y Sarukhán, 2005; Flora Mesoameri- 
cana, 2006), los ejemplares identificados 
fueron depositados en el Herbario Eizi 
Matuda (hem) de la Universidad de Cien- 
cias y Artes de Chiapas. 


RESULTADOS 

Historia de aprovechamiento 

De acuerdo con las encuestas, el aprove- 
chamiento forestal dio inicio en 1970 
cuando los primeros pobladores prove- 
nientes de los municipios de Villaflores, 
Ocozocoautla, Jiquipilas y Tecpatán, 
Chiapas se establecieron en los terrenos 
selváticos en los que actualmente se 
ubica la comunidad. Estos primeros pobla- 
dores talaban y quemaban las especies 
maderables nativas para abrir campos de 
siembra de cultivos como el maíz (Z. 
mays .) y frijol (P. vulgaris ); aunado a ello, 
la madera sólo se destinaba para la cons- 
trucción de las casas y no se comerciali- 
zaba debido a que no existían caminos 
para transportarla. 

A inicios de 1975, el aprovechamiento 
forestal maderable se intensificó con la 
apertura de la carretera y el estableci- 
miento de un aserradero privado dentro de 
los terrenos de la comunidad. A partir de 
ello, los habitantes de la comunidad pudie- 
ron vender los árboles de maderas precio- 
sas como el cedro ( Cedrela odorata) y la 
caoba ( Swietenia macrophylla) a la 
empresa maderera. En 1980, después de 
haber cerrado el aserradero, los habitantes 
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de la comunidad comenzaron a explotar la 
madera como principal actividad, misma 
que presentó una alternativa de genera- 
ción de recursos económicos inmediatos. 
Dicha madera era comprada, principal- 
mente por introductores provenientes de la 
cabecera municipal de Ocozocoautla. 
Actualmente, la explotación de la madera 
está regulada por la misma comunidad a 
través de un reglamento interno y por las 
instituciones vinculadas con este sector; la 
madera es utilizada en la mayoría de los 
casos para reparación de las casas, leña y 


cercado; existen programas de reforesta- 
ción en zonas degradadas. 


Aprovechamiento forestal 

A través de las encuestas se registró un 
total de 35 especies maderables nativas, 
distribuidas en 31 géneros y 24 familias 
(Tabla 1). Ochenta por ciento (n=35) se 
encuentran distribuidas en selva alta 
perennifolia y el 20% restante en la selva 
mediana subperennifolia. 


Tabla 1. Número total de familias, géneros y especies forestales maderables nativas 
aprovechadas en el ejido Emilio Rabasa, Rebiso, Chiapas. 


Familias 

Géneros 

Especies 

Amaranthaceae 

1 

1 

Bombacaceae 

2 

2 

Boraginaceae 

1 

1 

Burseraceae 

1 

3 

Cochlospermaceae 

1 

1 

Combretaceae 

1 

1 

Euphorbiaceae 

1 

1 

Fabaceae 

3 

3 

Flacurtiaceae 

1 

1 

Lauraceae 

1 

1 

Meliaceae 

3 

5 

Moraceae 

1 

1 

Myrsinaceae 

1 

1 

Myrtaceae 

1 

1 

Ochnaceae 

1 

1 

Polygonaceae 

1 

1 

Ramnaceae 

1 

1 

Rubiaceae 

2 

2 

Rutaceae 

1 

1 

Sapindaceae 

1 

1 

Sapotaceae 

1 

1 

Simaroubaceae 

1 

1 

Tiliaceae 

1 

1 

Ulmaceae 

2 

2 

TOTAL 24 

31 

35 
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De acuerdo con los encuestados, en 
Emilio Rabasa, 70% (n=35) de las espe- 
cies maderables utilizadas son multipro- 
pósito (Tabla 2). Se registraron 10 usos 
que se les da a las especies: leña, postes, 
cercas vivas, construcción de casas, ela- 
boración de muebles, herramientas de tra- 
bajo, medicinal, comestibles, melíferas y 
forrajeras, en donde 38% (n=35) de las 
especies son utilizadas en la extracción 
de leña y postes para cercado, 17% para 
la elaboración de herramientas de trabajo, 
16% para la elaboración de muebles, 
construcción y reparación de casas, 12% 
son melíferas, 8% son utilizadas en la 
medicina tradicional, 5% como cercos 
vivos y 4% como comestible y forrajera. 

De 43% (n=35) de la especies made- 
rables nativas se utilizan sólo las ramas y 
troncos para la obtención de leña, cercado 
de corrales y herramientas de trabajo, de 
16% se usa la hoja, tal es el caso de la 
colita de borrego ( Picramnia sp.),la 
pimienta ( Pimienta dioica), cuyas hojas 
son utilizadas en la medicina tradicional y 
el mojó (Brosimum alicastrum) que se 
aprovecha como forraje; de 15% de los 
árboles como el bojón ( Cordia alliodora), 
el palo amarillo ( Terminalia obovata), la 
caoba (Swietenia macrophylla), el cedro 
( Cedrela odorata) y el hormiguillo (Platy- 
miscium dimorphandrum) se utiliza todo el 
árbol para diferentes fines; de 11% se 
aprovecha la corteza, como el palo mulato 
(. Bursera simaruba), el copalillo ( Bursera 
bipinnata) y el copal ( Bursera excelsa ); de 
8% de las especies se utiliza únicamente 
la flor, como el sospó (Pseudobombax 
ellipticum) y el pumposhuiti ( Cochlosper - 
mum vitifolium); y de 7% de los árboles 
maderables nativos se aprovechan los fru- 
tos y las semillas, tal es el caso del mojú 
(Brosimum alicastrum), el chicozapote 
( Manilkara zapota) y la pimienta ( Pimienta 
dioica). 

Las especies como el cedro ( Cedrela 
odorata), la caoba ( Swietenia macro- 


phylla), el bojón ( Cordia alliodora), el mojú 
(Brosimum alicastrum), el palo amarillo 
(Terminalia obovata), el jolocín (Heliocar- 
pus donnell-smithii) y el copalchi (Croton 
guatemalensis), son las que presentan 
mayor frecuencia de uso (100%, n=30), 
mientras que especies como el espino 
(Acacia penn atula L.), el frijolillo (Pithe- 
cellobium arboreum (L.) Urban), el 
pumposhuiti (Cochlospermum vitifolium 
Willd. ex Spreng) y el sospó (Pseudobom- 
bax ellipticum (H.B.K.) Dugand) son uti- 
lizadas en menor medida (10%) por los 
habitantes de la comunidad. 

De acuerdo con los encuestados, en 
el ejido Emilio Rabasa, los hombres inte- 
grantes del núcleo familiar realizan el apro- 
vechamiento de las especies maderables 
nativas. La madera es cortada in situ ; para 
la extracción de leña y postes, la madera 
es cortada mediante hacha y llevada a la 
comunidad por las mujeres, mientras que 
para la construcción de viviendas, fabrica- 
ción de muebles y puertas, se utiliza princi- 
palmente la motosierra y es transportada 
mediante arrastre por caballos. 

Los árboles son seleccionados 
según el uso; para leña, postes y herra- 
mientas de trabajo se seleccionan árboles 
de 2 años o más, con d.a.p. (diámetro a la 
altura del pecho) entre 0,10 m y 0,30 m, 
altura de 2 m a 5 m; para construcción de 
casas y muebles se seleccionan árboles 
de 5 años en adelante, con d.a.p. de 0,50 
m a 1 ,50 m y altura de 1 5 m a 40 m. 

Los encuestados también indicaron 
que en la comunidad la venta de madera 
es casi nula, ya que dicha actividad está 
regulada por la misma comunidad a través 
de un reglamento interno para la explo- 
tación forestal y por las instituciones vin- 
culadas con este sector (Profepa, Conanp 
y Semarnap). Sin embargo, también 
aclararon que aún existen pobladores que 
venden árboles en pie, con precios que 
van de MXP$1000 a MXP$4000, con 
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compradores de la cabecera municipal de 
Ocozo-coautla, Chiapas. Las principales 
especies arbóreas nativas que se comer- 
cializan se enlistan en la tabla 3. 


DISCUSIÓN 

En Emilio Rabasa, más de 70% de los 
árboles nativos se aprovechan para la 
extracción de leña, elaboración de herra- 
mientas de trabajo, postes para cercado, 
elaboración de muebles, construcción y 
reparación de casas. Estos resultados 
coinciden con Torres (2004), al sugerir 
que el uso de leña para el hogar repre- 
senta, en pleno siglo XXI, 40% de la 
demanda de energía doméstica nacional. 

Las especies arbóreas nativas 
encontradas en Emilio Rabasa como el 
cedro ( Cedrela odorata), la caoba (Swie- 
tenia macrophylla), el bojón ( Cordia al lío- 
dora), el chicozapote ( Manllkara zapota ), 


el mojó ( Brosimum alicastrum ) y el palo 
amarillo ( Termlnalla obovata), son men- 
cionadas en diversos estudios como 
especies multipropósito para reforesta- 
ción, restauración y plantaciones foresta- 
les. Entre estos estudios se pueden 
señalar el realizado en el área de Jurisdic- 
ción, cdmb (Corporación Autónoma Regio- 
nal para la Defensa de la Meseta de 
Bucaramanga), Santander, Colombia, 
(2002) en el que se llevó a cabo una eva- 
luación para determinar el potencial de 
especies forestales nativas para la refo- 
restación, identificándose algunas espe- 
cies forestales nativas promisorias en 
diferentes formaciones ecológicas. Entre 
estas especies se encuentran el frijolito 
( Schlzoloblum parahybum), cedro 
( Cedrela odorata), anime ( Protlum sp.), 
pardillo ( Cordia alllodora), balso ( Ochroma 
pyramldale), caimo ( Pouterla sp.), cedrón 
( Simaba cedrón), chingalé ( Jacaranda 
cópala) y ojite ( Broslmun alicastrum), 
entre otras. 


Tabla 3. Especies arbóreas nativas comercializadas, mencionadas por los habitantes de 

Emilio Rabasa. 


Nombre científico 

Nombres comunes 

Uso 

Cedrela odorata 

Cedro 

Muebles, puertas 

Coccoloba cozumelensis 

Carnero 

Leña, poste 

Cordia alliodora 

Bojón 

Muebles, construcción 

Croton guatemalensis 

Copalchi 

Poste 

Faramea occidentalis 

Huesito 

Leña 

Heliocarpus donnell-smithii 

Jolocín 

Leña 

Swietenia macrophylla 

Caoba 

Muebles, puertas 

Terminalia obovata 

Palo amarillo 

Muebles, construcción 

Zuelania guidonia 

Palo de paragüita 

Poste 
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Por su parte, Rivas (2003) explicó 
que en El Salvador se realizaron impor- 
tantes investigaciones sobre el estableci- 
miento y manejo de plantaciones 
forestales en las que se identificaron 49 
especies forestales con potencial produc- 
tivo; entre las especies que se mencionan 
en dichos trabajos están el almendro de 
río ( Andira inermis), balsa ( Ochroma lago- 
pus), bálsamo ( Miroxylon balsamun), 
barillo ( Calophylum brasiliense), caoba 
( Swietenia macrophylla) carao (Cass/a 
grandls), caulote ( Guazuma ulmifolia), 
cedro ( Cedrela odorata) y pajarito ( Cordia 
alliodora). De ellas, S. macrophylla, C. 
Odorata y C. alliodora, fueron especies 
con mayor número de menciones en la 
comunidad Emilio Rabasa. 

Por su parte, el catie (2006) llevó a 
cabo el proyecto llamado “Developing best 
practices for seed sourcing of planted and 
natural regeneration in theneotropics”, 
conocido como Proyecto Seedsource, en 
la región comprendida desde el sur de 
México hasta Panamá. En dicho proyecto 
se registraron un total de 50 especies 
forestales prioritarias del trópico ameri- 
cano, de las cuales la caoba ( Swietenia 
macrophylla), el cedro (Cedrela odorata), 
el pajarito ( Cordia alliodora) y el copalchi 
( Croton guatemalensls) resultaron tener la 
mayor demanda para fines de reforesta- 
ción comercial. Estas especies coinciden 
con las mencionadas por los habitantes de 
la comunidad aquí estudiada como espe- 
cies de maderables nativas que son comer- 
cializadas a escala local. 

Otros trabajos que concuerdan con 
las especies multipropósito encontradas 
en la comunidad en estudio, son los reali- 
zados por Wightman et al. (2006) en la 
Amazonia Peruana, en donde la produc- 
ción maderera es una de las actividades 
económicas más importantes. Por ello 
desarrollaron un manual en el que listan 
42 especies vegetales de importancia 
maderable, de las cuales se estudiaron 


varios aspectos de cultivo. Entre las espe- 
cies mencionadas están: cedro ( Cedrela 
odorata), amasisa ( Erythrlna ulel), bolaina 
negra ( Guazuma ulmifolia), cacao ( Theo - 
broma cacao), caoba ( Swietenia macro- 
phylla), copalchi ( Croton guatemalensls), 
copaiba ( Copalfera offlclnalls), cumala 
( Virola calophylla) y estoraque ( Myroxylon 
balsamum). 

Igualmente, Palacios et al. (2008) lle- 
varon a cabo un estudio sobre el conoci- 
miento del aprovechamiento de los 
recursos naturales, en las comunidades 
de Chintadó, Doña Josefa y Yuto, en la 
cuenca media del río Atrato, Chocó, en 
Colombia. En dicho estudio encontraron 
que las principales especies maderables 
extraídas son caoba ( Swietenia macro- 
phylla), chanú (Humlrlastrus procenum), 
cedro ( Cedrela odorata) y lirio ( Couma 
macrocarpa). Lo que concuerda con el 
estudio en la comunidad Emilio Rabasa 
en donde se encontró a S. macrophylla y 
Cedrela odorata como especies madera- 
bles nativas extraídas. 

En México, Rzedowski (1992) encon- 
tró la presencia de aproximadamente 
9000 especies de árboles y arbustos nati- 
vos en el país. De éstos, Vázquez y Cer- 
vantes (1993) realizaron una búsqueda 
bibliográfica de aquellos potencialmente 
valiosos para la restauración ecológica y 
la reforestación; durante dicha investiga- 
ción seleccionaron 233 especies leñosas 
nativas de México. El criterio para la 
selección de estas especies se centró 
básicamente en que fueran especies mul- 
tipropósito y con características para ser 
empleadas en programas de restauración 
y reforestación en las diferentes regiones 
ecológicas del país. Entre las 233 espe- 
cies se encuentran Broslmum allcastrum, 
Bursera blplnnata, Bursera slmaruba, 
Cedrela odorata, Coccoloba cozumelen- 
sls, Cordia alliodora, Croton guatemalen- 
sls, Faramea occidentales, Swietenia 
macrophylla, Termlnalla obovata, Manll- 
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kara zapota, Pimienta dioica, Platymis- 
cium dimorphandrum, mencionadas por 
los habitantes de Emilio Rabasa como 
especies multipropósito. 

Posteriormente, Benítez et al. (2004) 
realizaron estudios sobre árboles multiu- 
sos nativos de Veracruz para reforesta- 
ción, restauración y plantaciones, algunas 
de las especies encontradas fueron el 
canshán ( Terminalia obovata), molinillo 
( Quararibea gentlei), chicozapote ( Mani - 
Ikara zapota), caoba ( Swietenia macro- 
phylla), guapaque ( Dialium guianense), 
mojó ( Brosimum alicastrum), cedro rojo 
( Cedrela odorata), ceiba ( Ceiba pentan- 
dra), baqueta ( Ulmus mexicana) y sonza- 
pote ( Licania platypus). Las cuales 
también son reportadas en la comunidad 
de Emilio Rabasa. 

Ahora bien, en la comunidad Emilio 
Rabasa las ramas y los troncos son las 
partes más usadas de los árboles. El prin- 
cipal aprovechamiento de los árboles se 
realiza con el fin de obtener materiales 
para combustible de uso doméstico, cons- 
trucción y postes para cercados. Ello sig- 
nifica que el uso de la madera sigue 
siendo fundamental para la vida diaria de 
las comunidades rurales, proporcionando 
techo, alimento y apoyando las activida- 
des productivas, y explica el por qué las 
ramas y troncos son las partes más usa- 
das. Como lo mencionan Palacios et al. 
(2008), de los árboles nativos se utilizan 
distintas partes vegetativas, desde semi- 
llas, flores, hojas, hasta ramas en porcen- 
tajes que van de 1 0% a 40%. Sin embargo, 
la extracción de madera en donde se 
derriba todo el árbol se realiza en casi 
50% en especies nativas, lo que ha traído 
como consecuencia el deterioro del 
recurso forestal. 

El potencial de producción de bienes 
forestales es considerable, se utilizan 
cientos de productos (hojas, frutos, semi- 
llas, cortezas, gomas, ceras, fibras, colo- 


17 

rantes, entre otros), obtenidos de entre 
5000 a 7000 especies distribuidas en los 
diferentes ecosistemas del territorio nacio- 
nal (Torres, 2004). La falta de datos esta- 
dísticos no permite saber el número 
exacto de productos que se utilizan ni el 
valor de la producción. Sin embargo, 20% 
de la población más marginada del país, 
que coincidentemente habita en las zonas 
rurales, depende para su subsistencia de 
los productos forestales obtenidos de los 
ecosistemas en los que viven (Torres, 
2004). 

En síntesis, existe un conjunto de 
especies maderables nativas amplia- 
mente mencionadas tanto por la diversi- 
dad de usos que tienen, como porque son 
muy demandadas y porque se usan varias 
partes de la planta. A partir de estos datos, 
se sugiere que las especies de Swietenia 
macrophylla, Cedrela odorata, Brosimum 
alicastrum, Cordia alliodora, Manilkara 
zapota, Terminalia obovata, Heliocarpus 
donnell-smithii, Croton guatemalensis, 
Bursera simaruba, Coccoloba cozumelen- 
sis y Cupania glabra deben ser considera- 
das para ser utilizadas en programas de 
restauración, reforestación y plantación 
forestal. Dos especies que no fueron men- 
cionadas por los habitantes de la comuni- 
dad de Emilio Rabasa son Acacia 
pennatula y Pithecellobium arboreum, sin 
embargo, éstas han sido utilizadas en 
otras regiones como especies forrajeras y 
cercos vivos. Incluso, hoy en día, a nivel 
mundial se realizan estudios sobre la con- 
tribución de A. pennatula en la productivi- 
dad agroforestal sostenible (Nieto et al., 
2001), por lo que podrían aprovecharse 
de igual forma dentro de la Rebiso. 

Newton etal. (2009) menciona que la 
falta de información sobre los procesos 
ecológicos y características de las espe- 
cies dificulta el manejo sostenible de los 
bosques. Por ello, es recomendable que 
para realizar un aprovechamiento ade- 
cuado de las especies forestales aquí cita- 
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das, deben realizarse antes estudios 
sobre germinación, crecimiento y el esta- 
blecimiento de las especies, así como la 
evaluación del impacto del aprovecha- 
miento, ya que los cambios en los esque- 
mas de uso de la tierra y los usos y 
costumbres ligados a ello, pueden signifi- 
car la pérdida de especies y de la diversi- 
dad genética, especialmente de la biota 
tropical que fueron seleccionadas en fun- 
ción de los ambientes y microhábitats y 
que hoy perduran en los ecosistemas tro- 
picales (Nasi y Frost, 2009). 


CONCLUSIONES 

Las especies maderables nativas aprove- 
chadas en la comunidad de Emilio Rabasa 
representan 6% de las especies arbóreas 
presentes en el área, 23% del total de las 
familias de Angiospermas encontradas 
para la Reserva de la Biosfera Selva El 
Ocote por Ochoa-Gaona (1996). Dichas 
especies representan una opción para ser 
utilizadas en forma integral de una manera 
económica y ecológicamente viable 
debido a que se presentan especies multi- 
propósito que pueden ser manejadas 
en programas de reforestación, restaura- 
ción y plantaciones, tal como lo sugieren 
Benítez et al. (2004). 

De la lista de 35 especies encontra- 
das en la comunidad de Emilio Rabasa, 
seis coinciden con las especies de mayor 
demanda a nivel nacional e internacional: 
Cedrela odorata, Swletenla macrophylla, 
Cordla alliodora, Manilkara zapota, Brosi- 
mum alicastrum y Terminalia obovata 
(Rivas, 2003; catie, 2006; Wightman et 
al., 2006). 

El presente trabajo muestra que el 
aprovechamiento del recurso forestal 
maderable nativo por parte de la comuni- 
dad sigue siendo vigente para necesida- 
des básicas de la vida diaria, y que hay un 
conjunto de especies preferidas. El uso de 


leña y de postes son dos usos vigentes y, 
en amplios sectores, también la construc- 
ción, de tal modo que se requieren pro- 
puestas de reforestación que garanticen 
el abasto de madera. El uso de postes, 
impulsando las cercas vivas, puede traer 
asimismo beneficios ambientales. Sin 
embargo, es importante la realización de 
un estudio detallado sobre el impacto 
ambiental generado por la extracción 
forestal. 

Las selvas tropicales, como la Selva 
El Ocote, representan un ecosistema 
único con gran biodiversidad y hábitat 
exclusivo de muchas plantas y animales 
en peligro de extinción y son considerados 
como principales fuentes de madera sus- 
ceptibles de ser explotados por las comu- 
nidades locales (Smouts, 2003; 
Mittermeier et al., 2003). 
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ARTÍCULO DE INVESTIGACIÓN 


Respuesta de Cedrela odorata L. 
a diversos inoculantes rmcorrízicos 
procedentes de dos ecosistemas tropicales 

Response of Cedrela odorata L. to several mycorrhizal 
inoculants from two tropical ecosystems 

Heriberto Méndez-Cortés, 1 * José G. Marmolejo-Monsiváis , 2 
César Cantú-Ayala , 2 Víctor Olalde-Portugal , 3 
Eduardo Estrada-Castillón 2 y César Posadas-Leal 1 

RESUMEN 

El presente estudio muestra el efecto de algunos inoculantes de hongos micorrízicos arbuscula- 
res (HMA) sobre el desarrollo inicial de Cedrela odorata L.; para ello, se realizó un bioensayo con 
diferentes fuentes de inoculo procedentes de dos ecosistemas tropicales donde se distribuye esta 
especie en el estado de Veracruz. Se utilizaron cuatro tratamientos: TI, inoculación directa con un 
consorcio de esporas; T2, inoculación con raíces colonizadas de campo; T3, sin inoculación (testigo); 
T4, suelo nativo no estéril (TI, T2 y T3 en suelo nativo esterilizado). Las variables evaluadas fueron: 
biomasa seca aérea y radicular, altura, diámetro, número de foliólos, número de esporas y porcentaje 
de colonización. Los resultados mostraron diferencia estadística en las variables de crecimiento de los 
inoculantes provenientes de la selva mediana subperennifolia (SMS), donde T4 arrojó los mejores 
resultados, seguido por T2, mientras que T3 presentó los valores más bajos. Aunque los inoculantes de 
selva alta perennifolia (SAP) no mostraron diferencias estadísticas, su comportamiento fue similar a 
aquellos de la SMS. Lo anterior muestra la fuerte relación micorrízica entre los HMA y C. odorata ; ante 
elfo, se sugiere reproducir y utilizarlos en la producción de plántulas a fin de lograr una mejor adapta- 
ción al establecerlas en condiciones naturales. 

PALABRAS CLAVE 

Cedro rojo, crecimiento, glomeromycota, selva alta perennifolia, selva mediana subperennifolia. 


ABSTRACT 

This study shows the effect of some inoculants of the arbuscular mycorrhizal fungí (AMF) on the 
initial development Cedrela odorata L; for this, a bioassay was conducted with different inoculum sources 
oftwo tropical forests where this species is distributed in Veracruz State. Fourtreatments were applied: TI, 
direct inoculation with a spore consortium; T2: inoculation with mycorrhizal roots from field; T3: without 
inoculation; T4: native soil not sterile. (TI , T2 and T3 in native sterilized soil). The analyzed variables were 
dry aerial and root biomass, height, diameter, leaflets number, spores number, and colonization percen- 
tage. Statistical analysis showed significant differences between inoculants from médium tropical forest 
(MTF), where T4 was the best results followed by T2, while T3 observed the lowest valúes. Although the 
inoculants of high tropical forest (HTF), do not have statistical significance, the results were similar to those 
in the MTF. These results showed a good micorrhizal relationship between these fungí and C. odorata ; for 
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this reason, it is recommended to repro- 
duce this fungus to use in plants produc- 
tion in order to obtain plants best adapted 
to natural conditions. 

KEY WORDS 

Red cedar, growth, Glomeromycota, high 
tropical forest, médium tropical forest. 

INTRODUCCIÓN 

La alteración de ecosistemas debido a la 
deforestación sigue en muchas ocasiones 
un proceso de erosión muy acelerado, lo 
que ocasiona una baja productividad, 
menor calidad y cantidad de agua, así 
como una alteración en la diversidad bioló- 
gica microbiana (Kennedy y Smith, 1995). 

Cierta microbiota interactúa simbióti- 
camente con la mayoría de las plantas, 
siendo operada por un consorcio o grupo 
de microorganismos capaces de transpor- 
tar y suministrar los nutrimentos necesa- 
rios al interior de sus raíces (Fitter y 
Garbaye, 1994; Bonfante y Anca, 2009). 

Una de las interacciones mutualistas 
más sobresalientes en la naturaleza, la 
forman los hongos micorrízicos arbúscula- 
res (hma) con las plantas, dada su capaci- 
dad de: 1) facilitar la transferencia de 
fósforo y otros elementos importantes 
para el desarrollo celular (Dever et al., 
2001 ; Cuenca et al., 2007); 2) contribuir a 
la formación de la estructura del suelo 
(Piotrowski et al., 2004); 3) permitir una 
mayor tolerancia a periodos prolongados 
de estrés por sequía (Sieverding y Toro, 
1998); 4) actuar como antagonistas contra 
patógenos del suelo (Baath y Hayman, 
1983); 5) mantener la diversidad de plan- 
tas en los diferentes ecosistemas terres- 
tres (Lovelock et al., 2003). Por las 
características anteriormente citadas, se 
considera a esta asociación parte primor- 
dial en los procesos de restauración en 
zonas altamente degradadas (White etal., 
2008), como sucede en las selvas tropica- 
les de México (Mas et al., 2003). 


Se estima que 95% de las especies 
arbóreas que se distribuyen en estos eco- 
sistemas presentan simbiosis con este 
tipo de hongos (LeTacon etal., 1998), los 
cuales son parte esencial para dar funcio- 
nalidad a sus hábitats naturales (Salas, 
2004). 

El cedro rojo (C. odorata L.) es una 
especie de importancia económica que se 
distribuye en la zonas tropicales de 
México, donde sus poblaciones naturales 
han sido seriamente afectadas por los 
impactos antropogénicos (Niembro, 
1996). Ante ello, la Comisión Nacional 
Forestal (Conafor) en las últimas dos 
décadas ha implementado estrategias de 
reforestación y plantaciones forestales 
comerciales para mitigar a este tipo de 
problemas en estas zonas. 

El manejo de los hma, integrado a la 
producción de plántula en vivero, es indis- 
pensable para asegurar el éxito de adap- 
tación de las repoblaciones bajo 
condiciones naturales (Alien et al., 2003); 
por tal motivo, es importante integrar 
paquetes tecnológicos de diferentes espe- 
cies forestales de importancia económica 
que se encuentran distribuidas en los tró- 
picos. 

Los estudios de inoculación a C. odo- 
rata con hma son escasos (Zulueta et al., 
2000; Chable, 2007; Amador, 2010) y se 
han enfocado a utilizar cepas comercia- 
les. Aunque se ha documentado la impor- 
tancia de los hma nativos en la producción 
de planta en vivero (Trejo et al., 2011); en 
México poco se conoce sobre este tema 
en árboles de importancia forestal. 


OBJETIVO 

Determinar el efecto de los inoculantes de 
hma procedentes de dos ecosistemas tro- 
picales del estado de Veracruz sobre el 
desarrollo inicial de C. odorata. 
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MATERIALES Y MÉTODO 

Se muestrearon suelo y raíces en dos 
ecosistemas tropicales con vegetación 
primaria donde se distribuye C. odorata en 
el estado de Veracruz. Se seleccionaron 
nueve árboles por tipo de vegetación en 
los meses de mayo y junio. Según la carta 
de uso de suelo y vegetación serie III del 
inegi, las cartas temáticas de precipita- 
ción, temperatura y elevación de la Cona- 
bio, así como la Conafor (2010), estos 
sitios presentan las siguientes caracterís- 
ticas: el primer ecosistema corresponde a 
una selva alta perennifolia (sap), con 
muestras provenientes del municipio de 
Tezonapa (18°36’ N, 96°41’ O), la altura 
sobre el nivel del mar (asmn) es de 1 80 m, 
se caracteriza por presentar climas tipo 
cálido húmedo, precipitación media anual 
de 2885 mm, temporada seca de tres 
meses como máximo, suelos profundos y 
alto contenido de materia orgánica. El 
segundo comprende una selva mediana 
subperennifolia (sms) con muestras prove- 
nientes del municipio de Tantoyuca 
(21°22’ N, 98°15’ O), se caracteriza por 
presentar árboles con menor altura con 
respecto a la sap, la asmn es de 
100 m-120 m, clima cálido subhúmedo, 
precipitación de 1200 mm-1500 mm anua- 
les, temporada de sequía muy prolongada 
y suelos menos profundos. 

La evaluación de la respuesta de C. 
odorata a la aplicación de inoculantes 
micorrízicos nativos, se realizó mediante 
un bioensayo en invernadero del cinves- 
TAV-Unidad Irapuato, donde se utilizaron 
diferentes fuentes de inoculo procedentes 
de la sap y sms. Por cada ecosistema se 
establecieron cuatro tratamientos bajo un 
diseño completamente al azar: TI) suelo 
nativo estéril y 100 esporas tomadas al 
azar por ecosistema (se identificaron los 
géneros Funneliformis, Glomus, Rhi- 
zophagus, Acaulospora, Diversispora y 
Pacispora para sap; y Claroideoglomus, 
Funneliformis, Glomus, Rhizophagus, 


Sclerocystis, Paraglomus, Entrophospora, 
Scutellospora y Pacispora para sms con 
base en las características propuestas por 
Morton (1 988) y Walker (1983)); T2) suelo 
nativo estéril y cinco gramos de raíces 
como fuente de inoculo por ecosistema 
(con presencia de micorriza superior a 
60% previamente evaluada); T3) suelo 
nativo estéril por ecosistema (testigo); T4) 
suelo nativo sin esterilizar por ecosistema 
(con 8 esporas/g y 2 esporas/g de suelo 
previamente cuantificadas en sap y sms, 
respectivamente). Es importante señalar 
que los inoculantes del tratamiento 1 y 2 
fueron aplicados directamente a la semilla 
en la siembra. El número de repeticiones 
integradas a cada tratamiento estuvo 
constituido por seis plantas. 

La esterilización del suelo se realizó 
en autoclave por tres días a 103 kPa (15 
psi) y 121 °C, durante una hora. La semilla 
utilizada fue de árboles con adecuados 
aspectos fenotípicos (sanidad, altura, 
fuste limpio, etc.) provenientes del inifap 
(campo experimental El Palmar, Tezonapa, 
Veracruz). La siembra se realizó en bolsas 
de polietileno de 25 cm x 18 cm. Cada 15 
días y durante cuatro meses y medio, se 
evaluó la altura total de la planta, diámetro 
al cuello de la base y número de folíolos, 
en donde se realizaron curvas de creci- 
miento con sus respectivos errores están- 
dar. Al final, se determinó la biomasa seca 
aérea y biomasa seca radicular; relación 
peso seco aéreo/peso seco radicular; 
número de esporas (con base en el 
método de tamizado húmedo y 
decantación propuesto por Gerdemann y 
Nicolson (1963) y por el método de centri- 
fugación en gradiente de sacarosa de 
acuerdo con Daniels y Skipper (1982)); 
porcentaje de colonización por hifas, vesí- 
culas y arbúsculos (con base en el método 
de tinción propuesto por Phillips y Hay- 
man, 1 970), mientras que la evaluación de 
colonización se realizó por el método de 
intersección propuesto por McGonigle et 
al. (1990). Para evaluar el comporta- 
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miento de estas últimas variables, se rea- 
lizó un análisis de varianza (p <_0,05) para 
analizar el comportamiento de cada una 
de las fuentes de inoculo. 

Los datos del número de esporas 
fueron transformados logarítmicamente 
(log x+1) debido a la gran variación, mien- 
tras que los datos porcentuales de coloni- 
zación tuvieron una transformación 
angular o de Bliss por presentar distribu- 
ción binomial; esto con el fin de cumplir 
con los supuestos del análisis de varianza, 
en el cual se requiere que respondan a 
una distribución normal como lo sugieren 
Steel y Torrie (1980). 


RESULTADOS 

Se encontraron diferencias en las variables 
de crecimiento, porcentaje de colonización 
y número de esporas tanto en la sap como 
en la sms (Figs. 1, 2, 3 y 4). El análisis de 
varianza arrojó únicamente diferencias 
estadísticas (p<0,05) a la aplicación de ino- 
culantes nativos procedentes de la sms 
para todas las variables, excepto para el 
número de esporas (Tabla 1). 

Los tratamientos que involucraron 
inoculantes nativos procedentes de la SMS, 
fueron los que mostraron los máximos 
valores en el desarrollo inicial de C. odo- 
rata (Figs. 1A, 1 C, 1 E, 2Ay 2C). La mayor 
respuesta la aportó el tratamiento que invo- 
lucró suelo nativo sin esterilizar (T4), 
seguido del tratamiento inoculado con raí- 
ces (T2). El tratamiento que se inoculó con 
100 esporas aleatorias de este ecosistema 
(TI), promovió un escaso desarrollo pero 
siempre por arriba del tratamiento testigo 
(T3). La diferencia estadística entre los tra- 
tamientos la aportó el tratamiento testigo 
con respecto a aquellos que involucraron 
algún tipo de inoculante. 

Aunque los tratamientos con inocu- 
lantes de la sap no mostraron diferencias 


estadísticas en el análisis de varianza, se 
observó un comportamiento muy similar al 
reportado en la sms en las variables de 
crecimiento (Figs. IB, ID, 1 F, 2B y 2D). 

La relación peso seco aéreo/peso 
seco radicular mostró una mayor biomasa 
en hojas y tallos con respecto a la bio- 
masa radicular. El tratamiento que involu- 
cró suelo nativo sin esterilizar (T4), fue el 
que presentó una relación menor compa- 
rada con los otros tratamientos (Figs. 2E y 
2F); esto es atribuido a la cantidad de bio- 
masa radicular que estimuló los mayores 
desarrollos en las etapas tempranas de C. 
odorata. 

La colonización de hma al interior de 
la raíz, marcó una alta diferencia estadís- 
tica entre los tratamientos de ambos sue- 
los, especialmente por la cero inoculación 
en T3 (Tabla 1); sin embargo, cuando se 
compararon por separado a TI, T2 y T4, 
no se observaron diferencias estadísticas 
tanto en hifas, vesículas y arbúsculos. Los 
tratamientos que involucraron alguna 
fuente de inoculo fueron los que presenta- 
ron un mayor porcentaje de colonización 
por hifas con respecto a la colonización por 
vesículas y arbúsculos encontradas al 
interior de la raíz del cedro rojo. El porcen- 
taje de colonización de vesículas fue el 
más bajo en los tratamientos empleados, 
mientras que los arbúsculos presentaron 
altos porcentajes de colonización, pero 
siempre por debajo de la colonización por 
hifas (Fig. 3). 

La propagación del número de espo- 
ras, marcó una alta diferencia estadística 
entre los tratamientos de ambos suelos 
(Tabla 1), especialmente por la cero inocu- 
lación en T3; sin embargo, al igual que en 
los porcentajes de colonización, al 
comparar por separado TI , T2 y T4, no se 
encontraron diferencias estadísticas en 
ambos tratamientos. En la figura 4 se 
observa que el mayor número de esporas 
lo presentó el tratamiento de suelo nativo 
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Tabla 1 . Análisis de varianza que muestra la respuesta en crecimiento de Cedrela 
odorata, la población de esporas y los porcentajes de colonización al utilizar diferentes 
fuentes de inoculo procedentes de dos ecosistemas tropicales (p<0,05). 



Selva mediana subperennifolia 

Selva alta perennlfolia 


F 

P 

F 

P 

Altura (m) 

11,59 

0,001*** 

1,76 

0,212NS 

Diámetro (cm) 

8,98 

0,002*** 

1,04 

0,41 3NS 

Número de folíolos 

12,08 

0,001*** 

0,95 

0,451 NS 

Biomasa aérea (g) 

8,26 

0,003*** 

1,69 

0,227NS 

Biomasa raíz (g) 

4,74 

0,021** 

1,36 

0.306NS 

Número de esporas/100 g de 
suelo 

30,38 

0,000*** 

80,82 

0,000*** 

% de hifas 

16,38 

0,000*** 

11,25 

0,001*** 

% de vesículas 

3,89 

0,037* 

3,32 

0,06NS 

% de arbúsculos 

14,89 

0,000*** 

11,59 

0,001*** 


*Significativo **Muy significativo ***Altamente significativo NS No significativo 


sin esterilizar (T4), tanto en la sap como 
en la sms. Los tratamientos 1 y 2 no supe- 
raron las 200 esporas en 100 g de suelo; 
sin embargo, fueron suficientes para 
lograr un porcentaje de colonización entre 
50% y 80% para estos factores de evalua- 
ción (Fig. 4). 

DISCUSIÓN 

Los hma forman diferentes fuentes de 
inoculo que son indispensables para for- 
mar la nueva asociación micorrízica; estas 
fuentes pueden ser: esporas, hifas extra- 
radicales y raíces colonizadas por hifas, 
vesículas y arbúsculos; sin embargo, no 
todos los hongos tienen la misma capaci- 


dad de formar estas asociaciones con las 
estructuras antes mencionadas (Klirono- 
mos y Hart, 2002). Esta característica se 
corroboró con los resultados de esta 
investigación, demostrando que ambas 
fuentes pueden tener diferente potencial 
para la producción de plántulas en vivero. 

Generalmente, el uso más común 
que se le ha dado a estas fuentes de 
inoculo, ha sido para la propagación de 
especies micorrízicas en cultivos trampa 
(Bellgard, 1992). Básicamente, estas 
fuentes pueden ser útiles para replicar 
cepas nativas con alto potencial y que 
posteriormente pueden ser utilizadas en la 
producción de plántulas con fines de res- 
tauración o de plantación comercial. 
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Figura 1 . Crecimiento en altura, diámetro y número de folíolos de C. odorata, utilizando 
diferentes fuentes de inoculo procedentes de dos ecosistemas tropicales. Los datos 
graficados representan la media ± error estándar en los días posteriores a la siembra 

(n=48). 


Relación peso seco aereo/peso seco 

radicular Biomasa seca de raíz (g) Biomasa seca aerea (g) 
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Figura 2. Crecimiento en biomasa seca de raíz y biomasa seca aérea de C. odorata, 
utilizando diferentes fuentes de inoculo procedentes de dos ecosistemas tropicales. Los 
datos graficados representan la media ± error estándar (n=48). 
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Selva mediana subperennifolia Selva alta perennifolla 

■ HIFAS «VESÍCULAS «ARBÚSCULOS 


Figura 3. Porcentaje de colonización por hitas, vesículas y arbúsculos en Cedrela 
odorata, utilizando diferentes fuentes de inoculo procedentes de dos ecosistemas 
tropicales. Los datos graficados representan la media ± error estándar (n=48). 



TI T2 T3 T4 

Selva mediana subperennifolia 


TI T2 T3 T4 

Selva alta perennifolia 


Figura 4. Número de esporas presentes en Cedrela odorata, utilizando diferentes fuentes 
de inoculo procedentes de dos ecosistemas tropicales. Los datos graficados representan 

la media ± error estándar (n=48). 


Los resultados de este bioensayo 
demostraron que las mayores tasas de 
crecimiento inicial en C. odorata se 
observaron cuando interactuaron simbió- 
ticamente con cualquier fuente de inoculo 
de estos hma. Los tratamientos que invo- 
lucran suelo nativo estéril o ausencia de 
fuentes de inoculo, fueron los que desa- 
rrollaron una menor tasa de crecimiento; 
por ende, se concluye que la interacción 
planta-micorriza es esencial para lograr 
un mayor crecimiento en esta especie 


maderable, tal y como lo documentó 
Amador (2010) al inocular este tipo de 
plantas con propágulos de Rhizophagus 
intraradices (antes Glomus intraradices). 

Aún con la interacción planta-mico- 
rriza, se observó que las plántulas son 
capaces de responder de manera dife- 
rente cuando se utiliza uno u otro tipo de 
inoculo. Esta diferencia es atribuida a la 
capacidad que tienen ciertos inoculantes 
del suelo para lograr una pronta coloniza- 
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ción en un tiempo más corto. Lo anterior 
fue documentado por Bellgard (1992), 
quien encontró que los fragmentos de raí- 
ces estimulan una rápida inoculación con 
respecto a las esporas o ciertos fragmen- 
tos hifales, lo cual puede traducirse en un 
mayor crecimiento inicial de las plántulas 
establecidas en vivero. 

En general, se observó un mayor 
efecto en el crecimiento de C. odorata 
cuando se utilizó suelo nativo no esterili- 
zado; esto se le atribuye a la gran varie- 
dad de fuentes de inoculo (esporas, 
vesículas, arbúsculos, hifas extra e intra- 
radicales) que desarrollan naturalmente 
los ecosistemas tropicales, lo que le con- 
fiere una rápida colonización de acuerdo 
con lo comentado anteriormente. 

Los porcentajes de colonización al 
final de la evaluación no tuvieron relación 
entre las diferentes fuentes de inoculo. 
Esta situación es atribuida a la pronta 
colonización que se da en las etapas tem- 
pranas de la planta, la cual aumenta con- 
forme van desarrollando las nuevas 
raíces. Es por esta razón que los porcen- 
tajes observados en estos tratamientos 
mostraron un comportamiento muy similar 
entre uno y otro al final del bioensayo. 

La producción de esporas por gramo 
de suelo no tiene relación con los porcen- 
tajes de colonización, lo cual puede estar 
asociado en mayor medida a la especie o 
especies micorrízicas que interactuaron 
durante el bioensayo; sin embargo, la baja 
producción de estas esporas deberá ser 
analizada una vez que se determine la 
efectividad o no efectividad relacionada 
con cada una de las especies en estudios 
de asociación. 

Por último, los hma desempeñan un 
papel importante en la sucesión ecológica 
(White et al., 2008); por ello, el uso de ino- 
culantes nativos es ideal para fomentar la 
interacción entre estos hongos y C. odo- 


rata. Además, estas fuentes pueden ser 
utilizadas en el futuro para la producción 
de plántulas con fines de reforestación o 
plantaciones forestales comerciales en 
áreas deforestadas del trópico mexicano, 
lo que puede propiciar una pronta recupe- 
ración de estas áreas y, por ende, una 
mayor adaptación, tal y como lo sugiere 
Alien et al. (2003) en sus resultados con 
algunas especies tropicales. 


CONCLUSIONES 

El uso de cepas nativas de hma proceden- 
tes de ambientes naturales, estimulan un 
mayor desarrollo inicial en C. odorata ; 
además, su implantación en zonas alta- 
mente degradadas podría establecer una 
mayor adaptación como lo han sugerido 
otros autores; por ello, se recomienda pro- 
pagar a estas cepas para generar paque- 
tes tecnológicos viables en la producción 
de plántulas en vivero o invernadero. 

Las raíces colonizadas son una de 
las fuentes de inoculo poco utilizadas en 
estudios de asociación; sin embargo, pue- 
den ser útiles para propagar especies 
micorrízicas asociadas directamente de 
sus áreas naturales de distribución. 

El uso de suelo nativo procedente 
directamente de las áreas naturales de 
distribución, en combinación con otros 
sustratos utilizados en vivero o inverna- 
dero, podría servir como fuente de inoculo 
para la asociación simbiótica entre los 
hongos nativos y las plantas de interés 
comercial. 

Las fuentes de inoculo provenientes 
de la sms estimularon un mayor desarrollo 
inicial en plántulas de C. odorata, lo cual 
explica la función que cumplen los hma en 
el mantenimiento de ecosistemas que 
aparentan tener un mayor estrés ante la 
sequía. 
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Es necesario replicar las especies de 
hma que se asociaron durante el bioen- 
sayo en cada uno de los ecosistemas tro- 
picales. Posteriormente, se recomienda 
realizar pruebas de efectividad y no efecti- 
vidad individual y en consorcio para la pro- 
ducción de plántulas de C. odorata. Esta 
situación permitirá obtener cepas con alto 
potencial para la inoculación de esta 
especie maderable y que posteriormente 
serán útiles para el desarrollo de paquetes 
tecnológicos. 
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Influencia de variables físicas en la 
productividad de Pinus arizonica y Pinus 
engelmannii en el sur de Chihuahua, México 

Influence of physical variables in the yield of Pinus arizonica 
and Pinus engelmannii in the south of Chihuahua, México 

Martín Martínez-Salvador, 1 * Ricardo David Valdez-Cepeda 2 

y Marín Pompa García 3 

RESUMEN 

Se utilizó el índice de sitio (IS) para identificar la influencia de variables de suelo, clima y relieve 
sobre el crecimiento de P. arizonica y P. engelmannii en el sur de Chihuahua. Los modelos de IS fueron 
desarrollados a partir de información de análisis troncales usando la función de Schumacher y el 
método de curva guía. Adicionalmente, se desarrolló un muestreo sistemático usando parcelas de 
400 m 2 , donde se midieron altura y edad de árboles dominantes, así como variables de suelo, clima y 
relieve. Con los modelos de IS se identificó la calidad de estación de cada parcela muestreada. Se 
realizó un análisis de discriminante canónico (ADC), para identificar diferencias multivariadas entre 
calidades de estación, y un análisis de varianza, para identificar la influencia de las variables físicas en 
el crecimiento de las especies. El método de curva guía predijo con precisión el crecimiento de ambas 
especies (Pseudo R 2 > 0,90). El ADC mostró que existen diferencias multivariadas entre las calidades 
de estación, con alta influencia de la pendiente, disección vertical, profundidad de suelo y temperatu- 
ras. De acuerdo con el análisis de varianza (a < 0,05), los valores para las variables de mayor influencia 
sobre la mejor calidad de estación fueron la profundidad de suelo > 43 cm, pendiente menor a 12° ± 8 o , 
temperatura mínima promedio menor a 8 °C ± 3 °C y temperatura máxima promedio menor a 18 °C ± 

3 °C. Estos valores pueden ser de gran utilidad para construir modelos espaciales de productividad en 
los estados de Chihuahua y Durango, México. 

PALABRAS CLAVE 

Calidad de estación, crecimiento, índice de sitio, variables físicas. 

ABSTRACT 

In the south of Chihuahua México, site Índex (IS) models were used in order to identify the influ- 
ence of soil, climate and topography variables on P. arizonica and P engelmannii growth. The IS models 
were fitted by using the Schumacher equation and the guide curve with information from stems analysis. 
Additionally a systematic sampling method was carried out on 400 m 2 field plots. Dominant tree height 
and age, and physical variables of the sites were measured. Site índex models were adjusted with tree 
height and age data for both species. With site index models the productivity of each sampled plot was 
identified. An analysis of canonical discrimínate was carried on in order to identify multivariate differ- 
ences in the productivity levels, and an analysis of variance was conducted to detect differences per 
each variable studied. Site index models accuracy were statistically significant with a pseudo r 2 > 0,9. 

The ADC showed differences in the levels of productivity, where slope, vertical dissection, soil deep and 
temperatures were the most important variables. According to the analysis of variance (a < 0,05), the 
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valúes for places whit highest productiv- 
ity, and those variables which have more 
influence over the growth of both species 
are: soil depth > 43 cm, slope lower than 
12° ± 8 o , average of minimum tempera- 
ture of 8 °C ± 3 °C, máximum tempe- 
raturas < 18 °C ± 3 °C. These valúes 
could be highly useful to build spatial 
models of productivity in the northwest of 
Mexican forest. 

KEYWORDS 

Station quality, growth, site Índex, physi- 
cal variables. 


INTRODUCCIÓN 

En los bosques de la Sierra Tarahumara 
del sur de Chihuahua, México, crecen 
poblaciones nativas de Pinus arizonica y 
Pinus engelmannii, las cuales son aprove- 
chadas para la extracción de madera. El 
manejo sustentable de los bosques 
requiere del conocimiento de los factores 
ambientales que determinan el creci- 
miento de las especies forestales (Palahí 
et al., 2004). La influencia de estos facto- 
res sobre la producción y productividad 
del bosque generalmente se evalúa a tra- 
vés de la construcción de modelos de 
índice de sitio, como una medida indirecta 
de la calidad de estación (Madrigal, 2002); 
sin embargo, en pocas ocasiones se eva- 
lúa la influencia directa de las variables 
físicas sobre el crecimiento del arbolado 
(Heiberg y White, 1956; Martínez et al., 
2006). 

La calidad de estación se refiere a la 
capacidad inherente de un territorio para 
hacer crecer los árboles que sostiene y 
usualmente se utilizan modelos que rela- 
cionan edad-altura de árboles dominantes 
para calificarla. Estos modelos son cono- 
cidos como índice de sitio (Clutter et al., 
1983), los cuales están determinados por 
la suma de los factores como suelo, topo- 
grafía, clima y factores bióticos, que inte- 
rrelacionados afectan la productividad de 
un bosque (Spurr, 1952; Spurr y Barnes, 
1982; Clutter et al., 1983; Daniel et al., 


1982; Musalem et al., 1986; Pritchett, 
1986; Deo, 2013). 

En trabajos realizados en México se 
describen las relaciones entre la altitud y 
la distribución de especies de pino (Egui- 
luz, 1977; Eguiluz, 1982; Madrigal, 2002) 
y recientemente se han desarrollado 
modelos bioclimáticos para reconocer 
patrones de distribución potencial de 
especies de pináceas y para la cuantifica- 
ción de características climáticas donde 
las especies prosperan o podrían desarro- 
llarse con mayor productividad (Téllez et 
al., 2004; Peters et al., 2013). 

Existe un creciente interés por cono- 
cer la influencia de variables físicas sobre 
el crecimiento de especies y productividad 
de los bosques (Martínez, 2008; Medlyn et 
al., 2011), ya que el conocimiento de esta 
información es importante para generar 
cartografía de potencial productivo de 
diferentes áreas mediante el uso de herra- 
mientas de Sistemas de Información Geo- 
gráfica (sig) (Martínez etal., 2006; Latta et 
al., 2009). 

Estudios ecológicos muestran que se 
han encontrado relaciones entre el creci- 
miento de coniferas y la profundidad de 
suelo, altitud, pendiente y exposición 
(Castaños, 1962; Griffiths et al., 2009). 
Otros trabajos mencionan que las propie- 
dades físicas que más influyen en el creci- 
miento de algunas coniferas son pH, 
material muerto y contenido de limo en la 
capa de 50 cm de suelo (Rodríguez, 1982; 
Stoler y Relyea, 2011). Estos anteceden- 
tes pueden ser usados para construir 
modelos espaciales a partir de valores de 
óptimas condiciones para el desarrollo de 
las especies (Sánchez y Chacón, 1986; 
Martínez, 2008). 

Con la finalidad de caracterizar la 
influencia de variables físicas de los facto- 
res de clima, suelo y relieve en los bos- 
ques del sur de Chihuahua, México, se 
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clasificaron sitios con diferente productivi- 
dad (alta, media y baja), usando modelos 
de índice de sitio. Los resultados son de 
gran utilidad para generar y validar carto- 
grafía de los bosques en el norte de 
México y para implantar mejores prácticas 
de manejo en los predios en los que se 
realiza cosecha de madera. 


MATERIALES Y MÉTODO 

El estudio se realizó en los bosques del 
municipio de Guadalupe y Calvo, en la 
Sierra Tarahumara del sur de Chihuahua, 
México en los años 2006 y 2007 (Fig. 1). 

El municipio de Guadalupe y Calvo 
se encuentra entre los 26°41’ y 25°34’ de 
latitud norte y 106°25’ y 107°42’ de longi- 
tud oeste, tiene una superficie aproximada 
de 904 000 hectáreas. El clima que predo- 


mina en el municipio es templado subhú- 
medo con lluvias en verano (García, 
1981). Los principales suelos del munici- 
pio son litosol, regosol, xerosol, feozem, 
cambisol y planosol (cetenal, 1972); los 
cuales, se caracterizan por ser muy delga- 
dos y con alto porcentaje de pedregosi- 
dad. Los tipos de vegetación de la región 
son bosques templados, selva baja cadu- 
cifolia, pastizales y matorrales (Rze- 
dowsky, 1978). La principal actividad pro- 
ductiva en el municipio es la silvicultura 
con una cosecha anual de aproximada- 
mente 300 000 m 3 de madera en rollo, lo 
cual representa 25% de la madera que se 
produce en el estado y 3,7% de la produc- 
ción nacional. 

Para el ajuste de modelos de índice 
de sitios se realizaron análisis troncales a 
un grupo de 40 árboles dominantes de P 
arizonica y 42 de P engelmarmii siguiendo 




Figura 1 . Área de estudio. 
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las metodologías propuestas por Mas 
(1970) y el fpc (1998); las cuales con- 
sisten en examinar secciones transver- 
sales del fuste, seleccionadas a 30 cm, 
1 ,30 m y cada 2 m hasta alcanzar 10 cm 
de diámetro, con medición de la punta. 
Las secciones transversales fueron anali- 
zadas mediante el conteo y medición de 
anillos para integrar la base de datos con 
información de edad, diámetro y altura 
para el posterior análisis mediante pro- 
cedimientos de regresión no lineal. 

Para calificar la calidad de estación 
en los bosques del municipio se mues- 
trearon 220 sitios circulares de 400 m 2 dis- 
tribuidos en un diseño sistemático de 5 km 
entre sitios. En cada sitio se anotaron los 
datos de ubicación geográfica y la altura 
sobre el nivel medio del mar, y se tomó 
información de altura y edad de tres árbo- 
les dominantes de las especies P arizo- 
nica y P. engelmannii. Además se midieron 
variables físicas de los factores suelo, 
relieve y clima. 

Para el caso del suelo se midieron la 
profundidad de suelo (cm), pedregosidad 
(%), suelo desnudo (%), grosor de la capa 
de mantillo (grosor de mantillo) (cm), gro- 
sor de la capa de hojarasca (grosor de 
hojarasca) en centímetros. Para el relieve 
se midieron la pendiente (°), exposición 
(Azimut), disección horizontal (densidad 
de arroyos por kilómetro cuadrado) y 
disección vertical (diferencia en altura 
sobre el nivel del mar en m/km 2 ), las dos 
últimas variables se derivaron del análisis 
del Modelo Digital de Elevación esc 
1:50 000. 

Finalmente, se obtuvieron datos de 
variables de clima como la temperatura 
mínima, temperatura máxima y precipita- 
ción. Las variables de clima se generaron 
mediante modelos de interpolación en sig 
con información de estaciones meteoroló- 
gicas ubicadas en la región de influencia 
del área de estudio que comprende los 


estados de Durango, Chihuahua, Sonora 
y Sinaloa (imta, 2006). 


Análisis de la información 

Con la información de altura y edad prove- 
niente de los análisis troncales se ajusta- 
ron funciones de índice de sitio para las 
especies Pinus arizonica y Pinus engel- 
mannii utilizando el modelo de Schuma- 
cher con el método de curva guía. 

Modelo Schumacher 

H = (3 0 e^ E 

donde: 

H = Altura dominante (m) 

E = Edad (años) 

| 3 0 = Parámetros de los modelos 

Cada sitio de muestreo fue calificado de 
acuerdo con su calidad de estación utili- 
zando la información de edad y altura de 
los árboles dominantes de las parcelas de 
muestreo y los modelos de índice de sitio 
ajustados para ambas especies. 

Se establecieron tres niveles de pro- 
ductividad (calidad de estación) I. Alta, II. 
Media y III. Baja, y posteriormente se rea- 
lizó un análisis de discriminante canónico 
(adc) para detectar diferencias multivaria- 
das entre las tres calidades de estación 
a partir de la base de datos de variables 
físicas. 

Se consideró el método de la distan- 
cia cuadrada Mahalanobis y la prueba de 
Lamda Wilks como estadísticos de prueba 
para la toma de decisiones sobre las fun- 
ciones de discriminante canónico (fdc) 
(Hair et al., 1998). 

Una vez identificadas diferencias 
multivariadas entre calidades de estación, 
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y posterior a una estandarización de la 
base de datos, se realizó un análisis de 
varianza para identificar diferencias a 
nivel de variables individuales, con espe- 
cial atención en las variables que definie- 
ron las funciones de discriminante 
canónico (fdc). Finalmente, se realizó una 
comparación de medias usando el proce- 
dimiento de Tukey con un valor de signifi- 
cancia de 0,05. 

Para el análisis de varianza se usó 
un modelo completamente aleatorio. 

Y,-M+T¡ + x , 

donde: 

Y j = variable resultado del modelo 
M j = media generalizada 
T j = efecto de tratamiento (calidad de 
estación) 

X ikj = error experimental. 

RESULTADOS Y DISCUSIÓN 

Los modelos de índice de sitio para P. arizo- 
nica, y P engelmannii fueron estadística- 
mente aceptables (Tabla 1) y altamente 
predictivos con una pseudo R 2 de 0,98 para 
P. añzonica y 0,92 para P. engelmannii. 
Estos modelos permitieron clasificar la base 


de datos de sitios de muestreo en función 
de tres calidades de estación (niveles de 
productividad) (Figs. 2 y 3). Los valores 
observados que se muestran en las figuras 
2 y 3 corresponden a las alturas de los árbo- 
les dominantes medidos en campo en los 
sitios de muestreo, mientras que las curvas 
corresponden a los modelos de IS ajusta- 
dos con el modelo de Shumacher, usando el 
método de curva guía a partir de pares de 
datos edad-altura de los análisis troncales. 

Estos modelos son útiles para prede- 
cir el crecimiento de los árboles en el 
futuro y de esta manera realizar adecua- 
damente la planificación del manejo, pro- 
yecciones del inventario y determinación 
del potencial de cosecha (Clutter et al., 
1983; Palahí et al., 2004; Martínez et al., 
2007); sin embargo, con estos modelos es 
importante construir cartografía, la cual 
permitirá integrar con mayor facilidad el 
índice de sitio a la planeación y manejo de 
bosques, por lo que es necesario conocer 
la influencia de las variables ecológicas 
como suelo, pendiente y relieve sobre el 
crecimiento de las especies (Spurr y Bar- 
nes,1982; Martínez, 2008). 

Una de las pruebas que resultan de 
gran utilidad para identificar diferencias en 
la influencia de variables o grupos de 
variables sobre las diferentes condiciones 


Tabla 1 . Modelos de Schumacher de índice de sitio ajustados por el método de curva guía. 


Model 

Estimadores 

Int. Con f id ( 3 1 

Pseudo 

R 2 

Pr> F 

Po 

Pi 

Inf. 

Sup. 

Pin us arizonica 

26,98 

30,96 

20,40 

41,53 

0,98 

0,0001 

Pinus engelmannii 

26,52 

34,44 

22,03 

46,84 

0,92 

0,0001 
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Figura 2. Modelo de Schumacher con el método de curva guía en índice de Sitio de 
Pinus arizonica Engelm., en Guadalupe y Calvo, Chihuahua, México. 

Se analizó el crecimiento de los árboles y se estableció la clasificación de las curvas de 
índice de sitio (IS) a una edad base de 80 años. En función de la edad base el IS < 15,5 m 
representa la curva de los sitios de productividad baja que agrupa a los árboles que 
alcanzan una altura de 15,5 m, el IS > 15,5 < 19,5, representa la curva de los sitios de 
productividad media que agrupa a los árboles que alcanzan alturas entre 15,5 m y 19,5 m 
y el IS > 19,5 m, representa la curva de los sitios de productividad alta que agrupa a los 
árboles que alcanzan alturas superiores a los 19,5 m. 


de un sistema es el análisis de determi- 
nante canónico (adc). En este estudio al 
realizar el adc, la distancia cuadrada de 
Mahalanobis indicó que las calidades de 
estación son diferentes entre sí (P < 0,05). 
Por su parte, la prueba de Lambda Wilks 
también mostró diferencias multivariadas 
entre las tres calidades de estación 
(P < 0,05). 

La correlación canónica y el vector 
raíz indicaron que la primer correlación 
canónica fue mayor que la segunda 
(P < 0,05) con un coeficiente de correla- 
ción R = 0,50, y un valor del vector raíz de 
0,60, mientras que el porcentaje de la 
variación explicada con la primer función 
discriminante fue de 0,83. 


La estructura canónica total (Tabla 2) 
muestra que las variables pedregosidad, 
pendiente, disección vertical y exposición 
están negativamente relacionadas con la 
primera variable canónica, mientras que 
las variables altura sobre el nivel del mar 
(asnm), profundidad de suelo, tempera- 
tura máxima y mínima, mantillo y precipi- 
tación, están relacionadas positivamente. 
La primer función de discriminante canó- 
nico muestra que la separación entre nive- 
les de productividad está dada por la 
disección vertical (R = 0,60), precipitación 
(R = 0,48), profundidad de suelo 

(R = 0,43) y la pendiente (R = 0,28). 

El análisis de discriminante canónico 
indica que las calidades de estación son 
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Figura 3. Modelo de Schumacher con el método de curva guía en índice de Sitio de 
Pinos engelmannii en Guadalupe y Calvo, Chihuahua, México. 

Se analizó el crecimiento de los árboles, y se estableció la clasificación de las curvas de 
índice de sitio (IS) a una edad base de 80 años. En función de la edad base el IS<14 m, 
representa la curva de los sitios de productividad baja que agrupa a los árboles que 
alcanzan una altura de 14 m, el IS>14 <20 representa la curva de los sitios de productividad 
media que agrupa a los árboles que alcanzan alturas entre 14 m y 20 m y el IS>20 m, 
representa la curva de los sitios de productividad alta que agrupa a los árboles que 
alcanzan alturas superiores a los 20 m. 


multivariadamente distintas (P < 0,005), 
considerando la primer función canónica 
(Fig. 4). Por su parte, la figura 5 muestra 
gráficamente estas diferencias entre las 
tres calidades de estación. 

Con base en el DCA se observa 
mayor influencia de las variables topográ- 
ficas y climáticas sobre la productividad. 
Al respecto, estudios recientes muestran 
que las características climáticas influyen 
considerablemente sobre la productividad 
del bosque, donde las temperaturas jue- 
gan un papel importante en la generación 
de microambientes que aceleran la 


descomposición del material orgánico y la 
consecuente liberación de nutrientes; 
mientras que las relaciones entre la pre- 
cipitación y las densidades de arbolado 
incrementan la disponibilidad de humedad 
en bosques secos (Deo, 2013; Raz-Yas- 
eef etal., 2010; Karhu etal., 2010; Latta et 
al., 2010), como los de la región de 
estudio. 

Una vez identificadas diferencias 
entre las calidades de estación, es impor- 
tante conocer la diferencia de cada varia- 
ble entre las calidades de estación con la 
finalidad de identificar el grupo de varia- 
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bles que podrían apoyar al entendimiento 
espacial de la productividad. 

Los resultados del análisis de 
varianza realizado al grupo de variables 
físicas estima que la variable asnm no 
tuvo diferencias significativas (a < 0,05) 
entre los sitios con diferente productividad 
(Tabla 3). Al respecto, Spurr y Barnes 
(1982) y Eguiluz (1982) mencionan que 
los factores ambientales que más influyen 
sobre la productividad del bosque son los 
climáticos, edáficos y biológicos, mientras 


que la altura, aunque influye en la delimi- 
tación del área óptima de desarrollo de las 
especies, tiene una mayor influencia 
sobre su distribución. 

La pendiente, disección vertical, 
exposición, temperatura mínima, tempe- 
ratura máxima y profundidad de suelo pre- 
sentaron diferencias significativas (a < 
0,06) (Tabla 3). Estos resultados coinci- 
den con estudios realizados en bosques 
del centro de México y noroeste de Esta- 
dos Unidos, donde se encontró que la 


Tabla 2. Coeficientes Canónicos Estandarizados (cce) y variables que influyen sobre la 
calidad de estación en el sur de Chihuahua, México. 


Variable 

CAN1 

CAN2 

Asnm 

0,18443 

-0,34208 

Pendiente 

-0,26774 

0,23286 

Profundidad de suelo 

0,42830 

0,15509 

Pedregosidad 

-0,27038 

0,08543 

Suelo desnudo 

0,09848 

-0,10904 

Grosor de mantillo 

0,26814 

0,01415 

Hojarasca 

0,10773 

-0,02531 

Disección horizontal 

0,11665 

0,26618 

Disección vertical 

-0,59968 

0,27544 

Exposición 

-0,42088 

-0,10468 

Temperatura mínima 

0,22402 

0,59289 

Temperatura Máxima 

0,49133 

0,23916 

Precipitación 

0,48251 

0,10239 


Variable canónica 1 (CAN1) dada por pendiente, profundidad de suelo, disección vertical, temperatura máxima y 
precipitación y la variable canónica 2 (CAN 2) por temperatura mínima, disección horizontal. 


Madera y Bosques 19(3), 2013:35-49 


43 


4 -- 


2 + 
Ian2 
O 

-2 - 


-11 


12 


1 


1 


1 


3 

3 


21 12 3 1 12 3 13 2 3 3 

1 12 2 31 2 1 322 332 333 

-2 2-2—1 1 331 1-2^—22232-32-3 3 

2 111 2 1 1 323 | 21 12 3 232 

22 I 2222 2 32331 22 3 2 


2 3 

2 


-4 ■■ 


-3 


-2 


-1 


Canl 


Figura 4. Separación de las variables canónicas. 


Separación de grupos a partir de valores canónicos de calidad de estación. Los niveles 
corresponden a 1= Calidad de estación mala, 2 = Calidad de estación media y 3 = Calidad 
de estación alta. 


productividad se ve afectada por la pen- 
diente, exposición y profundidad de suelo 
(Griffiths et al., 2009; Gómez Tagle, 1985; 
Chávez y Gómez Tagle, 1985). 

Las variables pedregosidad, porcen- 
taje de suelo desnudo, grosor de mantillo, 
grosor de hojarasca y precipitación no 
mostraron diferencias significativas (a < 
0,05) entre las calidades de estación 
(Tabla 3). Esto puede deberse a la com- 
plejidad del territorio estudiado o a las 
actividades de intenso manejo que se 
desarrollan en el área de estudio, las cua- 
les podrían modificar las características 
del suelo, del material orgánico sobre el 
suelo y del material en descomposición 
(Nave et al., 2010; Poeplau et al., 2011; 
Vincent et al., 2010). Aunque diversos 
estudios han mostrado que los compo- 
nentes del suelo ejercen una marcada 
influencia en la productividad del sitio y el 
desarrollo del arbolado; por ejemplo, 
Mader (1976) menciona que las caracte- 
rísticas del suelo explican alrededor de 
60% de la productividad, mientras que 


Sánchez y Chacón (1986) encontraron 
diferencias notables en rocosidad y pro- 
fundidad de suelo en relación con la pro- 
ductividad de Pin us arizonica en el norte 
de Chihuahua, México, por su parte Nar- 
váez (1990) encontró que el pH, el grosor 
de mantillo y las características del hori- 
zonte A influyen sobre el potencial produc- 
tivo del bosque. Sin embargo, a diferencia 
de las variables de clima y relieve, las 
variables de suelo pueden ser afectadas 
fuertemente por disturbios naturales o por 
factores antropogénicos como las cose- 
chas y el manejo (Heiberg y Whithe, 
1956), lo cual dificulta la identificación de 
influencia de macrovariables, sobre todo 
en estudios de tipo regional. Ante esto, se 
han aplicado estrategias como el uso de 
variables que caracterizan al relieve para 
definir características del suelo. Ésta ha 
sido una práctica útil para caracterizar 
suelos en extensas regiones (Thiers y 
Gerding, 2007). 

La tabla 4 muestra los valores medios 
para cada variable en cada calidad de 
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Tabla 3. Análisis de varianza comparación de variables ecológicas entre sitios 

(productividad alta, media y baja) 


Variable 

Suma de 
cuadrados 

Cuadrado medio 

Fe 

Significancia 

Altura sobre el nivel del mar (m) 

461386,09467 

230693,04733 

2,23822 

0,109659 

Pendiente (grados)** 

0,21467 

0,10734 

3,99936 

0,020047 

Profundidad de suelo (cm)* 

527,35401 

263,67700 

0,69847 

0,012343 

Pedregosidad (%) 

0,15431 

0,07716 

2,33223 

0,100090 

Suelo desnudo (%) 

0,00240 

0,00120 

0,05643 

0,945152 

Grosor de mantillo (cm) 

4,31110 

2,15555 

2,03871 

0,133350 

Grosor de hojarasca (cm) 

0,07070 

0,03535 

0,64367 

0,526795 

Disección horizontal (m/km 2 ) 

665331,51108 

332665,75554 

0,59109 

0,554842 

Disección vertical (m/km 2 )** 

170835,01070 

85417,50535 

11,22476 

0,000026 

Exposición* 

5,71505 

2,85753 

2,82671 

0,061883 

Temperatura mínima (°C)** 

37,37838 

18,68919 

3,26125 

0,040670 

Temperatura máxima (°C)** 

86,15165 

43,07583 

5,20779 

0,006349 

Precipitación (mm) 

1219,02285 

609,51142 

0,04307 

0,0957852 


estación. Estos valores pueden ser utiliza- 
dos para crear modelos espaciales 
(mapas) para la elaboración de un modelo 
general de productividad en la región de 
estudio. La temperatura mínima y la tem- 
peratura máxima mostraron sensibles 
diferencias. Por su parte, no se identifica- 
ron diferencias significativas en el análisis 
para la variable precipitación. 

Al respecto, Anaya et al. (1980) 
señalan que las principales variables que 
influyen sobre el crecimiento en altura y 
diámetro de las coniferas son las climáti- 


cas, como las temperaturas extremas y la 
precipitación. No obstante, en territorios 
con homogeneidad en las variables climá- 
ticas, como en el caso de la región de 
estudio, es complicado encontrar diferen- 
cias significativas. 

Los valores de las variables climáti- 
cas (Tabla 4) también pueden ser usados 
para reclasificar modelos espaciales para 
definir regiones con diferente productivi- 
dad. Estos datos podrían sumarse a las 
variables de relieve y suelo para desarro- 
llar un modelo espacial de productividad. 
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Tabla 4. Medias y desviaciones estándar de las variables físicas por sitios 

(productividad alta, media y baja) 


Variable 

P. Baja 

P. Media 

PAIta 

Media 

Dev. Std. 

Media 

Dev. Std. 

Media 

Dev. Std. 

Pendiente (grados) 

37 

15 

30 

15 

12 

8 

Profundidad de suelo (cm) 

25,82 

17,49 

35,48 

18,73 

43,73 

11,05 

Pedregosidad (%) 

0,29 

0,2 

0,22 

0,17 

0,21 

0,18 

Suelo desnudo (%) 

0,12 

0,1 

0,13 

0,14 

0,13 

0,17 

Grosor de mantillo (cm) 

1,51 

0,84 

1,64 

0,84 

1,9 

1,27 

Grosor de hojarasca (cm) 

0,58 

0,22 

0,63 

0,22 

0,64 

0,25 

Disección horizontal (m/km 2 ) 

1833,65 

781,46 

1704,69 

730,78 

1828,79 

752,79 

Disección vertical (m/km 2 ) 

311,23 

90,86 

238,51 

84,51 

234,25 

87,89 




NW, 




Exposición (rumbo) 

S, SE 


SySW 


N, NWyZ 


Temperatura mínima (°C) 

5 

2 

4 

1 

8 

3 

Temperatura máxima (°C) 

23 

3 

23 

2 

18 

3 

Precipitación (mm) 

794 

43 

897 

92 

901 

116 


CONCLUSIONES 

El uso de modelos de índice de sitio per- 
mitió definir niveles de productividad en la 
zona de distribución de Pinos arizonica y 
P. engelmannii, así como identificar las 
variables con mayor relación sobre la pro- 
ductividad de los sitios en la región sur del 
estado de Chihuahua. 

Los análisis de discriminante canó- 
nico y de varianza identificaron diferencias 
significativas en las variables físicas (pro- 
fundidad de suelo, pendiente, disección 
vertical, temperatura mínima y tempera- 


tura máxima). Los promedios e intervalos 
de estas variables pueden ser de gran 
utilidad para el desarrollo de cartografía 
de potencial productivo en la región de 
estudio. 

Las variables pedregosidad, porcen- 
taje de suelo desnudo, grosor de mantillo 
y grosor de hojarasca, no mostraron dife- 
rencias significativas entre los estratos 
con diferente productividad. Esto puede 
deberse a que son variables fáciles de 
modificar por acciones de manejo de los 
recursos. 
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Figura 5. Diferencias multivariadas por calidad de estación 
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Efecto de la micorrización 
en el establecimiento de 
Pinus pseudostrobus en cárcavas 

Effect of mycorrhyzation in the establishment 
of Pinus pseudostrobus in gullies 

Mariela Gómez-Romero ,* 1 Javier Villegas , 2 
Cuauhtémoc Sáenz-Romero 3 y Roberto Lindig-Cisneros 4 

RESUMEN 

Se evaluó el establecimiento de Pinus pseudostrobus Lindl., en tres condiciones de propagación 
y trasplante en cárcavas desprovistas de vegetación en la localidad de Atécuaro, municipio de Morelia, 
Michoacán. La primera condición corresponde a plantas de procedencia conocida, germinadas y pro- 
pagadas en condiciones controladas (estratificadas a 4 °C, germinadas a 25 °C con 12 h luz y propa- 
gadas con riego cada tercer día), para contar con dos lotes de plantas, uno con inoculo y otro sin 
inoculo del hongo ectomicorrícico Pisolithus tinctorius (Pers.) Coker et Couch. Como una segunda 
condición, ambos lotes fueron trasplantados al sitio de plantación a los 15 meses de edad con adición 
de sustrato fértil (fibra de coco, corteza de árbol y agrolita) en las cepas. En otro sitio muy cercano, una 
tercera condición fue una reforestación convencional con plantas de la misma especie pero sin contro- 
lar las variables antes descritas. Se evaluó altura, diámetro a la altura de la base, cobertura, número de 
ramas y porcentaje de supervivencia. Los resultados muestran que los pinos que fueron plantados bajo 
las primeras condiciones presentaron mayor crecimiento y supervivencia que los pinos de la segunda 
reforestación (tercera condición) (P < 000,1). La supervivencia de los pinos inoculados fue de 86%, la 
de la misma plantación sin inoculo 62% y la de la segunda reforestación 30%. La inoculación se plantea 
como una estrategia a aplicar en proyectos de restauración de sitios severamente degradados y ero- 
sionados, planteándose la necesidad de una inversión que incremente el éxito en los programas de 
reforestación. 

PALABRAS CLAVE 

Erosión, Pisolithus tinctorius, propagación, reforestación, restauración ecológica. 


ABSTRACT 

The establishment of Pinus pseudostrobus Lindl. was evaluated in three conditions of propaga- 
ron and transplantation in gullies devoid of vegetation in the locality of Atécuaro, Municipality of Morelia, 
Michoacán. The first condition correspond to plants of known provenance, germinated and propagated 
under controlled conditions (stratified at 4 °C, germinated at 25 °C with 12 h of light and watered every 
other day), with the aim of having plants without mycorrhizae and with the ectomycorrhizal fungus Piso- 
lithus tinctorius (Pers.) Coker et Couch. As a second condition, both groups where transplanted at 15 
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months of age with the addition of a fertile 
substrate (coconut fiber, bark and agro- 
lite). In a nearby site, a third condition 
was a reforestation with plants of the 
same species but without controlling the 
aforementioned variables. The results 
show that the pine trees that were plan- 
ted in the first condition performed better 
and had a higher survival than plants in 
the second reforestation (third condition) 
(P < 000,1). The survival of inoculated 
pines was 86%, in the absence of the fun- 
gus was 62% and 30% in the second 
reforestation. Inoculation is a good stra- 
tegy for restoration projects in severely 
degraded and eroded sites, although it 
increases the costs it also increases sur- 
vival and growth. 

KEYWORDS 

Erosión, Pisolithus tinctorius, propaga- 
ron, reforestation, ecological restoration. 


INTRODUCCIÓN 

Por el creciente deterioro de los recursos 
forestales en México, existe la necesidad 
de reforestar y restaurar eficientemente, 
tanto con fines comerciales como para 
recuperación de suelos (Cetina-Alcalá et 
al., 1999). El estado de Michoacán posee 
importantes recursos forestales por su 
cantidad, diversidad e importancia econó- 
mica. Ocupa el tercer lugar nacional en 
producción de madera (aproximadamente 
1 x 10 6 m 3 /año, después de Chihuahua y 
Durango) y el primer lugar nacional en 
producción de resina (35 000 t/año), no 
obstante que ocupa el sexto lugar nacio- 
nal en existencias maderables (cofom, 
2001). En cuanto a riqueza de especies, 
Michoacán tiene el quinto lugar nacional 
en biodiversidad (cofom, 2001) con 4831 
especies vegetales. Destaca, entre otros, 
el género Pinus con 16 especies (Conabio 
et al., 2005). Se ha publicado que Pinus 
pseudostrobus Lindl., es una de las espe- 
cies forestales de mayor importancia en 
México. Es un árbol de 15 m a 25 m, se 
encuentra en laderas con elevaciones de 
1600 msnm a 3200 msnm con tempera- 
tura de 9 °C a 40 °C con media de 14 °C, 
con precipitaciones de 500 mm a 2000 


mm con media de 1000 mm (Cruz de 
León, 1991). Se utiliza para la restaura- 
ción de suelos degradados (Conafor, 
2006), es buen productor de resina, 
ampliamente aprovechada en los estados 
del centro y sur del país. La madera es de 
buena calidad. Es una especie recomen- 
dable para plantaciones comerciales, tam- 
bién para su uso ornamental (Eguiluz, 
1978). Por otro lado, en el estado de 
Michoacán existen serios problemas de 
deforestación; se calcula que se pierden 
entre 30 000 ha y 40 000 ha de bosque 
por año. Las principales causas de la 
deforestación son: cambio de uso del 
suelo, incendios forestales y tala ilegal 
(cofom, 2001). En un estudio realizado 
con la finalidad de hacer una evaluación 
sobre los programas de reforestación en 
Michoacán y proponer estrategias que 
aumenten el éxito de establecimiento 
(Sáenz-Romero y Lindig-Cisneros, 2004), 
se menciona que existen serios proble- 
mas en cuanto al proceso de plantación y 
mantenimiento de las reforestaciones 
(Lindig-Cisneros et al., 2005), disminu- 
yendo significativamente la supervivencia 
de las plantas. Entre los principales pro- 
blemas detectados se encuentran la téc- 
nica de plantación (falta de capacitación 
del personal de campo, falta de herra- 
mientas adecuadas y falta de preparación 
del terreno previo a la plantación), fecha 
inadecuada de plantación (temporada de 
lluvias avanzada, limitando a las plantas 
para establecer un sistema radical sufi- 
ciente, antes del inicio de la temporada de 
sequía) y a la pobre calidad de la planta, 
que representa hasta 15% de la mortali- 
dad (Nienstaedt, 1994; Sáenz-Romero y 
Martínez-Palacios, 2000). 

Por otro lado, se ha documentado 
que en la restauración ecológica se debe 
dirigir al ecosistema a través de una serie 
de fases hasta recuperar su estructura, la 
composición de especies y su función ori- 
ginal (Martínez, 1996). En las estrategias 
de restauración de suelos degradados, es 
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determinante el uso de la microbiota del 
suelo (Haselwandter, 1997; Requena et 
al., 2001; Álvarez, 2009), en especial de 
hongos formadores de micorrizas, ya que 
previenen la erosión del suelo al ser un 
factor en la formación de agregados y en 
la acumulación de materia orgánica. De 
esta manera, le dan estructura al suelo y 
reducen o evitan el efecto de la erosión 
tanto hídrica como eólica, así como la 
consecuente pérdida de nutrientes. Ade- 
más, la asociación micorrícica otorga ven- 
tajas para el establecimiento de las 
plantas ante condiciones estresantes, 
como la habilidad de la planta hospedera 
para resistir a las condiciones de sequía 
(Haselwandter, 1997; Álvarez, 2009). La 
micorrización de varias especies de 
angiospermas y gimnospermas con hon- 
gos ectomicorrícicos ha tenido éxito, 
sobre todo en experimentos in vltro y, en 
menor cantidad, en condiciones de vivero. 
Entre los géneros más utilizados se 
encuentran Hebeloma, Lacearía, Rhlzo- 
phogon, Scleroderma, Suillus y Pisolithus 
(Pera y Parladé 2005). Por otra parte, la 
inoculación en vivero para plantas desti- 
nadas a la repoblación, se ha realizado 
con especies cosmopolitas, como Sclero- 
dema sp. y Pisolitus tinctorius (Valdés et 
al., 2010). El interés en utilizar Pisolithus 
tinctorius, se deriva de su éxito en la ino- 
culación de plantas forestales en países 
de los cinco continentes. Las plantas con 
las que P. tinctorius establece simbiosis 
ectomicorrícica incluyen más de veinte 
géneros de gimnospermas y angiosper- 
mas con distribución mundial; incluyendo 
especies forestales de las familias Casua- 
rinaceae, Dipterocarpaceae, Pinaceae, 
Myrtaceae y Oleaceae (Pérez-Moreno y 
Read, 2004; Ambriz etal., 2010). El cuerpo 
fructífero de este hongo mide entre 5 cm y 
20 cm de altura y hasta 10 cm de diáme- 
tro, fibroso de color pardo amarillento, 
esporas de color marrón canela, globosas 
de 7 p a 12 p y espinas de hasta 2 p de 
largo. Fructifica al comienzo del otoño. Un 
aspecto muy importante es que puede 


sobrevivir en suelos empobrecidos o 
terrenos perturbados, así como en suelos 
de alta acidez, o con altas concentracio- 
nes de metales pesados y resiste perio- 
dosdeestrésporsequía (García-Rodríguez 
etal., 2009). 


OBJETIVOS 

Evaluar las plantas de Plnus pseudostro- 
bus en un experimento en condiciones de 
restauración ecológica en presencia/ 
ausencia de inoculación del hongo ecto- 
micorrícico Pisolithus tinctorius. Comparar 
las plantas en condiciones de restaura- 
ción con las de una reforestación conven- 
cional, con la finalidad de identificar las 
estrategias que mejoren la supervivencia 
y desempeño de las plantas para su esta- 
blecimiento en un sitio severamente 
degradado y de esa forma recuperar una 
cobertura vegetal. 

MATERIALES Y MÉTODO 
Área de estudio 

El área de estudio se ubica en el paraje 
conocido como Huertitas, ejido de Até- 
cuaro, municipio de Morelia en Michoacán 
(Fig. la), lugar desprovisto de vegetación 
en donde se dan altas tasas de erosión en 
suelos acrisoles (acrisol ócrico), y existen 
numerosas cárcavas en al menos 6% de 
la superficie total. Presenta pendientes de 
18% en promedio, por lo que es muy sen- 
sible a la erosión hídrica. El resto del pai- 
saje está calculado en 46% de pastizales 
y 28% de tierras de cultivo (Duvert et al., 
2010). La erosión en cárcavas, como con- 
secuencia de la deforestación y técnicas 
de cultivo inadecuadas, han sido identifi- 
cadas en más de 65% de la superficie de 
los suelos del país (Covaleda, 2008), por 
lo que las consecuencias de la degrada- 
ción acelerada ha dado como resultado 
suelos sin capacidad de mantener la vida 
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vegetal. Un ejemplo claro se presenta en 
los suelos acrisoles; según la fao-isric 
(2006), estos suelos han sido clasificados 
como suelos evolucionados, que se han 
desarrollado principalmente sobre rocas 
muy ácidas o arcillas muy meteorizadas, 
que sufren posterior degradación. Se 
caracterizan por presentar un contenido 
mayor de arcillas en los horizontes más 
profundos como resultado de procesos 
pedogenéticos. Los acrisoles tienen, a 
ciertas profundidades, una baja satura- 
ción de bases y arcillas de baja actividad. 
Por lo general, estos suelos son arcillosos 
y tienden a compactarse cuando carecen 
de cobertura vegetal, impidiendo la infiltra- 
ción del agua (Covaleda, 2008). 

La zona presenta un clima subhú- 
medo, con lluvias en verano, temperatura 
media anual de 13,8 °C y precipitación 
media anual de 1000 mm con una des- 
carga media durante la época de secas de 
0,05 ls _1 -1 ,0 ls _1 y en la época de lluvias de 
10 Is 1 - 30 Is 1 , una altitud de 2150 msnm 
- 2450 msnm (Duvert et al., 2010). 


Establecimiento de pruebas de 
restauración 

Plantación sitio A 

En el año 2009 se estableció un experi- 
mento con fines de restauración con Pinus 
pseudostrobus. Las semillas procedentes 
de la región (Nuevo San Juan Parangari- 
cutiro) se estratificaron a 4 °C (como trata- 
miento pregerminativo) en condiciones 
controladas, fueron desinfectadas previa- 
mente con hipoclorito de sodio comercial 
al 20% (NaCIO 20 mi: 80 mi H 2 0), y poste- 
riormente se colocaron en cajas de Petri 
con papel filtro como sustrato humedecido 
con 2 mi de agua destilada. Las cajas fue- 
ron selladas con parafilm (R) y permane- 
cieron durante 15 días a una temperatura 
de 4 °C. Una vez transcurrido este tiempo, 
se trasladaron a una cámara de creci- 


miento a una temperatura constante de 
25 °C con 12 horas luz, manteniendo la 
humedad constante durante el proceso de 
germinación que fue de 21 días. Después 
de este periodo, se realizó el trasplante 
una vez que la radícula ya estaba emer- 
gida. Las plántulas fueron propagadas en 
una casa de sombra en contenedores rígi- 
dos de plástico medianos (350 cm 3 ), se 
utilizó un sustrato comercial a base de 
fibra de coco, corteza de árbol y agrolita 
mezclado con arena en combinación 1 : 1 y 
esterilizado. A las cuatro semanas de 
edad, la mitad de las plantas de pinos fue- 
ron inoculadas con 500 000 esporas del 
hongo ectomicorrícico Pisolithus tinctorius 
(el inoculo fue adquirido de manera 
comercial Biosyneterra Solutions Inc. 
L’Assomption, Quebec, Canadá) por cada 
planta. Formando dos tratamientos, uno 
de plantas inoculadas y otro de plantas sin 
inoculación, se consideraron 90 réplicas 
por cada tratamiento en un diseño com- 
pletamente al azar. Posteriormente, a los 
tres meses de edad, se realizó un tras- 
plante a bolsa de vivero de 600 cm 3 , con 
mayor cantidad de sustrato y se mantuvie- 
ron en casa de sombra durante 11 meses 
más, con riego constante cada tercer día. 
Posteriormente, las plantas se sacaron de 
la casa de sombra para su endurecimiento 
durante quince días antes de ser llevadas 
a campo. 

Una vez que inició la época de llu- 
vias, las plantas fueron transportadas a 
campo en cajas perforadas de plástico 
para facilitar el traslado hasta el sitio de 
estudio, donde permanecieron durante un 
día antes de ser trasplantadas. En la cepa 
de 40 cm 3 se adicionó !4 litro del medio de 
cultivo antes descrito y se realizó el tras- 
plante a los 15 meses de edad de las plan- 
tas, con el cuidado suficiente para que las 
raíces no quedaran dobladas o expuestas 
y se cubrió bien la cepa con sustrato del 
sitio en un diseño completamente al azar 
con 90 réplicas por tratamiento. 
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Plantación sitio B 

En el mismo año se localizó otra planta- 
ción en una cárcava aledaña muy cercana 
(177 m) con condiciones muy similares de 
pendientes, en el mismo tipo de sustrato y 
plantadas al mismo tiempo (Fig. Ib). La 
plantación con Pinus pseudostrobus se 
trataba de una reforestación convencio- 
nal, de las que se hacen anualmente en la 
región. Sin embargo, se desconoce la 


edad de las plantas al tiempo del tras- 
plante, así como las condiciones de alma- 
cenamiento de la semilla, germinación, 
propagación o endurecimiento de la planta 
previo a la reforestación, incluso, la proce- 
dencia de las semillas (mismas condicio- 
nes de los programas de reforestación 
anuales). Se evaluaron también 90 plan- 
tas para que fueran comparables con las 
del experimento de restauración. Se sabe 
que, dadas las condiciones de las planta- 



Figura la. Ubicación del área de estudio. Ib Sitios de estudio en la localidad de 
Atécuaro, municipio de Morelia, Michoacán. A) Sitio donde se encuentran los 
experimentos de restauración. B) Sitio donde se encuentra la reforestación convencional. 
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ciones en dos sitios diferentes, ocurre que 
si se consideraran todos los datos como 
un solo experimento, un tratamiento com- 
pleto (Sitio B) se encuentra en un sitio 
diferente, lo que corresponde a una pseu- 
dorréplica sensu Hurlbert (1984). A pesar 
de lo anterior se incluyeron los datos del 
sitio B por ser una comparación útil (como 
testigo), pero se debe reconocer que las 
diferencias podrían deberse a factores 
inherentes a los dos sitios y no a otras 
diferencias. 


Evaluaciones 

Se realizó la evaluación de altura total, 
diámetro a la altura de la base (dab) con 
un vernier tomada en centímetros, cober- 
tura en cm 2 obtenida a partir de la evalua- 
ción de dos diámetros de la copa e inferida 
mediante la fórmula de la elipse (r/^), 
número de ramas y supervivencia a las 
plantas al tercer año de edad de estableci- 
das las plantaciones de restauración en 
presencia/ausencia de inoculación con 
Pisolithus tinctorius (sitio A), y a las plan- 
tas de la reforestación convencional (sitio 
B). Los datos fueron analizados mediante 
Análisis de Varianza (anova) en el paquete 
estadístico R (R Development Core Team, 
2011). El análisis principal se realizó 
tomando en cuenta inicialmente sólo los 
datos del experimento del sitio A, al ser 
éste un diseño estadísticamente robusto 
y, de manera adicional, un segundo análi- 
sis, consistió en comparar una plantación 
de reforestación convencional con la de 
restauración. Para analizar la superviven- 
cia se siguió la misma estrategia, pero a 
través de análisis de modelos lineales 
generalizados (glm) para datos con distri- 
bución binomial. 


RESULTADOS 

Las plantas se desempeñaron de manera 
diferente dependiendo de la forma en que 


fueron propagadas y trasplantadas en el 
sitio. Al comparar las plantas de la refores- 
tación en el sitio A, en función de la presen- 
cia/ausencia de inoculación con el hongo 
ectomicorrízico Pisolithus tinctorius, no se 
obtuvieron diferencias significativas para 
variables de crecimiento (altura F (1 267) = 
0,27, P = 0,604 , diámetro a la altura de la 
base F (1267) = 2,45, P = 0,119, cobertura 

F (1,267) = °’ 782 ’ P = °’ 779 y ramaS F ( 1 , 267 ) = 

2,42, P = 0, 12). Sin embargo, para la super- 
vivencia sí se presenta un incremento sig- 
nificativo en las plantas inoculadas. Por 
otro lado, al hacer la comparación de las 
plantaciones de los sitios A y B se presen- 
taron diferencias estadísticamente signifi- 
cativas en todas las variables de 
crecimiento y en supervivencia. 

En términos de la altura, cuando se 
analizan sólo los tratamientos del experi- 
mento de inoculación no hay diferencias 
significativas, pero cuando se incluyen los 
datos de las plantas de la reforestación 
convencional (sitio B), sí se presentan 
diferencias estadísticamente significativas 
( F (2,267) = 87 ’ 87 ’ p < 0,0001). La altura de 
los tratamientos de restauración en pre- 
sencia/ausencia de hongos micorrícicos 
(80 cm), contrasta con las plantas de la 
reforestación convencional (40 cm) (Fig. 

2) . Para el dab se presentó un patrón simi- 
lar, entre tratamientos de inoculación no 
se presentaron diferencias, pero sí entre 
las plantas del sitio A y las del sitio B 
[P (2 267 ) ~ 233,75, P < 0,00001). Los valores 
de las plantas de la restauración son de 
alrededor de 3 cm, mientras que para las 
de la reforestación es menor a 1 cm (Fig. 

3) . Para la cobertura, se presentan dife- 
rencias altamente significativas entre los 
tratamientos del experimento de inocula- 
ción con los de la reforestación en el sitio 
B ( F (2267) = 147,07, P < 0,0001), con ten- 
dencias muy similares a las de la altura, 
con valores similares entre las plantas de 
los tratamientos de restauración (40 cm 2 ) 
pero muy diferentes a las plantas de la 
reforestación (10 cm 2 ) (Fig. 4). Para el 


Cobertura (cm 2 ) 
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Figura 2. Altura de Pinus 
pseudostrobus en condiciones de 
restauración ecológica con inoculación 
de hongo ectomicorrícico ( Pisolithus 
tinctorius) (he), restauración sin 
inoculación y reforestación 
convencional. Las barras muestran el 
error estándar y las letras grupos de 
acuerdo con la prueba de Tukey. 


Figura 3. Diámetro de Pinus 
pseudostrobus en condiciones de 
restauración ecológica con 
inoculación de hongo ectomicorrícico 
(Pisolithus tinctorius) (he), 
restauración sin inoculación y 
reforestación convencional. Las 
barras muestran el error estándar las 
letras grupos de acuerdo con la 
prueba de Tukey. 
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Figura 4. Cobertura de Pinus 
pseudostrobus en condiciones de 
restauración ecológica con inoculación 
de hongo ectomicorrícico ( Pisolithus 
tinctorius) (HE), restauración sin 
inoculación y reforestación 
convencional. Las barras muestran el 
error estándar las letras grupos de 
acuerdo con la prueba de Tukey. 
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número de ramas, se presentan tenden- 
cias muy similares, con valores mayores 
en las plantas del experimento en el sitio 
A en ausencia de hongos micorrícicos, no 
siendo diferente de los que se encuentran 
en presencia de hongos micorrícicos, 
pero sí estos dos de las plantas de la refo- 
restación en el sitio B (P (2267) = 232, 78, P < 
0,0001 ) (Fig. 5). 

El resultado más interesante se da 
en términos de los porcentajes de supervi- 


vencia; para el experimento de inocula- 
ción, las plantas inoculadas presentaron 
86% de supervivencia, en ausencia de 
inoculo 62%, siendo una diferencia en 
supervivencia estadísticamente significa- 
tiva (X 2 = 13,053, g.l=1; P=0,0003). Si se 
incluye en el análisis a las plantas presen- 
tes en el sitio B correspondientes a las de 
una reforestación convencional, éstas 
presentaron 30% de supervivencia al ter- 
cer año del trasplante. Al comparar con los 
tratamientos de inoculación, existe una 



Figura 5. Número de ramas de Pinus 
pseudostrobus en condiciones de 
restauración ecológica con inoculación 
de hongo ectomicorrícico ( Pisolithus 
tinctorius) (he), restauración sin 
inoculación y reforestación 
convencional. Las barras muestran el 
error estándar las letras grupos de 
acuerdo con la prueba de Tukey. 


Figura 6. Supervivencia de Pinus 
pseudostrobus en condiciones de 
restauración ecológica con 
inoculación de hongo ectomicorrícico 
(Pisolithus tinctorius) (he), 
restauración sin inoculación y 
reforestación convencional. Letras 
diferentes indican diferencias 
significativas con p<0,05. 
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diferencia altamente significativa entre los 
tres (X 2 =61 ,366, g.l=2; P<0,0001 ; Fig. 6). 

DISCUSIÓN 

Los resultados del experimento sugieren 
que, en condiciones de restauración, con 
estrés nutricional e hídrico, en cuanto a 
las variables estudiadas, las plantas ino- 
culadas con el hongo ectomicorrícico 
Pisolithus tinctorius, no responden de 
manera diferente que las plantas sin inoculo. 
Sin embargo, sí les confiere beneficios 
para la adquisición de nutrientes que les 
provee de mayor supervivencia (86%), lo 
cual implica una gran ventaja para el 
establecimiento de las plantas ante estas 
condiciones. Las plantas en ausencia del 
inoculo presentaron 62%. Por otro lado, 
aun cuando la asociación simbiótica 
puede representar un gasto energético 
para los pinos, debido a que tienen que 
intercambiar hidratos de carbono (Choi et 
al., 2005), la relación es positiva, ya que 
los hongos micorrícicos proveen de 
nutrientes y agua a la planta, además de 
que le ayudan a amortiguar cambios 
bruscos de temperatura (Smith y Read, 
2000 ). 

La disponibilidad de agua tiene el 
efecto directo de disminuir la turgencia, e 
indirectamente limita el incremento de bio- 
masa, por lo que una cantidad adecuada 
de agua disponible es esencial para man- 
tenerla productividad de las plantas (Mon- 
toya, 2000), de modo que la presencia de 
micorrizas puede influir de manera impor- 
tante para tener disponibilidad de agua e 
incrementar la biomasa, sobre todo en la 
parte radical. El sustrato, rico en materia 
orgánica que se adicionó a las plantas del 
experimento de restauración, muy proba- 
blemente también contribuyó a generar 
los patrones observados, ya que a las 
plantas de la reforestación convencional 
no se les suministró ningún sustrato adi- 
cional y la materia orgánica del suelo con- 


tribuye a un mejor desempeño de las 
micorrizas. En estas condiciones de refo- 
restación convencional las plantas pre- 
sentaron 30% de supervivencia. 

Arias (1992), menciona que P. pseu- 
dostrobus plantado bajo buenas condicio- 
nes de fertilidad, presenta alta tasa de 
crecimiento, pero también se ha publicado 
que no puede prosperar en sitios con con- 
diciones edáficas adversas (Domínguez- 
Caballero etal., 2001). El presente estudio 
sugiere que aplicando estrategias ade- 
cuadas, esta especie puede mostrar un 
buen desempeño temprano, aun sobre 
cárcavas. Pera y Parladé (2005) sugieren 
que la utilización de plantas micorrizadas 
pueden jugar un papel importante en la 
reforestación de zonas fuertemente alte- 
radas. La inoculación con hongos micorrí- 
cicos puede incrementarsignificativamente 
la supervivencia de Pinus pseudostrobus 
como se muestra en este experimento. 

Valdés et al. (2010) mencionan que 
la inoculación con Pisolithus tinctorius 
(cepa 202) a los 18 meses puede incre- 
mentar significativamente el volumen y 
biomasa (2,93 cm 3 y 3,68 g) de P. pseu- 
dostrobus con respecto al control (2,3 cm 3 
y 3,39 g). Por otro lado, se han realizado 
experimentos con esta especie de hongo 
ectomicorrícico (P. tinctorius) inoculando 
plantas de Clstus albldus en condiciones 
de campo y mencionan que la inoculación 
tiene un efecto positivo sobre el creci- 
miento de las plantas, obteniendo hasta 
30% mayor crecimiento que en las plantas 
control (Alguacil et al., 2003). 

En experimentos con Pinus greggll 
inoculados con los hongos ectomicorríci- 
cos Lacearla laccata, Sulllus pseudobrevl- 
pes, Boletus plnophllus y Pisolithus 
tlncorlus, Carrera-Nieva y López-Ríos 
(2004) encontraron diferencias estadísti- 
camente significativas a los nueve meses 
de inoculación en la altura y biomasa en 
peso seco de la planta, independiente- 
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mente del hongo que utilizaron, para la 
variable diámetro, el hongo que presentó 
un mayor efecto fue B. pinophilus. En el 
mismo estudio, para la especie Pinus 
patula utilizaron las especies Lacearía 
laccata, Suillus pseudobrevipes, Boletus 
clavipes, Amanita fu Iva y Amanita caesa- 
rea, encontrando que L. laccata es el 
hongo que presenta un mayor efecto en 
las variables altura, diámetro y biomasa 
en peso seco. Por otro lado, en estudios 
con Pinus ponderosa inoculados con los 
hongos micorrícicos Suillus luteus, Hebe- 
loma mesophaeum y Rhizopogon 
roseolus, Barroetaveña et al. (2012) no 
encontraron diferencias estadísticamente 
significativas en el diámetro ni en el largo 
del vástago, comparativamente con las 
plantas control. 

La respuesta puede deberse tanto a 
la especie del hongo ectomicorrícico, 
como a la especie de la planta hospedera 
debido al reconocimiento en la asociación. 

Debido a que ambos sitios se 
encuentran a una corta distancia (177 m 
en línea recta), los resultados sugieren 
que las diferencias entre las plantas del 
sitio A, el experimento de inoculación y del 
sitio B se deben a los tratamientos aplica- 
dos, aun cuando es posible que se deban 
a diferencias microambientales como la 
temperatura, la luz disponible o la hume- 
dad (Prince et al., 2001), que no fue posi- 
ble medir en el presente estudio. 

Tomar en cuenta todos los aspectos 
que se mencionan tanto para la germina- 
ción como para la propagación, implica 
una inversión adicional de recursos eco- 
nómicos que, sin embargo, se compensa 
por la mejor supervivencia y desempeño 
general de la plantación y que incluso 
puede ser calculado (García-Frapolli y 
Lindig-Cisneros, 2011). 

Si no se cuenta con recursos adicio- 
nales, se puede reducir el área refores- 


tada y asignar los recursos ahorrados a 
estas medidas adicionales; esta estrate- 
gia incluso permitiría contar con una 
mayor área reforestada en el largo plazo 
(Sáenz-Romero y Lindig-Cisneros, 2004). 


CONCLUSIONES 

La inoculación con el hongo ectomicorrí- 
cico Pisolithus tinctorius incrementa la 
supervivencia de Pinus pseudostrobus en 
sitios desprovistos de vegetación, con 
numerosas cárcavas y graves problemas 
de erosión. 

Por otro lado, el conjunto de condi- 
ciones a las que son sometidas las plan- 
tas durante la germinación (tratamiento 
pregerminativo), propagación (sustrato, 
endurecimiento) y trasplante (edad de las 
plantas, fecha y método de plantación) 
son de vital importancia para lograr el 
éxito en su establecimiento, supervivencia 
y crecimiento. El desempeño de las plan- 
tas (altura, diámetro, cobertura y número 
de ramas) independientemente de si las 
plantas fueron inoculadas o no, fue supe- 
rior a las de la reforestación convencional. 

Lo más sobresaliente del presente 
ensayo fue que las plantas inoculadas, 
presentaron un porcentaje de superviven- 
cia mayor, seguidas de las plantas no ino- 
culadas pero propagadas bajo las mismas 
condiciones que las primeras y muy por 
debajo los pinos de la reforestación con- 
vencional. Pinus pseudostrobus inoculado 
con Pisolithus tinctorius tiene potencial 
para reforestar sitios severamente degra- 
dados con problemas de erosión y forma- 
ción de cárcavas, como es el caso en este 
sitio de estudio. El conjunto de estrate- 
gias representan un gran potencial de 
restauración ecológica y, si las tasas de 
crecimiento a mediano plazo fueran 
satisfactorias, también para el estableci- 
miento de plantaciones con fines de pro- 
ducción. 
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ARTÍCULO DE INVESTIGACIÓN 


Tableros de partículas fabricados 
con residuos industriales de madera 

de Pinus patula 

Partióle boards manufactured from industrial residues 

of Pinus patula wood 

Luis Zambrano , 1 Pablo Moreno, 2 * Freddy Muñoz , 3 Jorge Duran , 2 

Darío Garay 2 y Styles Valero 2 

RESUMEN 

Este estudio es una propuesta para dar un aprovechamiento y uso a los residuos industriales del 
aserrío de la madera de Pinus patula para la elaboración de tableros aglomerados de partículas. Por lo 
tanto, los objetivos fueron manufacturar tableros aglomerados de partículas homogéneos a dos niveles 
de densidad, determinar sus propiedades físico-mecánicas y evaluar los resultados bajo las normas 
alemana y venezolana. La absorción de agua promedio a 2 y 24 horas de los tableros a ambas densi- 
dades cumplió con lo estipulado en las normas. La variación del espesor (hinchamiento) promedio por 
inmersión en agua a 2 y 24 horas para los tableros con densidad de 767 kg nr 3 cumplió con lo exigido 
en ambas normas. Para el caso de los tableros con densidad de 672 kg nr 3 la variación de espesor a 
24 horas cumplió con las normas, mientras que la variación del espesor a 2 horas no. Respecto de las 
propiedades mecánicas, tanto el módulo de ruptura en flexión estática como la tensión perpendicular 
(adhesión interna) para ambos niveles de densidad superaron el valor mínimo aceptado en las normas; 
demostrando que los residuos industriales de la madera de P. patula presentan buenas características 
tecnológicas para la fabricación de tableros de partículas de densidad media. 

PALABRAS CLAVE 

Flexión estática, densidad, propiedades físico-mecánicas, tableros de partículas de densidad media, 
tensión perpendicular, urea-formaldehído. 


ABSTRACT 

This study is a proposal to generate uses of industrial residues from sawmills of Pinus patula 
wood for the manufacturing of particleboards. The aims are to manufacture homogeneous particle- 
boards with two density levels, determine their physical and mechanical properties and evalúate the 
results with Germán and Venezuelan standards. The average water absorption at 2 and 24 hours of 
both density levels of boards meets the requirements for both standards. The average thickness varia- 
tion (swelling) by immersion en water at 2 and 24 hours for boards with average density of 767 kg nr 3 
meets the requirements for both standards. For the case of the panels with average density of 672 kg 
nr 3 thickness variation at 24 hours meets the requirements of both standards, while the swelling at 2 
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hours does not. Regarding the mechani- 
cal properties both the modulus of rup- 
ture in static bending and tensión 
perpendicular (internal bond) for both real 
density levels exceed the minimum 
accepted by standards, demonstrating 
that industrial residues from P. patula 
wood exhibit good technological charac- 
teristics for médium density particle- 
boards manufacture. 

KEYWORDS 

Static bending, density, physical and 
mechanical properties, médium density 
particleboard, perpendicular tensión, 
urea-formaldehyde. 


INTRODUCCIÓN 

Para un buen desarrollo de la industria 
forestal y maderera venezolana es impres- 
cindible el aprovechamiento integral de la 
madera, optimizando el rendimiento de los 
recursos maderables, es decir, utilizando 
en los procesos productivos el mayor por- 
centaje posible de la madera cosechada 
del bosque. Tradicionalmente, en Vene- 
zuela se aprovechan solamente las trozas 
provenientes del fuste principal hasta una 
altura comercial aproximada de 15 m, lo 
cual implica la pérdida de un volumen con- 
siderable de madera que incluye trozas de 
altura superior, copa y ramas. De igual 
manera, los residuos producidos por el 
procesamiento primario en las industrias 
madereras no son aprovechados, gene- 
rando un volumen considerable de 
subproductos de material lignocelulósico, 
que particularmente en Venezuela no han 
sido aprovechados, debido a la existencia 
de recursos maderables abundantes y 
falta de políticas de aprovechamiento inte- 
gral de los subproductos generados del 
procesamiento forestal industrial. Sin 
embargo, estos subproductos pueden ser 
utilizados en la fabricación de otros pro- 
ductos de ingeniería con múltiples aplica- 
ciones. 

Según lo señalan Alma et al., (2004), 
la población del mundo consume actual- 


mente más de 3,5 millones de toneladas 
de madera verde al año. Si la tasa de con- 
sumo de fibra de madera y la tasa de 
crecimiento de la población (de aproxima- 
damente 90 millones de personas por 
año) se mantienen constantes, la 
demanda de fibra de madera se incremen- 
tará en más de 60 millones de toneladas 
cada año, lo que podría aumentar la tasa 
de deforestación (Zheng et al., 2006), pro- 
piciando el agotamiento de algunas 
especies forestales de bosque natural. 
Aprovechar de manera óptima e integral el 
recurso forestal, los residuos industriales, 
las especies de rápido crecimiento y las 
fibras no leñosas, como los residuos 
agrícolas (bagazo y paja de trigo), es una 
buena opción para mitigar y disminuir pre- 
siones políticas y sociales, sobre la explo- 
tación de maderas procedentes de bosque 
natural (Latibari y Roohnia, 2010). 

Un producto de ingeniería que se 
adapta a los requerimientos anteriormente 
mencionados son los tableros de partícu- 
las, pues el suministro de materia prima 
para su elaboración potencia el uso de 
subproductos originados por las operacio- 
nes forestales (podas y raleos) hasta los 
del proceso primario y secundario de la 
madera (despuntes, descuadres, costane- 
ras, corteza, entre otros). De esta manera 
se generan nuevas alternativas para pro- 
mover el desarrollo de la industria made- 
rera venezolana con nuevos productos de 
ingeniería con diferentes aplicaciones. 

Existen diversos antecedentes biblio- 
gráficos referidos al uso exitoso de resi- 
duos provenientes de la industria forestal 
y maderera en la elaboración de tableros 
aglomerados de partículas, como ejemplo 
están las investigaciones realizadas por 
(Ginzel y Peraza, 1966; Moslemi, 1974; 
Peredo y Deppe, 1984; Suchsland y Xu, 
1989; Poblete y Sánchez, 1991; Peredo y 
Torres, 1991; Van Niekerk y Pizzi, 1994; 
Blanchet et al., 2000; Moreno et al., 2005; 
Pedieu et al., 2008; Latibari y Roohnia, 
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2010; Bazyar et al., 2011). Por otro lado, 
Garay et al. (2009), estiman que las pro- 
piedades de un tablero están determina- 
das por los siguientes factores: tipo, 
tamaño, distribución y orientación de las 
partículas, así como también el tipo y la 
calidad del proceso, contenido de hume- 
dad durante el prensado y tratamiento 
post-manufactura del tablero. 


con dimensiones de 30 mm x 30 mm x 1 00 
mm para la determinación de la densidad 
en la condición verde, seca al aire y seca 
al horno, según la norma din (1976). El 
cálculo del contenido de humedad se rea- 
lizó con la siguiente ecuación: 


%CH = 


( P -P\ 

i c 

P 


100 


( 1 ) 


El presente estudio es una propuesta 
sobre el aprovechamiento y uso de los 
residuos industriales del aserrío de la 
madera de P. patula para la elaboración 
de tableros aglomerados de partículas. 


OBJETIVOS 

Los objetivos fueron: a) manufacturar 
tableros aglomerados de partículas homo- 
géneos a dos niveles de densidad, b) 
determinar sus propiedades físico-mecáni- 
cas y c) evaluar los resultados de acuerdo 
con lo establecido en la norma alemana din 
(1981) y venezolana covenin (1991). 

MATERIALES Y MÉTODOS 

Como materia prima se utilizaron costane- 
ras derivadas del proceso de aserrío de 
árboles de P. patula de 25 años de edad, 
procedentes del bosque La Carbonera- 
San Eusebio, Mérida-Venezuela. La ela- 
boración de los tableros, la preparación de 
probetas y la ejecución de los ensayos, se 
realizaron en el Laboratorio Nacional de 
Productos Forestales: Sección tableros 
aglomerados de la Facultad de Ciencias 
Forestales y Ambientales de la Universi- 
dad de Los Andes, Venezuela. 


Determinación de la densidad y contenido 
de humedad de la materia prima 

De los segmentos de costaneras seleccio- 
nadas al azar, se prepararon 30 probetas 


donde: 

% CFI = contenido de humedad. 

p= peso inicial de la muestra (g). 

p 0 = peso de la muestra secada al horno 

(g). 

Viruteado de los residuos de P. patula y 
secado de partículas 

Las costaneras fueron procesadas en una 
viruteadora de cuchillos para obtener par- 
tículas con una dimensión promedio de 20 
mm de longitud, ancho variable y espesor 
de 0,25 mm. Desde la viruteadora las par- 
tículas fueron transferidas por un sistema 
neumático de ciclones a la secadora de 
suspensión Keller, donde fueron secadas 
a un contenido de humedad final de 4%. 


Coeficiente de esbeltez 

De manera aleatoria se tomaron 30 g de 
las partículas secas para determinar su 
geometría y tamaño, éstas se clasificaron 
de manera visual en tres categorías de 
tamaño: grandes, medianas y pequeñas. 
Cada categoría fue pesada y se tomaron 
30 muestras de cada una. Con un vernier 
digital de 0,02 mm (0,0001 pulgadas) de 
resolución, se midió largo, ancho y espe- 
sor de las partículas y se determinó su 
coeficiente de esbeltez (E), con la ecua- 
ción (2). La proporción de participación de 
cada categoría de tamaño de partícula 
para conformar el tablero, se determinó 
mediante su peso. 
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donde: 

E= coeficiente de esbeltez. 

L= longitud de las partículas (mm). 
e= espesor de las partículas (mm). 


Componentes para la fabricación de 
los tableros 

Se utilizaron los siguientes insumos: partí- 
culas de P. patula a 4% de contenido de 
humedad, adhesivo urea formaldehído 
(uf) a 60% de sólidos con viscosidad = 
250 cP, pH = 8,5 y como catalizador una 
solución de sulfato de amonio al 30% (p/p) 
de concentración, la cantidad de cataliza- 
dor corresponde al 8% (p/p) de la cantidad 
de adhesivo aplicada a cada tablero. Se 
fabricaron cinco tableros homogéneos 
para cada nivel de densidad. 


Encolado de las partículas 

El encolado de las partículas secas a 4% 
se realizó por el método de aspersión a 
0,203 MPa (2 atm) por un tiempo de 8 
minutos en una encoladora Drais usando 
una carga de adhesivo UF catalizado de 12 % 
respecto al peso seco de las partículas. 


Formación del colchón, proceso de 
pre-prensado y prensado del tablero de 
partículas 

Luego de encoladas las partículas fueron 
distribuidas, manualmente, sobre una 
lámina de aluminio previamente parafi- 
nada y humectada con agua. Para rete- 
ner las partículas en todos los costados 
del colchón se utilizó un marco o forma- 
leta de madera. Una vez conformado el 
colchón de partículas se aplicó una carga 


de aproximadamente 785 N (80 kg) a 
temperatura ambiente con la finalidad de 
compactar y disminuir el espesor, evi- 
tando la deformación del tablero en el 
proceso de prensado. Finalmente, para 
obtener los dos niveles de densidad en 
los tableros de partículas se aplicaron los 
respectivos ciclos de prensado (Fig. 1) en 
una prensa Becker van Huller bajo las 
siguientes condiciones operacionales: 
presión específica sobre el tablero: 3,05 
MPa (31,14 kg crrr 2 ), presión manomé- 
trica de la prensa simple: 14,71 MPa (150 
kg crrr 2 ), temperatura de los platos de la 
prensa: 180 °C a 200 °C, tiempo de pren- 
sado total: 9 min para densidad teórica 
650 kg nr 3 y 18 min para densidad teórica 
750 kg rrr 3 . 


Acondicionamiento y dimensionado de 
los tableros 

Los tableros fabricados fueron almacena- 
dos en un cuarto de acondicionamiento a 
20 °C ± 1 °C y una humedad relativa de 
65% ± 2%. Los tableros permanecieron 
bajo esas condiciones hasta que alcanza- 
ron un contenido de humedad de equili- 
brio de 12% aproximadamente. Luego los 
tableros fueron dimensionados en una 
sierra circular con escuadría de 500 mm x 
500 mm, conservando el espesor de pren- 
sado de 19 mm. 


Ensayos físicos y mecánicos de los 
tableros 

A partir de los tableros dimensionados se 
obtuvieron 10 probetas para cada uno de 
los ensayos, según las normas din (1 965 a - 
b c dye ). Los ensayos físicos fueron: con- 
tenido de humedad, densidad del tablero, 
absorción por inmersión en agua y varia- 
ción de espesor por inmersión en agua a 
2 y 24 horas. Los ensayos mecánicos 
fueron: resistencia a flexión estática 
(módulo de ruptura, mor) y tensión per- 
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— t m Densidad teórica 650 kg.rrr 3 -i- Densidad teórica 750 kg.rrr 3 



Tiempo de prensa (min) 


Figura 1. Diagrama del ciclo de prensado aplicado para densificar los tableros 
homogéneos de partículas de Pinuspatula a las densidades teóricas propuestas. 


pendicular a la superficie del tablero 
(adhesión interna) usando una máquina 
universal Instron y una velocidad de 
ensayo de 5 mm/min para ambos ensa- 
yos. La evaluación de las propiedades 
físico-mecánicas de los tableros se rea- 
lizó según lo estipulado por la norma 
venezolana covenin (1991) y alemana 
din (1981), utilizadas para tableros de 
partículas de madera de uso general con 
prensado plano. 

Diseño experimental 

Se utilizó un diseño con un sólo factor 
(densidad) con dos niveles categóricos y 
seis variables respuesta (absorción de 
agua a 2 h y 24 h, variación de espesor 
o hinchamiento a 2 h y 24 h, módulo de 
ruptura en flexión estática y tensión 
paralela a la superficie del tablero). 

El análisis de varianza (andeva) se 
realizó con un nivel de confianza «= o, 05 
y con el software estadístico Design 
Expert® V 8.0. 


RESULTADOS 

Análisis de varianza (andeva) de un 
solo factor 

El andeva (Tabla 1) indica que el factor 
densidad del tablero es significativo para 
todas las variables respuesta. La probabi- 
lidad normal de residuales de los datos 
obtenidos en todas las variables respuesta 
cumple con el supuesto de normalidad, no 
se observaron fluctuaciones significativas 
de las muestras ensayadas ni datos atípi- 
cos. La proporción de variabilidad de cada 
parámetro respuesta es explicada por el 
modelo estadístico mediante el valor de 
R 2 . Por otra parte, el valor R 2 ajustado 
indica que el número de factores presen- 
tes en el modelo es adecuado, por lo 
tanto, ambos confirman la validez del 

ANDEVA. 

En la figura 2 se modela el efecto del 
factor densidad sobre las propiedades 
físico-mecánicas determinadas para los 
tableros aglomerados de P. patula. 
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Tabla 1 . Valor-p para las variables respuesta. 





Valor-p* 

(a = 0,05) 



Factor 

ABS 2b 

ABS 24b 

VE 2h 

VE 24h 

MOR 

TPST 

Modelo 

< 0,0001 

< 0,0001 

0,0002 

0,0022 

< 0,0001 

0,0002 

Densidad 

< 0,0001 

< 0,0001 

0,0002 

0,0022 

< 0,0001 

0,0002 

R2 

0,9964 

0,9881 

0,8371 

0,7103 

0,9768 

0,8347 

R2 ajustado 

0,9959 

0,9867 

0,8167 

0,6741 

0,9739 

0,8140 


Fuente: software estadístico Design Expert® V 8.0 

ABS= absorción de agua a 2 y 24 h; VE= variación de espesor a 2 y 24 h; MOR= módulo de ruptura en flexión está- 
tica; TPST= tensión perpendicular a la superficie del tablero 

* Valor-p < 0,05 indica que los términos del modelo son significativos. 

* Valor-p > 0,05 indica que los términos del modelo no son significativos. 
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Figura 2. Modelamiento con software Design Expert® del efecto del factor densidad 
sobre las propiedades físico-mecánicas de los tableros homogéneos de Pinus patula. 
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Densidad de la madera de P. patula y 
tableros de partículas 

Los valores de la densidad promedio obte- 
nidos para madera sólida utilizada como 
materia prima para la fabricación de table- 
ros aglomerados de partículas, se presen- 
tan en la tabla 2. Los intervalos de densidad 
y coeficientes de variación para los table- 
ros se muestran en la tabla 3. 

Una comparación entre la densidad 
teórica y la real de los tableros de partícu- 
las y madera sólida de P. patula se pre- 
senta en la figura 3. 


Coeficiente de esbeltez 

Los valores promedio de participación de 
cada intervalo de partículas de P. patula 
seleccionado y su respectivo coeficiente 
de esbeltez se muestran en la tabla 4. 

El valor promedio del coeficiente de 
esbeltez para las partículas de P. patula 
(148,92) se encuentra dentro del intervalo 
de esbeltez comprendido entre 120 y 200, 
intervalo que es considerado por Moslemi 
(1974), Medina (1975), Peredo (1988), 
Garay (1997) y Moreno etal. (2002), como 
óptimo para la elaboración de tableros de 
partículas con propiedades físicas y 
mecánicas aceptables. 


Tabla 2. Densidad promedio para la madera sólida de Pinus patula utilizada para la 
fabricación de tableros aglomerados de partículas. 


Condición 

Densidad promedio (kg m 3 ) 

Coeficiente de 
variación (%) 

Contenido de 
humedad (%) 

Verde 

684 

4,4 

30,50 

Seca al aire 

668 

6,2 

13,31 

Seca al horno 

618 

5,3 



Tabla 3. Densidad promedio y coeficiente de variación de los tableros aglomerados de 

partículas de Pinus patula. 


Densidad real promedio 

Coeficiente de variación 

(kg m-3) 

(%) 

672 

3,25 

767 

2,88 


72 


Tableros de partículas fabricados con residuos Industriales 


850 


6 800 


e 



700 


i 


650 


600 


• densidad real 1 

• densidad teórica 1 

• densidad real 2 

• densidad teórica 2 
densidad seca al aire P. patula 




Figura 3. Comparación entre densidad teórica y real obtenida para los tableros 
homogéneos de partículas y la densidad de la madera sólida seca al aire de Pinus 

patula. 


Tabla 4. Coeficiente de esbeltez (E) de las partículas de Pinus patula. 


Rango de tamaño en ancho 
(mm) 

Proporción en 
peso (%) 

Largo 

(mm) 

Espesor 

(mm) 

E 

> 18 

60,00 

41,07 

0,23 

178,56 

9-18 

27,50 

29,62 

0,22 

134,63 

< 9 

12,50 

29,39 

0,22 

133,59 

Promedio 

100 

26,69 

0,22 

148,92 


Absorción de agua por inmersión 

En las figuras 4a y 4b, se presenta el 
comportamiento de la absorción de agua, 
por inmersión a 2 y 24 horas, de los 
tableros a ambas densidades y de 
acuerdo con la norma covenin (1991) 


que solicita 25% y 60% de absorción para 
2 horas y 24 horas, respectivamente. Los 
tableros con densidad de 672 kg nr 3 y 
absorción de agua a 2 y 24 horas no cum- 
plen con lo solicitado por la norma. Mien- 
tras que, los tableros con densidad de 
767 kg nr 3 sí lo cumplen. 
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Variación de espesor por inmersión en 
agua 

En las figuras 4c y 4d, se presenta el 
comportamiento de la variación de espe- 
sor por inmersión en agua a 2 y 24 horas 
de los tableros con ambas densidades y 
de acuerdo con las normas covenin 
(1991) y din (1981) que solicitan una 
variación no mayor a 6% y 15 % a 2 horas 
y 24 horas, respectivamente. Los tableros 
con densidad de 672 kg nr 3 no cumplen 
con la exigencia de la norma covenin 
(1991) para la variación de espesor por 
inmersión en agua a 2 horas, pero sí cum- 
plen los requerimientos para la variación 
de espesor a 24 horas. 

» densidad 672 k 9 . n 1 3 
• densidad 767 kg'm3 
-•-COVENIN 847-91 

70 n 



< 

i* 10 


0 -I 1 I « 

1 2 3 4 5 

N 5 labler o 


• densidad 672 kgim3 
densidad 767 kgim3 
-•-COVENIN 847-91 / DIN 68761 
0 1 



** 1 

o -I , , , 

1 2 3 4 5 

N' tablero 


Los tableros con densidad de 767 
kg nr 3 sí cumplen con lo solicitado por las 
normas para la variación de espesor por 
inmersión en agua a 2 y 24 horas. 


Flexión estática (mor) 

En la figura 5a, se presentan los valores 
del mor en flexión estática de los tableros 
de P. patula para ambas densidades. Los 
tableros con densidad de 672 kg nr 3 y 767 
kg nr 3 superan en 23% y 52 %, respectiva- 
mente, al valor mínimo de 17,16 MPa (175 
kg cnr 2 ) establecido por las normas cove- 
nin (1991) y din (1981). 
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Figura 4. Absorción de agua por inmersión (a) 2 h (b) 24 h y variación de espesor o 
hinchamiento por inmersión en agua (c) 2 h (d) 24 h. Para tableros homogéneos de 
partículas de Pinus patula. De acuerdo a las normas covenin (1991) y din (1981). 
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Tensión perpendicular a la superficie 
del tablero 

En la figura 5b se presentan los valores 
obtenidos para la tensión perpendicular a 
la superficie de los tableros de P. patula 
para ambas densidades. Los tableros con 
densidad de 672 kg nr 3 y 767 kg nr 3 supe- 
ran en 19 % y 30 %, respectivamente, al 
valor mínimo de 0,34 MPa (3,5 kg cnr 1 2 ) 
solicitado por las normas covenin (1 991 ) y 
din (1981). 


DISCUSIÓN 

El ciclo de prensado permitió obtener 
tableros con valores similares a la densi- 
dad de la madera seca al aire de P. patula. 
Sin embargo, aun cuando en la conforma- 
ción del producto participe una sola espe- 
cie, la variabilidad de la densidad dentro y 
entre tableros es posible, situación que 
menciona en su investigación Moreno et 
al. (2005) con tableros de partículas de 
ramas de Fraxinus americana. No 
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obstante, los coeficientes de variación 
determinados para la densidad de los 
tableros de P. patula indican una buena 
distribución de las partículas en el pro- 
ducto final, de acuerdo con lo propuesto 
por Durán (1981). 

La densidad del tablero es una 
medida de la compactación individual de 
una partícula dentro del tablero y es 
dependiente, principalmente, de la densi- 
dad de la madera y la presión aplicada 
durante el prensado (Rocha et al., 1974). 
Un incremento en la densidad del tablero 
resulta en una mejor eficiencia de la 
resina, debido a la mejora de uniones 
adhesivas entre partícula y adhesivo (Gat- 
chell et al., 1966). Por ejemplo, Poblete y 
Peredo (1990) lograron establecer, para 
tableros de partículas con especies nati- 
vas chilenas, ecuaciones de regresión 
lineal entre la propiedad de flexión y den- 
sidad, logrando producir tableros de partí- 
culas con densidad de 650 kgnr 3 que 
cumplieron con el valor mínimo de 15,99 
MPa (163 kg cnr 2 ) solicitado por la norma 
din (1965 c ). 

De acuerdo con Rocha et al. (1974), 
la densidad de la madera de la especie 
influye en el consumo de adhesivo y la 
rugosidad superficial del tablero. El usar 
especies que tengan una alta densidad 
provoca una disminución de la resistencia 
mecánica del tablero ya que habría poca 
tasa de compresibilidad de las partículas, 
provocando el aumento de la expansión 
lineal e hinchamiento del espesor en el 
tablero. 

Independiente del tiempo de inmer- 
sión, el incremento de la densidad del 
tablero produce una disminución en la 
absorción de agua, comportamiento que 
señalan (Rocha et al., 1974; Durán, 1982; 
Poblete y Peredo, 1990; Moreno et al., 
2002) para tableros de partículas. Esta 
tendencia se debe a que una mayor den- 
sificación produce una disminución en la 


porosidad del tablero, menor cantidad de 
espacio vacío entre partículas y adhesivo, 
lo que impide la penetración del agua. 
Así lo demuestran las investigaciones de 
Durán (1981 y 1982); Poblete y Peredo 
(1990); Maloney (1993) y Moreno et al. 
(2005), quienes indican que la facilidad de 
absorción de agua está en función de la 
relación de compresibilidad del tablero, 
eficiencia del adhesivo y la facilidad de 
penetración del fluido (permeabilidad del 
tablero) por los espacios existentes entre 
partículas. Relaciones de compresibilidad 
bajas permiten al tablero absorber más 
humedad en comparación con tableros 
con alta tasa de compresibilidad. 

Los tableros encolados con adhesivo 
uf no son aptos para uso a la intemperie 
(Peredo y Torres, 1991) una exposición 
prolongada a la humedad lleva a la rup- 
tura de los enlaces que forman la adhe- 
sión entre los grupos activos de la madera 
(partículas) y el adhesivo (Stark et al., 
2010) afectando sus propiedades físico- 
mecánicas. Un calor excesivo también 
produce una degradación química del fra- 
guado de la resina uf, por lo que, típica- 
mente los tableros son enfriados 
rápidamente luego de finalizado el ciclo de 
prensado (Stark et al., 2010). 

La variación de espesor para los 
tableros de baja densidad se debe a la 
facilidad de penetración del agua en los 
capilares y espacios entre partículas, ori- 
ginando un mayor hinchamiento de la 
masa de partículas a este nivel de densi- 
dad. En los tableros de mayor densidad 
existe un mayor grado de compresibilidad, 
contribuyendo a disminuir el grado de 
higroscopicidad de las partículas de 
madera. 

Los porcentajes de variación en 
espesor después de 24 horas de inmer- 
sión obtenidos para los tableros de P. 
patula son similares a los encontrados por 
Moreno etal. (2005) en tableros con partí- 
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culas de F. americana a dos niveles de 
densidad (600 kg-nrr 3 - 800 kg nr 3 ) y por 
Guler et al. (2008) en tableros de partícu- 
las de Pinus nigra con una densidad de 
700 kg rrr 3 . 

Por otra parte, los resultados aquí 
obtenidos para la variación de espesor a 
24 horas de inmersión están en concor- 
dancia con lo que exponen Durán (1981) 
y Moreno et al. (2005 y 2010) que, a 
mayor densidad del tablero, mayor será 
el hinchamiento; pero es contarlo a lo 
mencionado por Poblete y Peredo (1990) 
quienes no encontraron una correlación 
adecuada entre hinchamiento y densidad 
del tablero. 

Aumentar el tiempo de inmersión a 
24 horas permite saturar paredes celula- 
res que por su bajo contenido de hume- 
dad tienen capacidad de absorción 
produciendo un hinchamiento en las partí- 
culas, manifestándose finalmente en la 
variación del espesor del tablero. Como 
se ha mencionado anteriormente existe 
variación de la densidad dentro del 
tablero, por lo que el hinchamiento del 
tablero tampoco será uniforme, es decir, el 
tablero se deforma y pierde estabilidad 
dimensional. 

El hinchamiento del espesor está 
directamente afectado por muchas varia- 
bles de proceso, pero las interacciones 
entre las variables pueden cambiar este 
comportamiento, por ejemplo, la durabili- 
dad del tablero depende de hinchamien- 
tos bajos siendo una práctica aplicar un 
alto contenido de resina adhesiva para 
controlar el hinchamiento del tablero. Sin 
embargo, Halligan (1970), menciona que 
varias investigaciones sobre la eficiencia 
de los tratamientos para el control del hin- 
chamiento han dado resultados contradic- 
torios y propone que un vaporizado 
post-tratamiento parece ser una opción 
prometedora. Al respecto, Peredo y Torres 
(1991), indican en su investigación que la 


aplicación de agentes hidrófobos reduce 
el efecto negativo de la humedad sobre la 
resina uf, disminuyendo significativa- 
mente la absorción de humedad e hincha- 
miento del tablero de partículas y en 
consecuencia mejorando las propiedades 
mecánicas del mismo. 

Gatchell et al. (1966), así como 
Roffael y Rauch (1972), exponen que 
incrementando la densidad del tablero 
podrían esperarse mayores hinchamien- 
tos; sin embargo, de acuerdo con Leh- 
mann (1960), el beneficio de obtener 
aumento de la resistencia mecánica es 
razón suficiente para compensar la ten- 
dencia de hinchamiento. 

El módulo de ruptura es aceptado 
como una medida de esfuerzo y repre- 
senta la carga máxima de un miembro que 
sometido a flexión llega a su fractura (Cai 
y Ross, 2010). Sin embargo, la resistencia 
al esfuerzo en flexión debe ser considera 
hasta el límite elástico, para efectos de 
seguridad (McNatt, 1973). Mientras que, 
la adhesión interna o tensión perpendicu- 
lar determina la resistencia del tablero a la 
delaminación o separación de sus capas 
(Cai y Ross, 2010). 

El mor en flexión estática se incre- 
menta a mayor densidad del tablero, esto 
se debe entre otros factores a una mayor 
eficiencia en la distribución de la resina y 
a un mejor contacto entre partículas-adhe- 
sivo, así como también a una mayor den- 
sificación de las partículas (mayor razón 
de compresibilidad), disminuyendo el 
volumen de espacios vacíos. Esto con- 
lleva a que el tablero de mayor densidad 
trasmita los esfuerzos de partícula a partí- 
cula de manera más uniforme traducién- 
dolos en un aumento de la resistencia 
mecánica (Durán, 1981; Maloney, 1993 y 
Bazyar et al., 2011). 

En general, el mor en flexión obte- 
nido en este estudio con los tableros de 
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ambas densidades (Fig. 5a) supera al 
intervalo de mor de 15,20 MPa - 21,97 
MPa (155 kg cnr 2 - 224 kg cnr 2 ) encon- 
trado por Cai y Ross (2010) para tableros 
de partículas tipo industrial con densidad 
de 690 kg rrr 3 . 

En cuanto a la tensión perpendicular 
a la superficie del tablero (tpst), aumenta a 
medida que se incrementa la densidad del 
tablero y la relación de compresibilidad. 
Esto se puede atribuir al mayor contacto 
relativo entre las partículas y la mayor efi- 
ciencia de la resina, así como también a 
que se logra una densificación mayor en 
el centro del tablero debido al mayor volu- 
men de partículas a ser prensadas. Por lo 
tanto, el gradiente de densidad en el espe- 
sor del tablero es menos variable que en 
el caso de los tableros de menor densidad 
(Durán, 1982; Maloney, 1993 y Moreno et 
al., 2002). 

En general, la adhesión interna obte- 
nida en este estudio con los tableros de 
ambas densidades está dentro de los 
intervalos mencionados por Cai y Ross 
(2010), de 0,314 MPa - 1,177 MPa (3,2 kg 
cm 2 - 12,0 kg cnr 2 ) para tableros de partí- 
culas de tipo industrial con densidad de 
690 kg nr 3 , y el intervalo de 0,314 MPa - 
0,490 MPa (3,2 kg cnr 2 - 5,0 kg cnr 2 ) 
encontrado por Guler et al. (2008) para 
tableros de partículas con una densidad 
de 700 kg rrr 3 fabricados con residuos de 
Pin us nlgra. 

Finalmente, el modelamiento con el 
software estadístico Design Expert con 
una confiabilidad de 95% indica que, a 
mayor densidad, disminuye la absorción 
de humedad a ambos tiempos de inmer- 
sión, aumenta el hinchamiento o varia- 
ción del espesor a 24 horas de inmersión. 
Los esfuerzos en flexión estática y ten- 
sión perpendicular (adhesión interna) 
también aumentan. Estos resultados son 
coincidentes con lo presentado en las 
figuras 4 y 5. 


CONCLUSIONES 

Como se demostró en el análisis de 
varianza y en los resultados obtenidos, la 
densidad del tablero es un factor que 
influye de manera directa sobre las pro- 
piedades físico-mecánicas de un tablero 
aglomerado de partículas. 

El aprovechamiento de los residuos 
industriales de la especie Pinus patula, 
mostró buenos resultados y propiedades 
tecnológicas adecuadas para su trasfor- 
mación en un producto de ingeniería 
como el tablero homogéneo de una capa 
de partículas con densidad media. Los 
tableros de este estudio cumplieron con 
los valores mínimos solicitados por las 
normas, en particular, para las propieda- 
des mecánicas de mor y adhesión 
interna. Sin embargo, los porcentajes de 
absorción de humedad y variación de 
espesor en algunos casos no cumplieron 
con los valores mínimos aceptados por 
norma. 

De acuerdo con los resultados obte- 
nidos, el tablero homogéneo de partículas 
de P. patula con una densidad de 767 kg 
rrr 3 presentó las mejores características 
en cuanto a propiedades mecánicas de 
mor, adhesión interna y mejor comporta- 
miento en la absorción de humedad a 2 y 
24 horas y variación de espesor a 2 horas. 

Finalmente, los tableros aglomera- 
dos de partículas de P patula homogé- 
neos de una capa, encolados con 
urea-formaldehído y de ambas densida- 
des poseen buena resistencia mecánica y 
adhesión interna, por lo que se reco- 
mienda su uso en partes internas de mue- 
bles y gabinetes, como material decorativo 
superpuesto sobre otros materiales, 
como recubrimiento interno de casas pre- 
fabricadas, entre otros usos; siempre y 
cuando el tablero preste su servicio en 
sitios con poca o nula exposición directa a 
la humedad. 
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Opciones reales para determinar el turno 
óptimo en sistemas silvopastoriles: 
análisis de inversión 

Real options for determining optimum t me 
Systems sylvopastoral: investment analysis 

Gastón Milanesi , 1 Diego Ricardo Broz 2 
y Guillermo Bruno Woitschach 3 

RESUMEN 

El trabajo propone analizar y valorar, desde la perspectiva económica-financiera, el momento 
óptimo de cosecha en un sistema silvopastoril, utilizando las Opciones Reales (OR) como herramienta 
para la toma de decisiones estratégicas. Ésta se basa en un proceso estocástico para la proyección de 
precios, mientras que el Valor Actual Neto (VAN) utiliza una tasa fija en todo el proyecto. Los resultados 
mediante OR indican talar a los 18 años con un VOR (Valor Opción Real) $AR 1 088 339 frente a la 
alternativa VAN de talar a los 11 años y $AR 511 914. 

PALABRAS CLAVE 

Flujo de fondos, momento óptimo, opciones reales, rejilla binomial. 


ABSTRACT 

ln this work we assess, from an economic and financial perspective, the optimal time of timber 
harvest in joint forestry-pasture Systems. For this we apply a Real Options (RO) approach as a tool of 
strategic decisión making. It assumes a stochastic process of price projections and is compared with the 
present valué (NPV) computed on the basis of a fixed rate all along the project. The results obtained by 
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yields $AR 511 914. 

KEY WORDS 

Cash flow, optimal timing, real options, binomial lattice. 


1 Departamento de Ciencias de la Administración, Universidad Nacional del Su, CONICET, Bahía Blanca, 
Argentina. C. e.: milanesi@uns.edu. ar 

2 Doctorando en Ingeniería, Departamento de Ingeniería, Universidad Nacional del Sur, CONICET, Bahía 
Blanca, Argentina. C. e.: diego. broz@uns.edu. ar 

3 Doctorando en Cs. de la Administración, Departamento de Cs. de la Administración, Universidad Nacional del 
Sur, CONICET, Bahía Blanca, Argentina. C. e.: gwoitschach@facfor.unam.edu. ar 



82 


Opciones reales para determinar el turno óptimo en sistemas silvopastoriles 


INTRODUCCIÓN 

Los sistemas agroforestales o mixtos, 
específicamente los sistemas silvopastori- 
les (ssp) están teniendo auge en Argentina 
debido a la flexibilidad en el flujo econó- 
mico, permitiendo la generación de ingre- 
sos anuales a partir del cuarto año del 
establecimiento; producto del componente 
ganadero y del ingreso final del compo- 
nente forestal (Fassola et al., 2004). Son 
varios los autores que afirman bondades 
del ssp. Musálem (2002) menciona que los 
sistemas pueden estar en zonas fértiles, 
como también en zonas áridas, mejorando 
la productividad de estas últimas; Blasco et 
al. (2005) lo plantean como una estrategia 
económica para combatir el vinal ( Prosopis 
ruscifolia ) frente a métodos tradicionales; 
Cerruto et al. (2007) plantea que los ssp 
permiten abatirla erosión en comparación a 
una pastura degradada, aportando aumento 
de la viabilidad financiera; para Soares et 
al. (2003) los ssp son una alternativa de 
d ¡versificación de rentas y recuperación 
ambiental. Se concluye que los ssp tienen 
importancia desde el punto de vista econó- 
mico y ecológico, por ello se requieren 
herramientas para su gestión como una 
manera de tomar decisiones óptimas. 

En este tipo de inversiones se asocia 
a la maduración del proyecto con la rota- 
ción, éste es un término relacionado con 
el tiempo que lleva una población en ser 
cortada (Schneider, 2002), convirtiendo el 
capital en renta. El mismo autor menciona 
que la duración de la rotación es definida 
por la meta técnica de la empresa, siendo 
los criterios tanto económicos como bioló- 
gicos. Bettinger et al. (2009) enumeran 
siete tipos de rotación: física, técnica, silvi- 
cultura!, biológica, por generación de 
ingresos, económica-financiera y valora- 
ción del crecimiento porcentual. 

En general, el método empleado 
para valorar activos reales es el descuento 
de flujo de fondos (dff) (Copeland et al., 


2000; Pratt y Grabowski, 2008). No obs- 
tante, este método no incorpora en el 
valor los múltiples cursos de acción o fle- 
xibilidad estratégica, propios de la vida 
proyectada e incierta de la inversión en 
cuestión (Mun, 2005). Valorar la flexibili- 
dad estratégica de una inversión se logra 
adaptando la Teoría de Opciones para ins- 
trumentos financieros (Black y Scholes, 
1973; Merton, 1973) a problemas de 
inversión en activos reales, dando lugar a 
los modelos de Valoración de Opciones 
Reales (vor) (Dixit y Pindyck, 1994; Tri- 
georgis, 1997). En la literatura financiera 
el primer modelo corresponde a valorar la 
opción de crecimiento (Myers, 1977), a 
partir de este trabajo se han desarrollado 
diferentes propuestas analíticas para el 
tratamiento de categorías específicas de 
opciones. A continuación se enumeran 
algunos trabajos seminales para cada 
tipo: a) opción de diferimiento (Me Donald 
y Siegel, 1981; Paddock et al., 1988); b) 
opción de crecimiento (Pindyck, 1988; Tri- 
georgis, 1988; Smit, 1997); c) opción de 
abandono (Myers y Majd, 1990); d) opcio- 
nes de expandir-contraer o extensión de 
la vida útil (Trigeorgis y Masón, 1987); e) 
opción de cierre temporario o corte del 
proceso productivo (Brennan y Schwartz, 
1985); f) opción de intercambio (Mar- 
grabe, 1978; Kulatilaka, 1988); g) opcio- 
nes financieras de insolvencia (Masón y 
Merton, 1985). Si bien los clásicos mode- 
los plantean ecuaciones en diferencia 
estocástica en tiempo continuo, para el 
análisis estratégico de activos reales se 
emplea el método binomial (Cox et al., 
1979), utilizando rejillas en lugar de árbo- 
les, dada su capacidad de recombinar 
nodos intermedios (Smith, 2005; Milanesi 
y Vigier, 2010). 

En el campo forestal los pioneros en 
la aplicación de componentes estocásticos 
en cosecha fueron Malliaris y Brock 
(1982), Miller y Voltaire (1983), entre 
otros. Thomson (1992), así como Yin y 
Newman (1997) estudiaron el turno óptimo 


Madera y Bosques 19(3), 2013:81-98 


83 


de cosecha, asumiendo que los precios 
siguen un movimiento geométrico brow- 
niano, estos últimos autores mencionan 
que los enfoques tradicionales, como la 
Metodología de Faustmann (Grainger, 
1 968), se limitan a un mundo determinista. 

Para Sant'Anna y Nogueira (2010), 
la vor es una herramienta más realista de 
evaluación que el valor actual neto (van), 
ya que la tasa de descuento permanece 
inalterable en la vida del proyecto al reali- 
zar el análisis. El autor recalca que, en 
valoraciones de concesiones forestales, la 
vor presenta valores significativamente 
mayores que van, esto coincide con 
Moreira et al. (2000), quien menciona que 
los resultados obtenidos por vor fueron 
140% superiores a los calculados por van 
y más realistas. 

Soares et al. (2003) mencionan que 
la tasa de descuento afecta significativa- 
mente el van, siendo más rentable el pro- 
yecto cuanto menor sea la tasa. La tasa 
de descuento no se elabora mediante un 
procedimiento cuantitativo riguroso, sino 
mediante un análisis cualitativo integral, 
no riguroso o de poca consistencia, con- 
dicionado por costo de oportunidad, 
grado de riesgo del proyecto, costo del 
dinero, entre otros; afectando directa- 
mente en la toma de decisión sobre pro- 
yecto. Entre los trabajos que proponen 
alguna tasa encontramos a Monicault y 
Delvalle (2009) con una tasa de des- 
cuento de 12% obteniendo un van de 
$AR 9 357 y un turno de 12 años en el 
estudio de un ssp en la Provincia del 
Chaco, Argentina; Chidiak et al. (2003) 
plantea un estudio de factibilidad de 
inversiones forestales en la Patagonia 
con tasas de 8% a 12%, además men- 
ciona que por debajo de 8% el van es 
negativo; Colcombet et al. (2010) con- 
cluye que un proyecto forestal a 152 km 
de la industria resiste holgadamente 
tasas de descuento de 8% anual en la 
Provincia de Misiones; Fassola et al. 


(2004) mencionan Ingresos Anuales 
Equivalentes (iae) mayores a $AR 500 
con una tasa de 5%, siendo superior el 
beneficio de los ssp frente a forestación 
pura, invernada de novillos, maíz y soja. 

En este trabajo se plantea la necesi- 
dad de aplicar el enfoque de la vor en el 
ejercicio de la tala en un sistema de pro- 
ducción mixto forestal-ganadero, compa- 
rando cada año, a partir del décimo 
(madurez tecnológica) hasta el vigésimo 
(máximo horizonte de planificación), el 
valor obtenido, producto de la cosecha en 
los nodos de la rejilla binomial (año t) ver- 
sus el valor actual neto esperado de talar 
en el próximo año (año t+1). La compara- 
ción periodo a periodo es relevante por- 
que estos valores se caracterizan por ser 
cambiantes en función del tiempo y riesgo 
de las variables. El caso de análisis se 
corresponde con un escenario forestal 
simulado y un nivel de producción gana- 
dera acorde con el patrimonio forestal, el 
cual permitirá determinar la viabilidad eco- 
nómica del proyecto de inversión para el 
noreste de la República Argentina, 
mediante la aplicación de la vor y la com- 
paración con el van. 


OBJETIVOS 

El primero objetivo de este trabajo fue la 
determinación de la madurez del pro- 
yecto de inversión para un sistema de 
producción mixto forestal-ganadero para 
el noreste argentino, el cual se corres- 
ponde con la determinación del turno 
óptimo de cosecha mediante la aplica- 
ción de or. Como segundo objetivo se 
planteó la comparación de resultados 
entre un escenario evaluado mediante la 
vor y por otro lado el van. Como tercer y 
último objetivo se planteó fundamentar el 
uso de la vor como metodología de eva- 
luación económica-financiera de un sis- 
tema de producción mixto frente al uso 
del van. 
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MATERIALES Y MÉTODOS 


Valoración del componente forestal en 
el SSP 

Para evaluar el proyecto de inversión, se 
proyectaron a partir de los años 1 0 y hasta 
un horizonte de planificación de 20 años, 
los ingresos y egresos de cada una de las 
componentes del sistema para luego apli- 
car la vor y el van y, consecuentemente, 
tomar las decisiones estratégicas que 
generaron mayor beneficio. La evolución 
de la masa forestal fue obtenida mediante 
el software Simulador Forestal (Crechi et 
al., 2005) y la producción ganadera (carga 
por hectárea) se supuso constante desde 
el inicio hasta el final del horizonte de pla- 
nificación. 

La volatilidad se estimó a partir de la 
serie histórica de precios correspondientes 
a los subproductos de la explotación fores- 
tal, considerados estos como cartera de 
productos y aplicando los conceptos de la 
Teoría Moderna de la Cartera (tmc) 
(Markowitz, 1 991 ). 4 Por otro lado, para el 
componente ganadero se utilizó solamente 
la serie histórica, ya que a diferencia del 
forestal, existe un sólo producto comercial. 
Para calcular el precio promedio y la varia- 
bilidad de la mezcla de subproductos obte- 
nidos, aplicando los conceptos de la tmc, 
estimando el precio medio se utilizó la 
ecuación (1), el riesgo del precio es deter- 
minado por la ecuación (2). 

/■ -V.v/' í 1 » 


4 La Teoría de la Cartera (Modern Portfolio 
Theory) es aplicable a todo tipo de carteras de 
activos. En este caso, el precio promedio de la 
masa forestal es el promedio ponderado del 
promedio observado de cada subproducto. La 
volatilidad es representada por la suma ponde- 
rada de los desvíos observados ajustados por 
la covarianza (correlación) precio correspon- 
diente a la familia de subproductos. 


donde: 

P i : precio promedio de cada subpro- 
ducto 

X t : participación en el total producido 
P o : precio promedio 


= Yxfof + x^oj + 2(.v,.v, !/',<’,<>. (2) 

í = t7=i 

donde: 

o : riesgo del precio de la cartera de 
productos 

x yx : participaciones de cada subpro- 
ducto 

a- y cr : varianzas de los precios de los 
subproductos 

cr. y cr. : desvíos de los precios 
p 1 : coeficiente de correlación de precios 


Usando el precio promedio de la ecuación 
(2), se calcularon los flujos de fondos 
libres con la ecuación (3). 

FF¡ = Ló xvo/,]x(l -C,)x(l-/) O) 
donde: 

FF '■ flujos de fondos libres 
vol t producción esperada en cada pe- 
riodo del horizonte de planificación 
obtenido mediante simulación 
C t : porcentaje que representan los cos- 
tos de producción (variables y de 
estructura) sobre los ingresos por 
ventas, calculado mediante la ecua- 
ción (4). 

/ : tasa marginal de impuesto a las 
ganancias. 

La tasa marginal de impuesto a las ganan- 
cias / fue de 35%. 

l p <> x vo/J 
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VAN ft = FF t x e~ k (5) 

El valor intrínseco de la explotación sin fle- 
xibilidad estratégica se estimó con la 
ecuación (5), siendo su expresión el valor 
actual de los flujos de fondos derivados de 
la tala en cada uno de los años, desconta- 
dos a la tasa ajustada por riesgo sobre 
inversiones equivalentes. 5 La estructura 
de la tasa para mercado emergente, poco 
diversificado y segmentado (Damodaran, 
2002), surgió de ajustar el rendimiento de 
mercado por la volatilidad total, en este 
caso dada por el desvío del precio. La 
ecuación (6) muestra la formulación de la 
tasa ajustada por riesgo. 

k = (l + rm)x (l + cr.VA/') (®) 

donde: 

rm : adicional por riesgo de mercado 
o '■ variación del precio 
A t ’■ relación entre Tin , siendo esto la 
longitud de cada paso del modelo 
k : relación porcentual 


La proyección del precio de venta prome- 
dio de los subproductos forestales sigue 
un proceso estocástico del tipo Geomé- 
trico Browniano. Se trabajó en tiempo dis- 
creto y, para la elaboración de una rejilla 
binomial, 6 se utilizaron los parámetros del 


5 La tasa por riesgo equivalente se calcula apli- 
cando los clásicos modelos como el CAPM 
{Capital Asset Prlcing Model), APM {Arbitrage 
Pricing Model): MPM {Multifactor Pricing Model) 
y sus adecuaciones a mercados emergentes. 

6 Las rejillas binomiales son empleadas debido a 
su capacidad de recombinar nodos e incorporar 
las propiedades propias del triángulo de Pascal. 
Los árboles binomiales explicitan gráficamente 
el problema de manera más eficiente. No obs- 

tante, no permiten recombinar, lo que torna su 
cálculo computacionalmente más complejo que 
las rejillas. Si n representa los periodos o inter- 

valos de tiempo, la cantidades de nodos finales 
de una rejilla es de n+l y de nodos intermedios 

n(n+l)/2. En los árboles binomiales, en cambio, 
estos valores son 2n-l y 2n-l, respectivamente 
(Smith, 2005). 


modelo crr (Cox et al., 1979) mediante 
las ecuaciones (7) y (8). 


d = (8) 

Los movimientos de ascenso u y des- 
censo d se calcularon utilizando el des- 
vío de la cartera de precios calculado con 
la ecuación (2). El precio de cada nodo 7 * 
padre p de la rejilla en t - 1 genera los 
precios ascendentes i y descendentes j 
en t por medio de la ecuación (9). 

P u{t) = Lfy-1) x u * p j(t- 1) x d \ ( 9 ) 


Una vez proyectado el proceso estocás- 
tico del precio por medio de las ecuacio- 
nes (7) y (8), se calculó nuevamente el 
flujo de fondos para cada periodo t y cada 
nodo (i,j) mediante la ecuación (10). 

FF t (,) - !/>,(,) X vol, Jx (l - c, ) (1 0) 


Con los flujos de fondos para cada nodo 
de la rejilla se construyeron las or, para 
ello fue necesario definir el precio de ejer- 
cicio de la opción. La decisión a que se 
enfrenta el agente es seleccionar el 
momento de madurez de la inversión, que 
se corresponde con la determinación del 
turno óptimo de cosecha del ssp, para lo 
cual se debe seguir un criterio de maximi- 
zación de beneficios, conjugando el factor 
tiempo y el riesgo que emerge de las 
variables precio de los subproductos. Los 
flujos de fondos proyectados en cada 


7 Los nodos pueden interpretarse como posibles 
escenarios futuros; a mayor intervalo de tiempo 
mayor cantidad de nodos debido a la mayor 
incertidumbre. Así, para el periodo t = 1 la canti- 
dad de nodos de la rejilla es de 2 (dos escena- 
rios) mientras que para t= 10 la cantidad es 11. 
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escenario de la rejilla fueron comparados 
con el valor esperado del próximo año, 
según la expresión (11). 


n 



( 11 ) 


Viger, 2010; Milanesi, 2011) dado en la 
ecuación (15). 


V r 


0 1 



i=\J=\,t 


X 



(15) 


Donde x es el promedio de los flujos de 
fondos en cada nodo del próximo periodo 
FF v(t + \) ponderado por la probabilidad aso- 
ciada <lij(t+ 1 ) , siendo las probabilidades 
asociadas a cada nodo dadas por la ecua- 
ción (12). 


( 12 ) 

En este caso, los éxitos m estuvieron 
dados por los movimientos de ascenso 
(subas) del precio sobre el total de obser- 
vaciones n y condicionados por los coefi- 
cientes equivalente-ciertos P , estimados 
con la ecuación (13). 

(1 + rf-d 03) 

p = 

u — d 


Siendo v ot el valor actual de la opción de 
talar en el momento t . El ejercicio de la 
opción en cada nodo (escenario) en t se 
expresa por V m . La probabilidad de 
alcanzar el nodo asociado es ^(0 ■ El fac- 
tor de actualización desde el momento (t) 
hasta el inicio (f = 0 ) es e~ rAt - 

Valoración del componente ganadero 
en el ssp 

El análisis del componente ganadero fue 
análogo al del componente forestal. La 
particularidad de este componente es la 
generación de flujos de cajas anuales, 
donde hubo egresos por compra de terne- 
ros y costos de producción y un ingreso a 
los doce meses producto de la venta. El 
van se determinó con la ecuación (16) 


VAN 0 = (lG t -EG r )x e~ kGxt ( 1 6) 


Donde T representa la tasa libre de 
riesgo y los movimientos (u\ d), surgen de 
las ecuaciones (7) y (8) y A t se definió en 
la ecuación (6). El valor esperado v m sur- 
gió de comparar en cada periodo, el valor 
correspondiente a cada nodo FF m con el va- 
lor esperado en el próximo año;x ?+1 xe~ rAt 
matemáticamente el problema se resume 
en la expresión (14). 


donde 

IG '■ ingreso por venta del componente 
ganadero al final del período t 
EG t '■ egresos debido a costos operativos 
y compra del ganado al inicio de t 
kG '■ tasa ajustada por riesgo obtenida 
para el componente ganadero me- 
diante la ecuación (17). 




’ SÍFF v( t ) 

si no 


>X, +1 xe 


-rAt 


(14) 


El valor de la opción en el momento inicial 
{t = 0) requiere de utilizar el procedi- 
miento recursivo sobre la rejilla (Milanesi y 


kG = (l + rm)x (l + cr g VAr) (17) 


El precio inicial surgió de aplicar la ecua- 
ción (18). 


j = 



x 


n 


(18) 
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donde: 

x¡ : observación de precio mensual 
* : media de la serie de precios 
n : observaciones. 


El riesgo del precio cr u se obtuvo a partir 
de la serie de precios mensuales apli- 
cando la ecuación (19) 



donde: 

N : número de observaciones 

(jc — jc . ) / jc . : variación relativa del precio 


Análisis de sistema mediante la vor y el 

VAN 

En esta etapa convergen los resultados 
con el objeto de un análisis integral del 
sistema. Para la obtención de un VAN total 
(van t ) se sumó en t los van f y van g - de 
la misma manera se procede con el s VOR 
total ( vor ) ya que los valores son aditi- 
vos (Rocabert, 2007). En las ecuaciones 
(22) y (23) se puede ver la expresión 
matemática. 


VAN Tt = VAN ft + VAN gt (22) 

VOR Jt = VOR ft + VOR gt (23) 


Para el cálculo de las tasas de ascenso y 
descenso ( u , d) se aplicaron las ecuacio- 
nes (7) y (8). Se realizaron las proyeccio- 
nes de los precios estocásticos mediante 
la ecuación (9). Se calculó el nuevo flujo 
de fondo mediante la ecuación (20). 


El turno óptimo (to) de cosecha forestal 
en ssp se da cuando vor T{ n > vor tu+1) t S u 
expresión matemática es la (24); mientras 
que en el caso del VAN, la cosecha se da 
cuando se obtiene el máximo VAN Tt . 


FF, á ,^(p, m xQ,)-EG, (20) 


En este caso: 

P ¡j{t) : precio promedio 
O : cantidad total (en ~¡N,) de venta de 
ganado 
EG t : egreso 

Posteriormente, se calculó la probabilidad 
para cada nodo mediante la ecuación 
(12). Se obtuvo el flujo de fondo en fun- 
ción al equivalente cierto (21) que permiti- 
ría obtener V (t) mediante la ecuación (14). 


< 21 ) 

i 

El valor de la opción v 0t se obtuvo apli- 
cando la ecuación (15). 


VOR, 


VOR 


no’ 


si \ OR TU) > J ’OR T y +Y) 
si no — » No Talar 


Talar 


(24) 


Si el valor de los flujos de fondos libres 
FF ij{i} supera al valor actual esperado del 
próximo año X t+l xe~ rAt ; se ejerce la tala, 
siendo su valor de V lJ = FF ij{¡) . Caso contra- 
rio, si el va-lor esperado para el próximo 
año supera los flujos de fondos libres del 
nodo, el resultado es de = 0 difiriendo la 
tala al año próximo (Plantinga, 1998). La 
decisión en cuestión se asimila a una 
opción exótica barrera o barrier option 
(Hull, 2005). 


Valoración de un sistema silvopastoril: 
caso de aplicación 

Se propuso analizar el turno óptimo de 
cosecha de un proyecto de inversión mixto 
mediante la aplicación de opciones reales 
a un predio de 39,8 ha sobre el cual se 
reforestaría con Pinus taeda L. La producción 
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Tabla 1. Egresos e ingresos del componente ganadero para un ciclo de 12 meses. 


Compra 

94.656 

MO 

91 000 

Alimento y sanidad 

20 776 

Costo Total 

206 432 

Ingreso por venta 

336 240 


Fuente: Elaboración propia. Sistema Integrado de Información Agropecuaria y Subsecretaría de Agricultura. Moneda: $AR 


se estimó mediante la simulación del cre- 
cimiento de la masa boscosa. La franja de 
edad evaluada fue a partir del año 
10 (f = 0), tomando un horizonte de plani- 
ficación de 10 años. Los objetivos del cul- 
tivo y manejo del rodal estaban alineados 
con aquellos tendientes a obtener subpro- 
ductos de alto valor; para ello, los destinos 
eran la industria de laminado y la industria 
de aserrío de grandes dimensiones 
(0 > 30 cm); el resto de los subproductos 
eran para aserrío de menores dimensio- 
nes (0 < 30 cm) y pulpa celulósica. 

En la simulación del manejo de la 
masa forestal se previo el inicio del ciclo 
con 1600 plantas por hectárea, 3 interven- 
ciones de aclareos (raleos) con una inten- 
sidad de 50% a 55% y periodicidad de 2 
años entre intervención, lo que permite 
obtener madera de alto valor favoreciendo 
el volumen individual (Crechi et al., 2005). 
Se planificó la ejecución de podas, que 
serían realizadas de acuerdo con la obten- 
ción de un cilindro defectuoso menor a 12 
cm y una altura total de poda de 6,5 m. 

La simulación se llevó a cabo 
mediante el software Simulador Forestal 
(Crechi et al., 2005). Se consideró un 
índice de sitio de 22 m (edad índice 15) y 
una densidad final de 190 árboles por hec- 
tárea. Los parámetros iniciales (año 10) 
del bosque fueron: altura media 16,4 m; 
diámetro medio 29,9 cm; área basal 13,5 
m 2 /ha, un volumen de 108,7 m 3 /ha y un 
idr (Prodan et al., 1 997) de 256. 


El componente ganadero estuvo 
constituido por la raza Braford. Se estimó 
alrededor de 2 animales por hectárea 
requiriendo un total de 80 individuos con 
un peso inicial de aproximadamente 120 
kg y un precio de compra de $AR 9,86/kg, 
el peso final para comercialización se esti- 
pula en 450 kg y un precio de venta de 
$AR 9,34/kg. El manejo correspondió al 
tipo invernada con un ciclo de 12 meses. 
Los costos e ingresos se detallan en la 
tabla 1. 


RESULTADOS 

Estimación de los flujos de fondos 

El precio promedio inicial surgió de aplicar 
la ecuación (1), arrojando un valor de $AR 
189 por tonelada y aplicando la ecuación 
(17) un valor de $AR 6,16/kg de peso vivo 
para el periodo t = 0. Para el riesgo de pre- 
cio forestal se aplicó la ecuación (2), utili- 
zando cuatro series históricas de precios 
correspondientes a cuatro clases 
diamétricas; 8 para el riesgo de precio 
ganadero se aplicó la ecuación (19), utili- 
zando una serie de precios históricos de 
ganado en pie 9 , los valores obtenidos para 


8 Fuente: Colegio de Ingenieros Forestales de 
Misiones (COIFORM). Datos mensuales para el 
Norte de la Provincia de Misiones entre los 
períodos julio de 2002 a diciembre de 2011 . 

9 Fuente: Sistema Integrado de Información Agro- 
pecuaria y Subsecretaría de Agricultura. Datos 
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Tabla 2. Valor Actual Neto del sistema entre el año 10 y 20, en gris se señala el turno. 


Año 

10 

11 

12 

13 

14 

15 

16 

17 

18 

13 

20 

VAN Forestal 

378 224 

382 106 

360 674 

324 062 

283 374 

241 659 

202 137 

167139 

136 462 

110395 

88 657 

VAN Ganadero 

129 808 

129808 

129 808 

129 808 

129 808 

129808 

129 808 

129808 

129 808 

129808 

129 808 

VAN Total 

508 032 

511 914 

490 482 

453 870 

413182 

371 468 

331 946 

296947 

266 270 

240203 

218.466 


Fuente: Elaboración propia. Moneda: $AR 

la cartera de subproductos forestales (<7 s ) 
fueron de 13,57% y 38,3% (cr g ) para el 
componente ganadero. 


Descuento de flujos de fondos 

Empleando los flujos de fondos se deter- 
minó el valor actual de las potenciales 
talas (ecuación (5)). La tasa ajustada por 
riesgo (ecuación (6)) fue de 28,39% 10 para 
el componente forestal y 34, 58%^ para el 
componente ganadero. Debido al creci- 
miento de la masa forestal, los mayores 
flujos de fondos se obtuvieron en el último 
periodo. Por la sensibilidad de los benefi- 
cios al factor de actualización, el criterio 
del descuento de flujos de fondos indicó 
que el mayor valor actual se produjo en el 
instante t = 1 (año 11), el resultado sugiere 
que la tala se realice de manera inminente 
con un van de $AR 511 914. En la tabla 2 
se pude ver la evolución del van en el hori- 
zonte analizado. Se aprecia que el turno 
de cosecha del componente forestal en 
forma independiente es igual al turno ana- 
lizado como sistema. 


Opciones reales 

Para la construcción del proceso estocás- 
tico de la variable de riesgo a los flujos de 


mensuales para la Provincia de Misiones entre 
los períodos julio de 2008 a diciembre de 2011 

1 0 Para ello se supuso un rendimiento de mercado 
de 25% (suponiendo nivel esperado de infla- 
ción). 

11 Para ello se supuso un rendimiento de mercado 
de 25% (suponiendo nivel esperado de infla- 
ción). 


fondos se incorporó el riesgo de precio. 
Conforme fue expresado, este ascendió a 
13,57% y 38,30% para el componente 
forestal y ganadero, respectivamente. El 
movimiento de ascenso y descenso, cal- 
culado mediante las ecuaciones (7) y (8), 
fue de u = 1,1453 y d = 0,8731 para el 
componente forestal y u = 1,4667 y 
d = 0,6818 para el componente ganadero. 
Para la construcción de la rejilla se aplicó 
la ecuación (10) y se estimaron los flujos 
de fondos mediante la ecuación (11). 
Mediante la ecuación (13), se determina- 
ron las probabilidades binomiales para 
cada nodo siendo el precio de ejercicio 
determinado con la ecuación (12). 

Los coeficientes equivalentes ciertos 
(ecuación (13)) para el componente fores- 
tal fueron p = 0,6545 y 7-p = 0,3455 y para 
el componente ganadero p = 0,4707 y 
1-p = 0,5293. Conforme fue expuesto, se 
acepta la tala (o madurez de la inversión) 
en el nodo proyectado si los flujos de fon- 
dos (ecuación (1 0)) superan al valor actual 
esperado en el próximo ejercicio (ecua- 
ción (11)). Para estimar el valor actual, el 
precio de ejercicio se actualizó al tipo sin 
riesgo, Milanesi et al. (2012) sugiere una 
tasa r = 5% para la Argentina. Para este 
cálculo se aplica la ecuación (14). 

Para estimar el periodo conveniente 
para la tala desde el punto de vista econó- 
mico-financiero, la información de la rejilla 
se resumió en una medida de valor espe- 
rado (ecuación (15)). Los diferentes resul- 
tados obtenidos en cada nodo, deben 
ponderarse por sus probabilidades de 
ocurrencia (ecuación (12)), obteniendo de 
esta manera el valor del ejercicio (tala). 
Para llevarlos al mismo horizonte tempo- 
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Opciones reales para determinar el turno óptimo en sistemas silvopastoriles 



■ VAN ForeslaS 

* — + - — VAN Ganadero 
+ 4 *— Total 
«» — -«» — VOR Fonstaf 

J ' VOR Ganadera 

^ mn ** m * VOR Total 


Años 


Figura 1 . Comparación de la evolución del van y vor por subsistema y total en $ar 


ral, se actualizó por el factor libre de riesgo 
( r ). De la tabla 3 a la tabla 1 2, se pueden 
ver los principales pasos para el desarro- 
llo de la metodología para cada compo- 
nente del sistema, mientras que en la 
tabla 13 se integran los resultados de las 
componentes del sistema para obtener un 
análisis integral. La decisión en juego en 
el último periodo se balanceó entre talar o 
abandonar, puesto que no había camino 
estratégico alternativo (venta o diferi- 
miento). 


A diferencia del criterio tradicional, el 
enfoque de opciones reales incorpora 
variables estocásticas (riesgo de precio) y 
la flexibilidad estratégica del ejercicio de 
talar. Entre los periodos t = 0 hasta t = 9, el 
máximo valor de ejercicio de la tala surgió 
en el periodo t = 8 (año 18) con una vor 
de $AR 1 088 339. 

Estos resultados difieren del criterio 
tradicional de descuento de flujo de fon- 
dos en más de 100%. 


Tabla 3. Proyección de precios promedios del componente forestal. 


0 1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

189 216,6 

248,1 

284,2 

325,5 

372,8 

426,9 

489,0 

560,0 

641,4 

734,6 

165,2 

189,2 

216,6 

248,1 

284,2 

325,5 

372,8 

426,9 

489,0 

560,0 


144,2 

165,2 

189,2 

216,6 

248,1 

284,2 

325,5 

372,8 

426,9 



125,9 

144,2 

165,2 

189,2 

216,6 

248,1 

284,2 

325,5 




109,9 

125,9 

144,2 

165,2 

189,2 

216,6 

248,1 





96,0 

109,9 

125,9 

144,2 

165,2 

189,2 






83,8 

96,0 

109,9 

125,9 

144,2 







73,2 

83,8 

96,0 

109,9 








64 

73,2 

83,8 









55,8 

63,9 

Precio Promedio: ecuaciones (1) y (9). 








48,7 


Fuente: Elaboración propia. Moneda: $AR 
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Tabla 4. Flujo de fondos libres del componente forestal. 


0 1 2 

3 

4 

5 

6 

7 

8 

9 

10 

581 883 837 539 1 153 997 

1 538143 

2 005 086 

2 566176 

3 234 572 

4 033 541 

4 979 071 

6 097 691 

7 416 481 

638 486 879 733 

1 172 581 

1 528 549 

1 956 287 

2 465 829 

3 074 912 

3 795 723 

4 648 486 

5 653 847 

670 652 

893 900 

1 165 267 

1 491 347 

1 879 789 

2 344 114 

2 893 614 

3 543 706 

4 310 129 


681 452 

888 324 

1 136 907 

1 433 030 

1 787 002 

2 205 905 

2 701 493 

3 285 764 



677 201 

866 705 

1 092 450 

1 362 295 

1 681 640 

2 059 444 

2 504 855 




660 720 

832 813 

1 038 526 

1 281 974 

1 569 987 

1 909 540 





634 883 

791 705 

977 294 

1 196 857 

1 455 710 






603 545 

745 026 

912 407 

1 109 739 







567 960 

695 560 

845 994 








530 250 

644 931 

Flujo de Fondo: ecuación (10) 








491 654 


Fuente: Elaboración propia. Moneda: $AR 

Tabla 5. Probabilidades asociadas a cada nodo 


0 1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

100,00% 65,45% 

42,84% 

28,04% 

18,35% 

12,01% 

7,86% 

5,14% 

3,37% 

2,20% 

1 ,44% 

34,55% 

45,23% 

44,40% 

38,75% 

31,70% 

24,90% 

19,01% 

14,22% 

10,47% 

7,61% 

100,00% 

11,94% 

23,44% 

30,68% 

33,47% 

32,86% 

30,11% 

26,27% 

22,11% 

18,09% 


100,00% 

4,12% 

10,80% 

17,67% 

23,13% 

26,49% 

27,74% 

27,23% 

25,46% 



100,00% 

1,43% 

4,66% 

9,16% 

13,98% 

18,30% 

21,56% 

23,52% 




100,00% 

0,49% 

1,93% 

4,43% 

7,73% 

11,38% 

14,90% 





100,00% 

0,17% 

0,78% 

2,04% 

4,01% 

6,56% 






100,00% 

0,06% 

0,31% 

0,91% 

1 ,98% 







100,00% 

0,02% 

0,12% 

0,39% 








100,00% 

0,01% 

0,05% 









100,00% 

0,00% 

Probabilidades: ecuación (13) 








100,00% 


Fuente: Elaboración propia 


Tabla 6. Comparación del valor correspondiente a cada nodo ( FF mn) con el valor 

esperado en el próximo año ( x t+l xe~ rAt ) 


0 1 2 

581 883 768 763 972 255 

3 

1 189 488 

4 

1 423 260 

5 

1 671 958 

6 

1 934 388 

7 

2 214 120 

8 

2 508 712 

9 

2 820 042 

10 

3 148 299 

0 1 2 

3 

4 

5 

6 

7 

8 

9 

10 

1 153 997 

1 538 143 

2 005 086 

2 566 176 

3 234 572 

4 033 541 

4 979 071 

6 097 691 

7 416 481 

- 


- 

1 956 287 

2 465 829 

3 074 912 

3 795 723 

4 648 486 

5 653 847 



- 

- 

- 

- 

2 893 614 

3 543 706 

4 310 129 








- 

3 285 764 








- 

2 504 855 








- 

1 909 540 








- 

1 455 710 






• 

• 

- 

1 109 739 







• 

- 

845 994 








. 

644 931 

Comparación de cada período: ecuación (14) 










Fuente: Elaboración propia. Moneda: $AR 
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Opciones reales para determinar el turno óptimo en sistemas silvopastoriles 


Tabla 7. Valor de la Opción Real para el componente forestal 


0 1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

- 

494 319 

431 222 

367 908 

308 173 

254 229 

207 490 

167 633 

134 363 

106 958 





620 115 

613 883 

584 525 

539 708 

486 664 

430 448 








760 212 

783 426 

779 546 










836 601 










589 203 










284 547 










95 429 










21 946 









- 

3 312 









■ 

296 


447 279 

371 156 

301 218 

722 951 

643 113 

558 124 

983 730 

895 519 

1 909 532 

Valor de la opción: ecuación (15) 











Fuente: Elaboración propia Moneda: $AR 


Tabla 8. Proyección de precios promedios del componente ganadero 


0 1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

6,16 9,04 

13,26 

19,45 

28,52 

41,83 

61,35 

89,98 

131,98 

193,57 

283,90 

4,20 

6,16 

9,04 

13,26 

19,45 

28,52 

41,83 

61,35 

89,98 

131,98 


2,87 

4,20 

6,16 

9,04 

13,26 

19,45 

28,52 

41,83 

61,35 



1,95 

2,87 

4,20 

6,16 

9,04 

13,26 

19,45 

28,52 




1,33 

1,95 

2,87 

4,20 

6,16 

9,04 

13,26 





0,91 

1,33 

1,95 

2,87 

4,20 

6,16 






0,62 

0,91 

1,33 

1,95 

2,87 







0,42 

0,62 

0,91 

1,33 








0,29 

0,42 

0,62 









0,20 

0,29 

Precio Promedio: ecuaciones (18) y (19) 








0,13 


Fuente: Elaboración propia. Moneda: $AR 


Tabla 9. Flujo de fondos libres del componente ganadero 


0 1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

15 448 118 995 

270 864 

493 608 

820 301 

1 299 454 

2 002 218 

3 032 946 

4 544 693 

6 761 939 

10 013 924 

-55 151 

15 448 

118 995 

270 864 

493 608 

820 301 

1 299 454 

2 002 218 

3 032 946 

4 544 693 


-103 287 

-55 151 

15 448 

118 995 

270 864 

493 608 

820 301 

1 299 454 

2 002 218 



-136 106 

-103 287 

-55 151 

15 448 

118 995 

270 864 

493 608 

820 301 




-158 483 

-136 106 

-103 287 

-55 151 

15 448 

118 995 

270 864 





-173 740 

-158 483 

-136 106 

-103 287 

-55 151 

15 448 






-184 142 

-173 740 

-158 483 

-136 106 

-103 287 







-191 234 

-184 142 

-173 740 

-158 483 








-196 070 

-191 234 

-184 142 









-199 367 

-196 070 

Flujo de Fondo: ecuación 

(20) 








-201 615 


Fuente: Elaboración propia. Moneda: $AR 
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Tabla 10. Probabilidades asociadas a cada nodo del componente ganadero 


0 1 

2 

3 

4 

5 

6 

7 

8 

9 

1000,00% 

100,00% 47,07% 

22,16% 

10,43% 

4,91% 

2,31% 

1 ,09% 

0,51% 

0,24% 

0,11% 

0,05% 

52,93% 

49,83% 

35,18% 

22,08% 

12,99% 

7,34% 

4,03% 

2,17% 

1,15% 

0,60% 

100,00% 

28,01% 

39,56% 

37,24% 

29,22% 

20,63% 

13,60% 

8,53% 

5,16% 

3,04% 


100,00% 

14,83% 

27,92% 

32,85% 

30,93% 

25,48% 

19,19% 

13,55% 

9,11% 



100,00% 

7,85% 

18,47% 

26,08% 

28,65% 

26,97% 

22,85% 

17,93% 




100,00% 

4,15% 

1 1 ,73% 

19,33% 

24,26% 

25,69% 

24,19% 





100,00% 

2,20% 

7,24% 

13,64% 

19,26% 

22,67% 






100,00% 

1,16% 

4,38% 

9,28% 

14,56% 







100,00% 

0,62% 

2,61% 

6,14% 








100,00% 

0,33% 

1,53% 









100,00% 

0,17% 

Probabilidades: ecuación (13) 









100,00% 


Fuente: Elaboración propia 


Tabla 11 . Comparación de flujo de fondo libre en t con el valor presente esperado t+1 


0 

1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

15 448 

26 824 

38 784 

51 356 

64 573 

78 468 

93 075 

108 431 

124 574 

141 545 

159 386 

- 

1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

- 

118 995 

270 864 

493 608 

820 301 

1 299 454 

2 002 218 

3 032 946 

4 544 693 

6 761 939 

10 013 924 



• 

118 995 

270 864 

493 608 

820 301 

1 299 454 

2 002 218 

3 032 946 

4 544 693 




- 

- 

118 995 

270 864 

493 608 

820 301 

1 299 454 

2 002 218 





- 

• 


118 995 

270 864 

493 608 

820 301 






■ 

- 

- 

■ 


270 864 

Comparación de 

cada período: 

ecuación (14) 








. 


Fuente: Elaboración propia Moneda: $AR 


Tabla 12. Valor de la Opción Real para el componente forestal 


0 1 2 

3 

4 

5 

6 

7 

8 

9 

10 

56 014 60 019 

51 486 

40 277 

30 034 

21 784 

15 533 

10 956 

7 674 

5 349 

- 

41 867 

59 814 

64 137 

60 208 

52 379 

43 418 

34 830 

27 297 

- 



34 769 

55 883 

67 114 

70 002 

67 114 

60 848 






30 319 

51 979 

66 883 

74 745 








- 

48 563 

53 282 54 308 

80 350 

81 947 

100 419 

102 140 

116 516 

118215 

112 542 

1 31 497 

Valor de la opción: ecuación (15) 










Fuente: Elaboración propia. Moneda: $AR 
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Tabla 13. Valor Opción Real del sistema entre el año 10 y 20, en gris se señala el turno 


Ario 

10 

11 

12 

13 

14 

15 

16 

17 

18 

19 

20 

VOR Forestal 



447 279 

371 156 

301 218 

722 951 

643 113 

558 124 

983 730 

895 519 

1 909 532 

VOR Ganadero 


44 464 

45 889 

68 392 

70 531 

86 907 

89 240 

84 290 

104 609 

100 760 

117 660 

VOR Total 

- 

44 464 

493 168 

439 548 

371 749 

809 858 

732 353 

642 414 

1 088 339 

996 279 

2 027 193 


Fuente: Elaboración propia. Moneda: $AR 


Además, el riesgo en el modelo van se 
incorporó íntegramente en el factor de 
actualización y los periodos fueron trata- 
dos sin flexibilidad alguna (solamente tala) 
y el resultado se encontró en parte influen- 
ciado por el impacto de la tasa ajustada 
por riesgo en el flujo de fondos esperado 
estático. En la figura 1 se expone el perfil 
de beneficios en ambos enfoques. 

DISCUSIÓN 

Se expuso un modelo alternativo al clá- 
sico método de descuento de flujos de 
fondos, para estimar la conveniencia eco- 
nómica-financiera y temporal de ejercer la 
opción de talar una masa forestal sujeta a 
un manejo silvopastoril, siendo esto 
correspondiente con la madurez del pro- 
yecto de inversión. Se propone aplicar la 
teoría de opciones reales para incorporar 
en el marco de decisión la posibilidad de 
diferir la tala y tomar en cuenta el proceso 
aleatorio del precio. La decisión consiste 
en comparar cada escenario proyectado 
versus el valor actual promedio de los 
escenarios del próximo ejercicio conforme 
a Plantinga (1998); en términos de opción 
financiera la decisión se asemeja a una 
barrier option del tipo knock-in. 

En el caso analizado, el método de 
descuento de flujos de fondos sugiere 
talar en el año 11 obteniéndose un van de 
$AR 511 914 y utilizando una tasa ajus- 
tada por riesgo {k) de 28,39% 12 siendo 


12 Metodología CAPM ( Capital Asset Pricing 
Model). 


muy superior a las planteadas por Moni- 
cault y Delvalle (2009) y Fassola et al. 
(2004) los cuales no mencionan el criterio 
de selección y adicional por riesgo de 
mercado utilizado en el cálculo. 

El resultado se debió a la alta sensi- 
bilidad negativa de los flujos al factor de 
descuento, reafirmando lo mencionado 
por Soares et al. (2005), además de que la 
variable precio se supuso constante en 
todo el periodo de proyección. 

La aplicación del método de opcio- 
nes reales, por el contrario, sugiere que el 
máximo valor se obtiene en el último 
periodo. Descartada esta opción, que no 
posee flexibilidad, el segundo momento 
óptimo se encontró en el año 18 con un 
ingreso con una vor $AR 1 101 944, 
siendo superior al van en más de 1 00% ya 
que el mismo no permite el desarrollo 
futuro del bosque. 

Estos valores coinciden con lo afir- 
mado por Sant'Anna y Nogueira (2010) y 
Moreira et al. (2000). Además, se observó 
que la componente ganadera en el sis- 
tema no altera la decisión de tala ya que 
en forma independiente (sólo componente 
forestal) o como sistema se obtuvo el 
turno óptimo en el año 18. 


CONCLUSIONES 

Desde el punto de vista económico- finan- 
ciero el enfoque de or presenta importan- 
tes ventajas en contraposición al clásico 
van debido a su robustez, ya que incor- 
pora flexibilidad estratégica y comporta- 
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miento estocástico de los precios con 
base en series históricas; de esta manera 
elimina el uso de la tasa ajustada por 
riesgo. 

Esta cuestión no se tiene en cuenta 
en el van, siendo la tasa ajustada por 
riesgo constante en la vida del proyecto, 
por lo cual supone un escenario constante 
en todo el horizonte de planificación, 
siendo esto un enfoque muy limitado, en 
especial tratándose de proyectos foresta- 
les, dado que presentan horizontes tem- 
porales extensos. 


REFERENCIAS 

Bettinger, P., K. Boston, J.P. Siry y D.L. 
Grebner. 2009. Forest management 
and planning. la. ed. Academic 
Press, California. 360 p. 

Black, F. y M. Scholes. 1973. The pricing 
of options and corporate liabilities. 
Journal of Political Economy 
81 (3):637-659. 

Blasco, C., S. Carenzo y E. Astrada. 2005. 
Evaluación de un sistema silvopasto- 
ril sobre vinales en Formosa, Argen- 
tina. Revista Iberoamericana de 
Economía Ecológica 2(1 ):57-68. 

Brennan, M. y E. Schwartz. 1985. A new 
approach to evaluating natural 
resource investment. Midland Corpo- 
rate Financial Journal 1(3):37-47. 

Cerruto R., C. Leite, J. Gongalves y M. 
López da Silva. 2007. Estimativa do 
abatimento de erosáo aportado por 
um sistema agrossilvipastoril e sua 
contrubugáo económica. Revista 
Árvore 31(2):285-293. 

Chidiak, M., A. Moreyra y C. Greco. 2003. 
Captura de carbono y desarrollo 
forestal sustentable en la Patagonia 


Argentina: Sinergias y Desafíos. 
Documento cenit-cepal-udesa. 136 p. 

DRLhttp://www.territorioscentroamerica- 
nos.org/ecoagricultura/Documents/ 
capturadecarbonoydesarrollofores- 
tal.pdf (último acceso 08/08/2013). 

Colcombet, L., E. Crechi, A. Keller, N. 
Pachas, H. Fassola, S. Lacorte y J. 
Esquivel. 2010. Comparación preli- 
minar de resultados financieros y 
económicos de proyectos forestales, 
ganaderos y silvopastoriles en Misio- 
nes. Acta 14 as Jornadas Técnicas 
Forestales y Ambientales, unam - eea 
Montecarlo, INTA. 86 p. 

Copeland, T., T. Koller y J. Murrin. 2000. 
Valuation: measuring and managing 
the valué of companies. 3 a ed. Wiley. 
Nueva York. 490 p. 

Cox, J., S. Ross y M. Rubinstein. 1979. 
Option pricing: a simplified approach. 
Journal of Financial Economics 
7(3):229-263. 

Crechi, E., R. Fernández, A. Friedl y J. 
Falher. 2005. Efectos de raleo en 
Pin us taeda, Norte Misionero. 
Revista Idia 1 1 (8):41 -45. 

Crechi, E., H. Fassola y A. Friedl. 2005. 
Simulador Forestal. Revista Idia 
11 (8):46-50. 

Damodaran, A. 2002. Investment Valua- 
tion. 2 a ed. Wiley. Nueva York. 992 p. 

Dixit, A. y R. Pindyck. 1994. Investment 
under uncertainty. I a ed. Pricenton 
University Press. Nueva Jersey. 
476 p. 

Fassola, H., S. Lacorte, J. Esquivel, F. 
Colcombet, F. Moscovich, E. Crechi, 
N. Pachas y A. Keller. 2004. Siste- 
mas silvopastoriles en Misiones y ne 


96 


Opciones reales para determinar el turno óptimo en sistemas silvopastoriles 


de Corrientes y su entorno de nego- 
cio. Directorio de información Insti- 
tuto Nacional del Tecnología 
Agropecuaria (inta). 

L/R/_:http://anterior.inta.gob.ar/monte- 
carlo/info/documentos/forestales/ 
sistsilvopastoril.PDF (último acceso 
08/08/2013). 

Grainger, M. 1968. Problems affecting the 
use of Faustmann’s formula as a val- 
uation tool. New Zealand Journal of 
Forestry (1968) 1 3(2):1 68-1 83. 

Hull, J. 2005. Futures, options and other 
derivatives. 6 a ed. Prentice Hall. 
Nueva Jersey. 816 p. 

Kulatilaka, N. 1988. Valuing the flexibility 
of flexible manufacturing systems. 
ieee Transactions in Engineering 
Management 35(4):250-257. 

Malliaris, A. G. y W.A. Brock. 1982. Sto- 
chastic methods in economics and 
finance. I a ed. Holanda del Norte, 
Amsterdam. 317 p. 

Margrabe, W. 1 978. The valué of an option 
to exchange one asset for another. 
Journal of Finance 33(1 ):1 77-1 86. 

Markowitz, H. 1991. Portfolio selection: 
efficient diversification of invest- 
ments. 2 a ed. Blackwell Publishers. 
Cambridge, Massacchusetts. 384 p. 

Masón, S. y R. Merton. 1985. The Role of 
contingent claims analysis in corpo- 
rate finance. División of Research, 
Gradúate School of Business Admin- 
istration, Harvard University. Home- 
-wood Irwin, Nueva York. 65 p. 

Me Donald, R. y D. Siegel. 1981. Option 
pricing when the underlying asset is 


non-stored. Discussion Papers num- 
ber 512. Center for Mathematical 
Studies in Economics and Manage- 
ment Science, Northwestern Univer- 
sity. Evaston, Illinois. 12 p. 

Merton, R. 1973. The theory of rational 
options princing. The Bell Journal of 
Economics and Management Sci- 
ence 4(1 ):1 41 -1 83. 

Miller, R. A. y K. Voltaire. 1 983. A stochas- 
tic analysis of the tree paradigm. 
Journal of Economic Dynamics & 
Control 6(1):371-386. 

Milanesi, G. y H. Vigier. 2010. Árboles de 
decisión, opciones reales y enfoque 
integrado en mercados completos e 
incompletos. Anales Asociación 
Argentina de Economía Política. 
Buenos Aires. p:1-30. 

Milanesi, G. 2011 . Opciones reales: teoría 
y casos prácticos. I a ed. Induvio, 
Buenos Aires. 156 p. 

Milanesi, G., G. Woitschach y D. Broz. 
(2012). Aplicación de la teoría de 
opciones reales a la determinación 
del momento óptimo de cosecha 
forestal. Rev. fca uncuyo. 44(2) :65- 
78. 

Monicault, L. A. y P. Delvalle. 2009. Eva- 
luación Económico-Financiero de un 
Sistema Silvopastoril (ssp). Docu- 
meto técnico inta online. URLhttp:// 
inta.gob.ar/documentos/evaluacion- 
economico2013financiero-de-un-sis- 
tema-silvopastoril-ssp/ (último 
acceso 08/08/2013). 

Moreira A., E. Reis, K. Rocha y L. Carva- 
Iho. 2000. A valoragáo das conces- 
sóes ñas florestas nacionais da 
Amazonia: urna abordagem com 


Madera y Bosques 19(3), 2013:81-98 


97 


opgóes reais. Pesquisa e Planea- 
mento Económico 30(3):327-354. 

Mun, J. 2005. Real options analysis: tools 
and techniques for valuing strategic 
investment and decisions. 2 a ed. Wiley. 
Nueva York. 704 p. 

Musálem, M.A. 2002. Sistemas agrosivo- 
pastoriles: una alternativa de desarro- 
llo rural sustemtable para el trópico 
mexicano. Revista Chapingo 8(2):91- 
100 . 

Myers, S. y S. Majd. 1990. Abandonment 
valué and project lite. Advances in 
Futures and Options Research 4:1-21. 

Myers, S. 1977. Determinants of corporate 
borrowing. Journal of Financial Eco- 
nomics 5(2):147-175. 

Paddock, J., D. Siegel y J. Smith. 1988. 
Option valuation of claims on physical 
assets: the case of Offshore Petro- 
leum Lease. Quaterly Journal of Eco- 
nomics 103(3):479-508. 

Pindyck, R. 1988. Irreversible investment, 
capacity choice and the valué of the 
firm. American Economic Review 
78(5):969-985. 

Plantinga, A. 1 998. The optimal timber rota- 
tion: an option valué approach. Forest 
Science, Society of American Forest- 
ers 44(2): 192-202. 

Pratt, S. y R. Grabowski. 2008. Cost of cap- 
ital: applications and examples. 3 a ed. 
John Wiley & Sons Nueva Jersey. 
778 p. 

Prodan, M., R. Peters, F. Cox y P. Real. 
1997. Mensura Forestal. I a ed. Costa 
Rica. Instituto Interamericano para la 
Agricultura - Gesellschaft für Zusam- 
menarbeit (iica-gtz). 511 p. 


Rocabert, J. (2007). Los criterios valor 
actual neto y tasa interna de rendi- 
miento. Revista electrónica sobre la 
enseñanza de la Economía Pública 
2 ( 1 ): 1 - 1 1 . 

Sant'Anna, A. y J. Nogueira. 2010. Valora- 
gao económica dos servigos abientais 
de florestas nacionais. Revista de 
Administragáo e Negocios da Amazo- 
nia 2(1):82-108. 

Schneider, P. 2002. Manejo florestal: plane- 
jamento da produgáo florestal. Santa 
María, Universidade Federal de Santa 
María, Centro de Ciencias Rurais, 
Departamento de Ciencias Florestais. 
492 p. 

Smit, H. 1997. Investment analysis of off- 
shore concessions in The Netherlands. 
Financial Management 26(2):5-17. 

Smith, J. 2005. Alternative approach for 
solving real options problems. Deci- 
sión Analysis 2(2):89-102. 

Soares, T.S., R. M. Carvalho,A. Bartolomeu 
do Vale. 2003. Avaligáo económica de 
um povoamento de Eucalyptus gran- 
dis destinado a multiprodutos. Revista 
Árvore. 27(5): 689-694. 

Thomson, T.A. 1992. Optimal forest rotation 
when stumpage prices follow a difu- 
sión process. Land Economics 
68(3):329-342. 

Trigeorgis, L. y S. Masón. 1987. Valuing 
Managerial Flexibiliy. Midland Corpo- 
rate Finance Journal 5(1 ):1 4-21 . 

Trigeorgis, L. 1988. A conceptual options 
framework for capital budgeting. 
Advances in Futures and Options 
Research 4(3): 145-1 67. 


98 


Opciones reales para determinar el turno óptimo en sistemas silvopastoriles 


Trigeorgis, L. 1997. Real options: mana- 
gerial flexibility and strategy in 
resource allocations. 2 a ed. MIT 
Press, Cambridge. 427 p. 


Yin, R. y D. Newman. 1997. When to cut a 
stand of trees? Natural Resource 
Modeling 10(3):251-261. 


Manuscrito recibido el 8 de abril de 2012. 

Aceptado el 8 de agosto de 2013. 

Este documento se debe citar como: 

Milanesi, G., D.R. Broz y G.B. Woitschach. 2013. Opciones reales para determinar el turno óptimo en sistemas sil- 
vopastoriles: análisis de inversión. Madera y Bosques 1 9(3):81 -98. 


MADERA Y BOSQUES Vol. 19 Núm. 3 

Se terminó de imprimir en el mes de diciembre de 2013 
en los talleres de Tavera Hermanos, S.A. de C.V. 

Av. Lázaro Cárdenas 3052. Col. Chapultepec Sur. 
CP 58260, Morelia, Michoacán. 

La edición consta de 300 ejemplares 
más sobrantes para reposición. 



Madera y Bosques 19(3), 2013 
CONTENIDO 


3 Editorial 

Artículos de investigación 

7 Aprovechamiento del recurso maderable tropical nativo en la 
comunidad de Emilio Rabasa, Reserva de la Biosfera Selva 
El Ocote, Chiapas, México 

Carolina Orantes-García, Miguel Ángel Pérez-Farrera, 

Carlos Uriel del Carpio-Penagos y César Tejeda-Cruz 

23 Respuesta de Cedrela odorata L. a diversos inoculantes 
micorrízicos procedentes de dos ecosistemas tropicales 
Heriberto Méndez-Cortés, José G. Marmolejo-Monsiváis, 
César Cantú-Ayala, Víctor Olalde-Portugal, 

Eduardo Estrada-Castillón y César Posadas-Leal 

35 Influencia de variables físicas en la productividad de Pinus 
arizonica y Pinus engelmannii en el sur de Chihuahua, México 
Martín Martínez-Salvador, Ricardo David Valdez-Cepeda 
y Marín Pompa García 

51 Efecto de la micorrización en el establecimiento de Pinus 
pseudostrobus en cárcavas 
Mariela Gómez-Romero, Javier Villegas, 

Cuauhtémoc Sáenz-Romero y Roberto Lindig-Cisneros 

65 Tableros de partículas fabricados con residuos industriales de 
madera de Pinus patula 
Luis Zambrano, Pablo Moreno, Freddy Muñoz, 

Jorge Durán, Darío Garay y Styles Valero 

81 Opciones reales para determinar el turno óptimo en sistemas 
silvopastoriles: análisis de inversión 
Gastón Milanesi, Diego Ricardo-Broz 
y Guillermo Bruno Woitschachz 


Guía de autores 


